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INTRODUCTION 


J.  V.  Foa 

This  report  is  a  summary  of  work  done  at  Rensselaer  since 
1961,  as  part  of  a  program  of  research  on  the  dynamics  of  a  body 
njoving  in  a  tube,  under  a  three-year  grant  from  the  U.  S.  Army 
Research  Office  (Durham).  More  specifically,  this  report  covers 
the  portion  of  the  program  that  dealt  with  problems  relating  to 
a  transportation  system  in  which  a  vehicle  would  propel  itself  at 
high  velocity  (250  to  2000  mph)  through  a  tube. 

Other  research  activities  under  the  same  grant  have  been  re¬ 
ported  separately  (Refs.  1-3). 

THE  CONCEPT 

Land  vehicles  are  of  necessity  constrained  to  travel  along 
narrowly  defined  paths.  Therefore,  unless  the  prescribed  path  is 
a  straight  line,  a  land  vehicle  traveling  at  very  high  speed  must 
be  supported  from  all  sides.  This  kind  of  support  is  best  provided 
within  a  tube,  also  because  there  are  safety  features  that  are  in¬ 
trinsic  in  the  operation  of  a  high-speed  vehicle  within  a  tube. 

However,  in  conventional  tube  transportation  systems,  the 
vehicle,  in  traveling  within  the  tube,  acts  in  relation  to  the  ad¬ 
jacent  air  masses  either  as  a  driving  piston  (as  in  the  case  of  a 


train  in  a  tunnel)  or  as  a  driven  pisjtoh  (as  in  pneumatic  dis¬ 
patch  systems),  and  therefore'  its  ipotidn  is  aecoinpanied  by  large 
displacements  of  the  air  masses  ih.“t'he  tube.  Irt  contrast  to  the 

incompressible-flow  situation  (Ref.  4),  where  these  displacements! 

•< 

decrease  the  drag  of  the  vehicle,  their  effect  dn  a  compressible 
medium  is  found  to  be  an  increase  of  drag.  The  increase  is  so 
large  that  the  required  thrust  could  not  be  produced  by  wheel 
traction  at  speeds  much  above  those  of  present-day  land  transporta¬ 
tion  even  if  there  were  no  structural  limitations  to  the  use  of 
this  kind  of  propulsion  at  high  speeds.  Propulsion  by  linear  in¬ 
duction  could  extend  the  range  of  attainable  speeds  somewhat,  but 
only  at  a  prohibitive  cost  in  power  consumption.  To  overcome  thii* 
difficulty,  schemes  have  been  proposed  for  the  operation  of  vehici 
in  evacuated  tubes,  but  any  systems  of  this  sort  would  be  very  cos 
ly  to  operate  (despite  its  low  power  demands)  apd  highly  vulnerablt 
Furthermore,  at  speeds  at  which  the  lyte  of  wheels,  levapads,  or 

,  (p 

other  direct -contact  or  low -clearance  supports  must  be  ruled  out, 
the  evacuation  of  the  tube  aggravates  the  problem  of  supporting  th 

■p; 

vehicle  and  leaves  a  magnetic  suspension,  with  Its  own  pcicullar 

a 

difficulties  and  hasards,  aa  the  only' workable  solution. 

p 

9 

A  study  of  the  bropd  problem  of  high-speed  land  transportation 
with  consideration  of  the  intertwined  difficulties  that  this  probl 
Involvaa,  lad,  in  1947,  to  the  conception  of  the  scheme  that  Is 
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discussed  in  the  present  report.  The  work  that  was  done  on  this 
scheme  at  Rensselaer,  as  an  internal  research  project,  prior  to 
the  initiation  of  the  studies  that  are  summarized  here,  is  des¬ 
cribed  in  Refs.  5-7. 

The  system  under  discussion  overcomes  the  speed  limitations 
of  conventional  land-transport  systems  without  requiring  an  evacua¬ 
ted  tube.  Its  main  feature  will  now  be  described  in  some  detail. 

PROPULSION 

The  vehicle  propels  itself  by  the  fore-to-aft  transfer, 
through  the  space  which  separates  it  from  the  tube  wall,  of  the 
air  which  fills  the  remainder  of  the  tube.  The  transfer  is  the 
process  by  which  the  motion  of  the  vehicle  is  produced  and  main¬ 
tained,  not  the  consequence  of  a  motion  that  is  produced  by  other 
means.  Its  volumetric  rate  is  such  that  the  flow  perturbations 
in  the  tube  are  kept  at  the  lowest  level  that  is  compatible  with 
the  maintenance  of  the  required  motion  of  the  vehicle.  This  con¬ 
dition  would  be  ideally  and  automatically  satisfied  if  the  fluid 
filling  the  tube  were  incompressible  and  if,  in  addition,  its 
motion  were  fully  prevented  at  some  station  by  means  of  a  gate  or 
simply  by  keeping  one  end  of  the  tube  closed.  In  this  case,  the 
volumetric  transfer  rate  would  always  be  equal  to  the  product  of 
the  cross-sectional  area  of  the  tube  and  the  velocity  of  the 
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vehicle,  and  the  fluid  would  be  at  rest  everywhere  in  the  tube 
except  in  the  transfer  passages,  regardless  of  the  mode  of  pro¬ 
pulsion.  Because  of  the  compressibility  of  the  air,  the  flow 
perturbations  produced  by  the  pressure  disturbances  which  origi¬ 
nate  at  the  vehicle  may  extend  over  great  distances,  "owever, 
in  contrast  to  the  incon^jressible-flow  situation,  the  amplitude 
and  distribution  of  these  perturbations  depend  not  only  on  the 
boundary  conditions  but  also  on  the  shape  of  the  vehicle,  on  its 
blockage  ratio,  on  its  speed,  and  on  its  mode  of  propulsion.  In 
particular,  it  has  been  found  that  modes  of  propulsion  of  the  tj^ne 
specified  above  ("internal”  modes)  are  capable  of  producing  much 
Sioaller  flow  perturbations,  all  other  conditions  being  equal,  than 
the  conventional  ("external")  modes,  in  which  the  thrust  is  pro¬ 
duced  as  the  reaction  to  a  force  exerted  on  a  stationary  structure 
With  internal  propulsion,  if  the  tube  is  relatively  short  and  the 
motion  of  the  fluid  is  inhibited  at  any  one  station,  the  transient 
flow  velocities  that  are  generated  in  the  tube  even  during  acceler. 
tion  and  deceleration  of  the  vehicle  are  quite  low  (Section  A) .  Th 
analysis  also  predicts  remarkiibly  low  flow  velocities  with  the 
present  system,  during  accelerations  of  the  vehicle,  for  the  ideal 
zed  case  of  a  tube  of  Infinite  length  (Section  A).  This  means  tha 
no  gate  or  other  flow-inhibiting  device  is  necessary  when  the  ac¬ 
celerations  start  at  such  a  distance  from  the  nearest  opening  in 
the  tube  that  no  reflected  disturbance  can  reach  the  vehicle  befor 
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the  acceleration  phase  is  completed.  In  the  steady  state  that  is 
finally  established  by  the  reverberating  pressure  waves, the  flow 
velocities  are  higher,  but  still  very  much  lower  than  with  external 
propulsion  (Section  F). 

Since  the  slipstream  of  the  propulsion  system  is  here  sub¬ 
stantially  at  rest,  its  "propulsive  efficient’  "  is  close  to  unity. 
Furthermore,  the  wake  associated  with  the  dra  '  of  the  vehicle  is 
also  at  rest  in  a  stationary  frame  of  reference.  Thus  the  kinetic 
energy  losses  associated  with  fluid  motions  in  both  the  &>lipstream 
and  the  wake  are  almost  entirely  eliminated.  The  power  required 
for  propulsion  of  the  vehicle  at  ataady  speed  is^  accordingly, 
greatly  reduced.  At  supersonic  speed,  this  power  is  reduced  to 
that  necessary  to  make  up  for  flow  losses  inside  the  vehicle  and 
in  the  annular  space  surrounding  it  within  the  tube. 

In  fact,  at  supersonic  speeds,  Che  vehicle  itself  may  be  con¬ 
sidered  as  the  equivalent  of  the  centerbody  of  a  jet  engine,  the 
outer  shroud  of  which  is  represented  by  the  wall  of  Che  tube.  The 
power  required  at  steady  speed  can  then  be  calculated  by  analogy 
with  that  of  a  jet  engine  in  the  condition  of  zero  thrust.  Fig.  1 
Illustrates  this  equivalence  for  a  system  in  which  the  transfer 
flow  is  sustained  by  the  addition  of  heat  alone,  as  in  a  ramjet. 

It  should  be  noted  that  the  outer  boundary  of  the  airstream  passage 
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(the  internal  surface  of  the  tube  wall)  is  not  stationary  in 
the  vehicle-fixed  frame  of  reference  as  it  is  in  conventional 
jet  engines i  instead,  it  moves  in  this  frame  of  reference  in  the 
direction  of  the  relative  air  flow.  Energy  losses  due  to  viscous 
forces  are  therefore  substantially  reduced. 

Problems  relating  to  drag  and  propulsion  are  discussed  in 
Sections  A,  E,  F,  G,  H  And  J  of  this  report. 

SUPPORT 

It  has  already  been  noted  that  a  land  vehicle  traveling  at 
very  high  speed  must  be  supported  on  all  sides.  In  the  lower 
range  of  speeds  the  required  support  may  be  provided  by  wheels 
or  skids  or  low-clearance  devices,  such  as  the  "levapads".  Large- 
clearance  suspensions  are  preferable  at  high  speeds.  Over  the 
entire  range  of  speeds,  the  vehicle  may  be  kept  on  a  centered 
path  through  the  use  of  suitably  arranged  GEM  pads.  Since  the 
air  curtains  in  Chase  supports  would  be  blowing  against  the  solid 
walls  of  the  tube,  their  "rigidity”  would  not  be  limited  as  it  is 
in  conventional  ground  effect  machines  that  are  to  operate  over 
soft  ground  or  water.  Therefore,  this  kind  of  support  could  be 
made  to  operate  effectively  within  the  tube  at  speeds  far  above 
the  maximum  permissible  with  "conventional"  ground  effect  machines. 
In  forward  motion,  the  rearward  momentum  of  the  deflected  curtains 
contributes  Co  the  thrust  (Section  B).  Floating  suspensions  of 
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this  type  also  have  the  important  merit  that  a  vehicle  so  pro¬ 
vided,  if  its  center  of  gravity  is  appropriately  low,  will 
automatically  tilt  itself  to  the  correct  angle  of  bank  in  every 
turn,  independently  of  the  speed  at  which  the  turn  is  negotiated. 

Problems  of  stability  and  ground-effect  support  of  the 
vehicle  are  discussed  in  Sections  B,  C,  and  D  of  this  report. 

TR^',>'oONIC  ACCELERATION 

In  dealing  with  supersonic  vehicles,  a  variety  of  schedules 
of  acceleration  through  the  transonic  regime  may  be  considered. 
Some  aspects  of  two  of  these  schedules  are  discussed  in  Section  I. 
According  to  the  first  one,  the  transfer  flow  is  permitted  to 
choke,  with  the  formation  of  a  fbrwar dr- facing  shock  ahead  of  the 
vehicle;  and  the  vehicle  is  then  made  to  overtake  this  shock.  The 
second  sr.hedule  calls  for  surge  chambers  around  the  tube  walls, 
or  for  wall  perforations,  oyer  those  stretches  of  the  tube  within 
which  the  transonic  acceleration  Is  e;>cpec ted  to  take  place.  Others 

«  *  ■'  l.*  ' 

yet  are  based  on  the  utlliMtien  of  pi^oeesaes  in  the 

6  .  • 

transfer  flow.  '  •  .  : 


BRAKING  ;  :  ;•  _  ' 

♦  *  •  * 

The  presence  of  large  and  confined  masses  of  ,air  within  the 

•  ,  . 

tube  provides  a  safe  and  simple  solution  to  the  problem  of.  bringing 
the  vehicle  to  rest  from  speeds  at  which  it  wojild  not  be  practicable 
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to  dissipate  all  the  stored  kinetic  energy  through  sliding 
friction  between  solid  surfaces,.  As  soon  as  the  input  of  pro¬ 
pulsive  power  to  the  transfer  flow  is  shut  off,  the  behavior  of 
the  vehicle  in  relation  to  the  adjacent  fluid  masses  changes  to 
the  already  discussed  pistcnilke  action  of  conventionally- driven 
vehicles.  The  result  is  a  large  retarding  force,  w'hich  may  be 
increased  further  through  the  Insertion  of  spoilers  in  the  trans¬ 
fer  passage  or  may  be  decreased  through  the  opening  of  appropriate 
vents.  It  is  easy  to  verify  that  this  mechanism  is  capable  of 
producing,  at  high  forward  speeds,  any  braking  action  that  may  be 
desired  or  required  even  under  sevete  circumsstances .  As  the  vehi¬ 
cle  slows  down  this  braking  force  decreases,  until  the  velocity  cf 
the  vehicle  becomes  low  enough  that  it  can  be  brought  to  a  full 
stop  by  conventional  sliding-friction  means.  In  any  case,  the 
dissipation  of  ttie  kinetic  energy  initially  stored  in  the  vehicle 
does  not  entail  any  difficulty  with  local  heating  or  wear,  because 
the  mechanisni  just  described  causes  this  energy  to  become  distri¬ 
buted  over  long  stretches  of  the  tube. 

ILLUSTRATIVE  EMBODIMENTS 

Illustrative  embodiments  of  the  concept  are  shown  in  Figures 
2  through  6.  Fig.  2  shows  a  low-speed,  wheel-supported  vehicle 
with  a  propeller -energised  flow.  In  the  arrangement  of  Fig.  3  a 


compressor  energizes  the  flow  both  for  propulsion  and  for  ground- 
effect  support. 

Fig.  4  shows  a  vehicle  with  a  driving  mechanism  which  is 
essentially  an  ejector  pump.  The  primary  gas  is  discharged  from 
a  con5)ressor  or  other  suitable  gas  generator  into  a  mixing  duct 
behind  the  vehicle,  where  the  ejector  action  takes  place.  The 
mixing  duct  is  formed  by  the  tube  itself,  as  shown  in  the  lower 
portion  of  the  figure.  There  is  a  similarity  between  this  embodi¬ 
ment  and  that  of  Fig.  1,  in  that  the  transfer  flow  is  energized, 
in  both  cases,  in  a  region  where  the  addition  of  energy  to  the  flow 
would  have  no  useful  effect  if  the  tube  v?ere  absent.  The  tube  con¬ 
stitutes,  in  both  cases,  an  integral  part  of  the  propulsion  system. 
For  example,  when  the  primary  gas  generator  is  a  turbojet  (as  in 
the  embodiment  of  Fig.  4) ,  the  thrust  generator  may  be  said  to  be 
a  shrouded  turbojet  in  which  the  shroud  Is  replaced  by  a  stationary 
tube. 


In  the  embodiment  of  Fig.  5  the  driving  mechanism  is  again,  an 
ejector  pump,  but  the  transfer  flow  is  energized  upstream  of  the 
throat. 

Fig.  6  shows  a  supersonic  arrangement  in  which  the  transfer 
flow  is  energized  by  a  cryptosteady  energy  exchanger  (Ref.  8)  ahead 
3f  the  throat  and  by  heat  addition  downstream  of  the-  throat.  In  a 
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variation  of  this  embodiment,  the  addition  of  heat  is  effected 
behind  the  vehicle  by  sources  which  are  fixed  to  the  tube  and 
successively  activated  by  the  passing  vehicle  (Section  F). 
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"The  Flow  Induced  by  a  New  Type 
of  Self-Propelled  Vehicle  In  a  Duct" 


SUIffARY 


An  analysis  is  carried  out  to  determine  the  magnitude  of 
the  flow  velocities  that  are  Induced  Inside  a  tube  In  the  presence 
of  a  vehicle  which  propels  Itself  by  the  fore-to-aft  transfer,  at 
appropriate  rates,  of  the  fluid  In  the  adjacent  layers  within  the 
tube  itself.  A  fixed  schedule  of  accelerations  of  the  vehicle  is 
considered,  with  various  tube  end  conditions. 

The  greatest  induced  air  velocities  are  produced  if  the  tube 
is  open  at  both  ends.  If  the  tube  is  closed  at  one  or  both  ends 
the  induced  air  velocities  are  more  than  an  order  of  magnitude 
lower.  In  an  infinite  or  half -infinite  tube  the  induced  air  velo¬ 
cities  are  greater  than  those  produced  in  a  closed  tube,  but  still 
an  order  of  magnitude  lower  than  those  produced  in  a  tube  open  at 
both  ends. 
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2.V  refers  to  conditions  adjacent  to  rear  face  of  vehicle 

Superscripts 

o  stagnation  conditions 
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INTOODUCTION 

In  recent  years  a  number  of  reports  have  been  written 
(Refs,  1,  2,  3,  and  4)  concerning  a  method  of  high-speed  land 
transportation  involving  a  self-propelled  vehicle  operating  in 
a  tube.  This  concept  of  transportation,  first  proposed  by  Foa 
at  Cornell  Aeronautical  Laboratory  in  1947,  consists  of  a  vehicle 
propelling  itself  through  a  tube  by  the  fore-to-aft  transfer, 
through  or  around  itself,  of  the  air  which  fills  the  remainder 
of  the  tube,  This  is  in  contrast  to  conventional  systems  of 
transport  in  tubes,  or  of  pneumatic  dispatch,  where  the  column  of 
air  filling  the  entire  length  of  the  tube  (or  twice  that  length 
in  the  case  of  a  closed-loop  system)  mast  be  moved  abo*'t  the 
same  distance  as  the  vehicle,  at  about  the  same  velocity. 

The  proposed  system  en^jloys  the  air  in  the  tube  to  provide 
both  thrust  for  propulsion  and  drag  for  braking,  and  its  operation 
is  accompanied  by  disturbances  in  the  air  columns  adjacent  to  the 
vehicle.  The  present  report  is  concerned  only  with  the  accelera¬ 
tion  phase  and  more  particularly  with  a  determination  of  tne  order 
of  magnitude  of  the  induced  air  velocities  which  are  produced  by 
accelerating  vehicles  of  the  type  described  above,  operating  in 
tubes  with  various  end  conditions  (open,  closed,  and  half -open 
tubes,  finite,  infinite,  and  half-infinite  tubes). 
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ANALYSIS 


Statement  of  the  Problem 

A  rigorous  study  of  the  dynamical  situation  considered  here 
would  require  a  detailed  analysis  of  the  fluid  flow  around  or 
through  the  vehicle.  Great  simplifications  are  possible  if  vlscor^ 
stresses  are  assumed  to  be  negligible,  with  the  understanding  that 
the  validity  of  this  assumption  is  to  be  checked  in  the  light  of 
the  results  (end  particularly  of  the  flow  velocities)  so  obtained. 
Then,  the  flow  may  be  treated  as  one-dimensional  everywhere,  except 
for  a  limited  region  coaprlsing  the  fluid-transfer  passages  and 
their  inniediate  neighborhood.  Furthermore,  if  the  duct  is  suffi¬ 
ciently  long,  this  region  may  be  treated  as  a  discontinuous  "actuate 
disk",  and  the  flow  may  accordingly  be  treated  as  one -dimensional 
everywhere.  A  schematic  representation  cf  this  analytical  model 
is  shown  in  Fig,  1  with  the  vehicle  moving  to  the  right  and  the  flow 
relative  to  the  vehicle  moving  to  the  left.  All  velocities  and  dis¬ 
tances  relative  to  the  duct  ere  considered  positive  to  the  right. 

In  the  present  ansiysis  only  subsonic  vehicle  speeds  are 
considered.  To  obtain  a  connon  basis  for  comparison,  the  stipulated 
schedule  of  accelerations  of  tha  vahicle  is  the  same  for  all  con¬ 
ditions  considered. 

Vehicle  Ferformance 


The  thrust  acting  on  the  vehicle  is  equal  to  the  change  of 
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stream  force  (relative  to  the  vehicle)  of  the  flow  as  it  passes 
through  the  vehicle,  or 

0  -  Xljv)  -  (■■Ui.-Uv)}  +  wfllv  -  P„)  (1) 

v<fhere  9  is  the  thrust,  m  is  the  mass  flow  through  the  vehicle, 
relative  to  the  vehicle,  v^y  and  are  the  velocity  and 

pressure,  respectively,  of  the  flow  just  behind  the  vehicle,  and 
V,y  and  p)^  are  the  velocity  and  pressure,  respectively,  of 
the  flow  just  ahead  of  the  vehicle.  Since  the  velocity  change  a- 
cross  the  vehicle  is  small,  the  contribution  of  the  momentum  change 
to  the  thrust  is  small  and  will  consequently  be  neglected  in  the 
present  analysis.*  Let  the  pressure  rise  across  the  vehicle  be 

fi.  -  P.  -■  ‘P.(l  -  (2) 

where  is  a  reference  pressure,  Vy  the  vehicle  velocity,  Ut 

the  terminal  velocity  of  the  vehicle,  and  s  a  constant  less  than 
1.0.  Then 

e  "  -  P,)  ■  TT),  (3) 


The  absolute  value  of  the  thrust  contributed  by  the  momentum 
change  is  sm^,ll  for  all  cases  considered  in  the  Numerical 
l!xanxple,  but  the  percentage  of  the  total  thrust  which  the 
i33omentUQi  change  contributes  Increases  from  less  than  107.  when 
the  vehicle  first  starts,  to  207.  or  30%  at  higher  vehicle 
velocities  when  the  pressure  difference  drops  to  a  small  value. 
Therefore  the  approximation  of  the  present  analysis  is  poor  in 
the  latter  range  of  speeds,  which  is,  however,  the  range  of 
least  interest,  in  that  it  involves  the  lowest  flow  velocities. 


This  relation  Is  plotted  In  Fig.  2. 
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Air  Velocities  Induced  by  the  Vehicle 

The  air  velocities  which  are  induced  in  the  tube  by  the  accelera¬ 
ting  vehicle  depend  on  the  tube  end  conditions  on  both  sides  of  the 
vehicle.  The  end  conditions  considered  here  are:  (1)  fully  open  end, 
(2)  fully  closed  end,  and  (3)  no  end  (tube  extending  to  infinity). 

Two  of  the  nine  possible  combinations  of  these  conditions  at  the  two 
erids  admit  no  solution  for  the  analytical  model  used.  These  two 
cases  are  those  of  a  half -Infinite  tube  on  either  side  of  the  vehicle 
with  the  tube  end  closed  on  the  other  side.  They  require  a  gradual 
increase  of  the  pressure  increment  from  zero  to  the  value  given  by 
Eq,  2,  whereas  the  analytical  model  used  here  requires  a  finite  pres¬ 
sure  rise  across  the  vehicle  from  the  start. 

Each  of  the  other  seven  cases  must  be  solved  separately,  as  they 
require  different  assumptions  and  simplifications. 

1)  Both  ends  closed 

For  this  tube  configuration,  the  flow  processes  are  treated  as 
quasi -static.  The  air  temperature  is  stipulated  to  be  constant  with 
respect  to  time  and  uniform  throughout  the  tube,  so  that  the  density 
i“  dependent  only  on  the  pressure.  The  pressure  is  considered  constant 
throughout  each  of  the  two  regions  ahead  of  and  behind  the  vehicle 
(  «  P  *  P  ^  *  Pjy  *  Pt  total  mass  of 

air  within  the  tube  is  constant 

^  yL 


+  R(l-x,)  •  P.L 
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where  L  is  the  total  length  of  the  tube.  The  pressure  ahead  of 
vehicle  (6n  Its  right  side  in  the  figures)  is  then 

f?  =  B  -  (p.-r)  r 


R  «  P.{i  - 


and  from  Eq,  (2) 


P,  -  R  (  I  +  <1  -  WJ  -  l)} 


(10) 


Consider  the  layer  of  air  which  at  tine  t  is  at  a  distance 
y  from  the  left  end,  in  region  1.  The  volfmne  and  mass  of  air 
closed  between  this  layer  and  the  right  end  of  the  tube  are  (L  -• 
and  ,  respectively.  The  latter  is  constant,  whe 

the  former  chenges  with  time.  Denoting  the  velocity  of  this  laye 
by  Vm  , 


i3£:  ,  -^-^1 


(i-’in 


'Um 

(L-x)R 


For  «n  IdMl 


iv  .  RZ  dR 

dt  Ft^  dt 


Thus, 


■u,.  =  (l-<)  f, 


A 
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j  p 

Substituting  tor  R  and  ^  from  Eq,  9  yields; 


This  equation  gives  the  velocity  distribution  in  region  1  at  the 
time  when  the  instantaneous  velocity  of  the  accelerating  vehicle  is 
Uv 


Similarly,  in  region  2,  denoting  the  instantaneous  flow  velo¬ 
city  at  X  by  V-ix 


and  vttn  are  linear  functions  of  x  ,  zero  at  the  closed  ends  and 
maximum  at  the  vehicle.  Thus,  at  any  instant  the  maximum  induced 
air  velocities  are  obtained  froa  Eqs,  (11)  and  (12)  evaluated  at 

/  s  Xy 

2)  One  end  open  and  the  other  closed. 

When  the  tube  is  open  at  one  end  the  mass  of  air  in  the  tube 
does  not  remain  constant  during  the  acceleration  of  the  vehicle. 
However,  as  will  be  verified,  the  induced  flow  velocities  are  small 
in  this  situation  also.  Therefore  the  pressure  in  the  region  of  the 
tube  between  the  open  end  and  the  vehicle  remains  close  to  atmospheric 
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pressure,  and  the  flow  In  this  region  can  be  treated  as  tnconpr 
Ible.  If  the  pressure  Is  also  considered  constant  throughout  t 
region  between  the  closed  end  and  the  vehicle,  then  at  each  Ins 
the  pressure  is  uniform  throughout  both  regions,  or  : 

and  tei^erature  is  stipulated  to  be 

stant  with  time  and  uniform  throughout  the  tube,  then  the  flow 
velocity  may  be  considered  uniform  throughout  the  region  betwf.^ 
open  end  and  the  vehicle.  The  continuity  equation  across  the  \ 
gives 

w)  *  (1. 

For  the  case  In  which  the  closed  end  la  behind  the  vehicle 
(region  2)  the  flow  velocity  In  region  1  Is 

Pi  ' 

where  the  preaeure  ratio,  as  given  by  Eq.  (2),  Is 


The  air  velocity  In  the  closed  region  behind  the  vehicle  (  'U 
can  be  found  using  the  sasM  approach  as  In  case  1. 

« 


-J.  dh 

H  dt 
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-  fl  - 

he  {  Ut/ 

' "  ^  -  s:) 


(16) 


The  velocity  is  at  each  instant  a  linear  function  of  X  , 

zero  at  the  closed  back  end  and  maximum  at  the  vehicle.  Thus, 


.  Ui-  .  M  ^  ^  ] 

Vx*  ■  V,  ■  I  "tj  I  ,  ^  f  J 


(17) 


which  shows  that  the  maximum  induced  air  velocity  in  region  2  is 
dependent  only  upon  6  and  ^  . 


Combining  Eqa.  14,  15  and  17  yields 

s  ■  •  2  a) 


(18) 


Thus,  the  air  velocity  in  region  1  la  dependent  also  upon  6  and 


For  the  case  in  which  the  closed  end  is  ahead  of  the  vehicle 


(Region  1)  the  flow  velocity  in  the  Region  2  is 


where  again  the  pressure  ratio  is  found  from  Eq.  (2), 


(19) 


I  - 


s) 


f 
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The  air  velocity  in  Region  1  can  again  be  found  by  using  the  same 
approach  as  in  case  1, 


dr 


or 


W 


(20) 


Thus,  14, »  is  at  each  instant  a  linear  function  of  X  ,  maximur 

at  the  vehicle  and  zero  at  the  front  end.  By  substituting  the 
of  into  Eq,  (19),  the  velocity  in  the  open-end  region  (Reg-; 

2)  is  obtained  as: 


(21) 


It  can  be  verified  that  the  air  velocity  in  Region  2  is  maximum  whe 
‘Uy  «  0  ,  which  occurs  at  Xy  »  O  : 


H  L 

be  ■  n 


From  this  expression  it  can  be  seen  that  the  initial  flow-velocity 
a  function  of  tube  length  (  L  )  and  that,  for  very  long  tubes,  ver 
high  air  velocities  will  result.  This  result  is  a  consequence  of  t 
quasi-static  approach  used.  The  pressure  in  the  open-end  region  is 
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atmospheric  and  because  of  the  pressure  ratio  across  the  vehicle 
(which  is  maximum  when  s  o  )  the  closed-end  region  ahead  of  the 
vehicle  is  under  low  pressure.  As  the  vehicle  begins  to  nx)ve,  the 
pressure  (which  is  assumed  constant  throughout  the  length  of  the 
tube)  begins  to  rise.  In  order  to  satisfy  the  conditions  of  the 
quasi-static  approach,  the  inflow  velocity  must  be  very  high  to 
supply  the  large  amount  of  air  necessary  to  raise  the  pressure  in 
the  entire  tube.  The  in^lausibility  of  the  result  indicates  that 
the  present  analysis  for  this  condition  (tube  closed  ahead  of  the 
vehicle)  is  applicable  only  to  very  short  tubes, 

3)  Infinitely  long  tube 

The  effect  of  the  con^jression  and  expansion  waves  which  the 
vehicle  generates  as  it  accelerates  is  determined  in  the  following 
nonsteady-flow  analysis.  To  avoid  the  complicated  interactions  of 
outgoing  waves  with  reflected  waves,  an  infinite  tube  will  be  con¬ 
sidered.  In  this  case,  if  no  entropy  producing  disturbances  exist 
in  the  tube,  the  value  of  the  Riemann  parameters  (see  Ref.  5)  on 
the  characteristics  which  reach  the  vehicle  from  both  sldet.  remains 


p. 

--  A 

r-i 

*  Uu 

It 

A  iw 

f  ir„ 

It 

(23) 

Q, 

=  ^  A 

t-i 

- 

-  ^ 

A»v 

-  Uv 

(24) 

constant: 
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Since  the  flow  Is  assumed  to  be  hooentroplc,.  Eq«  2  yields 


Substituting  the  above  expressiohs  for  Uy  ,  U/v  and  17, y  int 

the  equation  for  the  conservation  of  mass  across  the  vehicle, 

(u,  - 14)  •  I;  f  c/.  -  u»v)  (26) 

and  using  the  isentroplc  relationship  between  (p  and  T  ,  yields 


This  equation,  which  can  be  solved  numerically,  gives  a  relation^ 
between  Aiy  and  Ay  .  Once  this  relationship  Is  established, 
the  corresponding  values  of  V,v  ,  U»y  and  lAy  can  be  deter 
allied  from  Eqs«  (23),  (24)  and  (25). 

This  solution  gives  the  waves  which  are  generated  as  the  vet 
accelerates  along  the  tube*  The  only  transient  this  solution  omi 
is  the  one  generated  when  the  power  is  first  applied  to  the  vehic 
which  is  at  rest  in  the  tube.  It  does,  however,  give,  under  the 
conditio:^  stipulated,  tfie  velocitv  of  the  flow  after  the  power  i 
been  applied. 


4)  Sfenl-inflnlte  tube 


As  nentioned  previously,  when  the  tube  is  closed  on  one  side  of 
the  vehicle  and  infinitely  long  in  the  other  direction  there  is  no 
solution  similar  to  those  outlined  above,  if  the  pressure  differ¬ 
ence  given  by  Eq.  (2)  is  considered  to  be  established  when  the 
vehicle  is  still  at  rest.  Such  a  solution  does  exist,  however,  if 
the  tube  is  open  on  one  side  of  the  vehicle  anl  infinitely  long  in 
the  other  direction. 


If  the  region  between  the  open  end  and  the  vehicle  is  considered 
to  remain  at  atmospheric  pressure,  then  for  the  case  where  the  open 
end  is  behind  the  vehicle,  Eq.  (2)  and  the  isentropic  relations  give 

1,[,-  A,f]  (28) 

The  im-ariance  of  the  Riemann  parameter  on  the  characteristics  which 
reach  the  vehicle  from  the  front  gives 


i. 


ef-i  ^ 


tx, 


- 


2 


(29) 


Substituting  these  two  equations  into  the  continuity  equation  across 
the  vehicle  gives 


(30) 


The  corresponding  e  /nations  for  the  case  when  the  open  end  is  ahead 
cf  the  vehicle  and  the  ti?,t.e  is  infinitely  long  behind  it  are; 


A-IA 


/  - 

€ 

0 

(28') 

A  - 
Y'f  “  r-/ 

Azx 

11 

Uv 

(29') 

K'  - 1 } 

■  - 

> 

> 

(30') 

The  speed  of  sound  A  in  the  region  adjacent  to  the  vehicle 
on  the  infinite  tube  side  can  be  obtained  from  either  Eq,  (28)  or 
(28')  for  any  value  of  ^  ,  Once  this  is  determined,  the 

corresponding  value  of  the  flow  velocity  on  both  sides,  'Uix 
Ut^  ,  can  be  obtained  from  either  Eqs.  (29)  and  (30)  or  Eqs*  ( 
and  (30*). 


5)  Both  ends  open 

When  both  ends  of  the  tube  are  open  the  static  pressure  at  th 
outflow  end  and  the  stagnation  pressure  at  the  inflow  end  are  cqua 
to  atmospheric  pressure.  If  the  inlet  and  exit  flows  are  consider 
incompressible  at  constant  temperature,  then 


Equatloitii  and  31  give 


I 


(31) 


(32) 
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and,  from  the  continuity  equation  across  the  vehicle, 


»  ~  ^1 _ 

=  V.  -  (  T,.  \  (33) 

Numerical  Example 

Because  of  the  differences  of  approach  in  the  above  analyses 
of  the  various  configurations  considered,  the  results  could  not  be 
compared  on  the  basis  of  a  common  set  of  equations  or  nondimens ional 
chares.  For  this  reason  a  numerical  exan^le  will  be  given  here  to 
provide  an  indication  of  the  order  of  magnitude  of  the  induced  air 
velocities  produced  xn  the  various  configurations.  Consider  the 
following  vehicle  and  tuber 


^  ) 

L 

cx 

6 


5  100,000  lbs 
3  500  ft /sec  8  340  nph 
3  100  miles 

=  177  ft^ 

3  0.1 


and  the  following  atmospheric  conditions: 


thus 


B  3  2100  psf 

To  3  540®R 

b  =  ^  =  0*0147 

C  =  ZM  =  0.269 


For  these  conditions  the  vehicle  velocity  and  the  time  to 
travel  to  a  given  location  are  plotted  on  Fig.  3  as  functions  of 
the  vehicle  location.  The  corresponding  induced  air -velocities 
the  various  end  conditions  are  presented  on  Figs.  4,  5,  and  6, 
the  maximum  induced  air  velocities*  which  always  occur  immediatf 
in  front  of  and  immediately  at  the  rear  o^  the  vehicle,  are  pre 
seated  on  these  figures. 

The  calculated  air  velocities  for  the  case  of  a  tube  whfch 
open  behind  the  vehicle  and  closed  ahead  of  it  are  not  presente 
on  the  curve  because  the  analytical  model  used  becomes  unrealist 
in  this  case,  for  the  reasons  discussed  previously.  Also,  onl’* 
of  the  half-inf Inlte-tubc  cases  is  shown  on  the  figures;  for  th 
other  case,  only  the  total  pressure  level  throughout  the  tube  wc 
be  changed.  The  pressure  ratio  and  velocity  change  across  the 
vehicle,  and  the  abaolute  value  of  the  air  valoclty  would  remaix 
aane. 


The  maxitaum  induced  air  velocity,  which  occurs  for  a  tube 
at  both  ends,  is  452. b  ft/sec,  which  is  a  speed  at  which  comprei 
blllty  effects  (as  reflected  in  the  ratio  between  total  and  stal 
pressure)  are  still  negligible.  Therefore  the  assunptlon  of  in* 
press ible  flow  is  valid  for  all  caaei  in  which  this  assumption 
isadc. 


A- 17 


It  is  also  noteworthy  that  for  the  esse  of  the  tube  closed  at 
both  ends  changes  of  the  parameter  be  up  to  0.2,  which  would  most 
likely  be  the  maximum  value  desirable  (see  Fig.  2),  the  change  in 
the  induced  air  velocities  is  so  small  that  it  cannot  -be  plotted 
on  the  same  graph.  Also  the  parameter  c  does  not  appear  in  any  of 


the  other  cases 
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DISCUSSION 

The  ressults  of  the  analysis,  as  applied  to  a  specific  numt 
cal  example,  Indicate  that,  as  could  be  expected,  the  lowest  1' 
flow  velocities  are  produced  when  the  tube  is  closed  at  one  or 
ends,  and  that  these  velocities  are  in  either  case  one  order  c 
magnitude  lower  than  those  that  would  be  produced  by  a  vehlcl'^* 
the  same  thrust  schedule  in  a  tube  open  at  both  ends* 

In  a  tube  that  is  infinitely  long  in  both  directions,  t r 
half-infinite  tube  with  an  open  end,  the  same  vehicle  will  ind 
air  velocities  much  more  nearly  equal  to  those  of  the  closed 
case  than  to  those  of  the  open  tube.  This  indicates  that,  alct 
placing  a  gate  behind  the  vehicle  to  produce  a  closed-end  cond. 
would  give  the  best  results,  very  little  would  be  lost  by  not  i 
such  a  gate  if  the  tube  la  ao  long  that  waves  reflected  from  ei 
end  do  not  return  to  the  v^iele  before  the  latter  has  reached 
cruising  speed. 

In  the  case  of  the  short,  open  tube,  the  velocities  induce 
the  vehicle  considered  here  in  the  acceleration  phase  are  of  tl 
same  order  of  magnitude  (but  of  the  opposite  sign)  as  those  ths 
would  be  Induced  by  a  wheel -tract ion  vehicle  -  or  as  those  that 
would  be  required  to  propel  a  pneumatic -dispatch  vehicle  -  at  t 
same  terminal  velocity. 
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Viscous  effects  j,  which  are  neglected  in  this  amilysis,  are 
uniniportant  in  all  cases  except  that  of  the  acceleration  in  proxi¬ 
mity  to  an  open  end,  because  the  maxinium  calculated  flow  velocitiea 
in  all  such  cases  are  actually  very  low,  and  also  because  these 
maximum  velocities  appear  only  as  brief  and  localized  transients* 
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SECTION  B 


THE  EFFECT  OF  SIMULATED 
HIGH  FORWARD  SPEED  ON  A 
TWO-DIMENSIONAL 
GROUND  EFFECT  SUPPORT 


by 

Anthony  A.  West 


SUMMARY 


Preliminary  tests  have  been  conducted  to  investigate 
the  feasibility  of  using  ground  effect  supports  to  provide 
suspension  and  propulsion  of  a  ”ehicle  at  high  forward 
speeds.  The  effects  of  simul^'^d  high  forward  speeds  on 
the  Ifft  and  drag  of  a  two-dimensional  ground  effect  support 
are  correlated  through  a  forward-speed  parameter  which  is 
derived  in  this  report. 
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SYMBOLS 


A  augmentation  ratio  for  GES  configuration 

c  chord  length  of  model  7.5  in. 

D 


i-H 

Mr/ 


Cj)  drag  coefficient  = 
lift  coefficient  = 


^  Q. 
L 


Cn  =  0.07 


Cl'  aerodynamic  lift  coefficient 
Cd  recovery  factor  - 

Mr/ 

jet  "omentum  coefficient  [ GESjj^ 

D  drag 

h  height  of  model  from  groundboard 

J  nozzle  length  -  length  of  exposed  model  span  ■  9.75  in. 


L  total  lift 

Ljj  lift  at  zero  forward  speed  for  GES 
L'  aerodynamic  lift  .  L  ~  for  GES 
M  momentum  flux/unit  length  of  nozzle  ■ 

Mx  total  momentum  flux/unit  length  of  nozzle  ■  Mp-f  for  GES 
ih  total  nozzle  ness  flow  rate  meeaured  from  venturi  meter 

I 

p  pressure 

q  dynamic  head  ■  (3  ^  ’ 

S  exposed  model  plan  area  »  Jc  •  0.508  £t^ 

V  velocity 

Vi  Jet  velocity  calculated  assuming  the  air  in  the  chambers 
^  is  expanded  Isoenergetically  and  isentroplcally  down  to 
the  free-atream  static  pressure 


ill 


V. 


wind  tunnel  free  streein  velocity  ■  simulated  forward  speed 
of  model 


0  jet  angle  (see  Fig.  4) 

p  air  mass  density 

^  forward  speed  parameter  -  defined  in  Appendix 
Subscripts 

c  critical  condition  at  which  Jets  are  swept  under  model 
(se?  Figs.  1(c)  and  !(£)•) 

F  front  j  et 

0  no  air  supplied  to  nozzles 

R  rear  j  et 

wind  tunnel  free-stream  condition 


Superscript 

o  stagnation  condition 
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I.  INTRODUCTION 

In  a  report  on  a  proposed  high-speed  self-propelled  air- 
breathing  vehicle  in  a  tube,  Foa  (Ref.  3)  has  suggested  that 
ground  effect  supports  could  be  used  for  suspension  of  the 
vehicle.  The  tube  wall  would  provide  the  kind  of  transverse 
guidance  and  control  that  is  absent  in  ground-effect  supports 
operating  over  open  terrain  or  water  surfaces.  Furthermore, 
the  rigidity  of  the  tube  wall  should  permit  much  higher  velo¬ 
cities  in  the  air  curtains  than  are  permissible  over  soft 
ground  or  water,  and  should  therefore  make  possible  the  utili¬ 
zation  of  ground-effect  support  at  very  high  forward  speeds. 

At  the  present  time,  however,  there  is  little  information  avail¬ 
able  on  the  effects  of  high  forward  speeds  on  the  air  cushion  of 
a  ground  effect  support  and  the  present  work  was  undertaken  to 
provide  a  fundamental  understanding  of  the  problem.  For  this 
purpose  it  was  believed  that  a  two-dimensional  rather  than  a 
thrae-dimensionel  study  would  be  more  fruitful.  The  specific 
aiins  of  this  preliminary  investigation  were  to  determine; 

(a)  the  feasibility  of  ground  effect  support 
at  high  subsonic  forward  speeds. 

(b)  the  possibility  of  utilizing  the  momentum 
of  the  ground  effect  curtain  in  the  genera¬ 
tion  of  thrust. 

The  work  on  this  project  has  been  carried  out  under  the 
auspices  of  the  Army  Research  Office  (Durham)  Grant  Mb.  DA-ORO- 
31-124-61-688. 


II.  DESCRIPTION  OF  APPARATUS 


Wind  Tunnel  and  Balance 

The  dynamic  expar-fments  on  the  model  were  carried  out  in 
the  RPI  4x6  ft.  subsonic  wind  tunnel,  whidh  is  of  the  closed- 
throat  return  type.  The  wind  tunnel  control  panel  and  test 
section  are  shown  in  Figs,  B-2  and  B-3  respectively.  The  maxt- 
mum  airspeed  in  the  test  section  was  200  fps.  The  aerodynamic 
forces  on  the  model  were  measured  on  a  six-coo^onent  mechanical 
balance. 

Model 

The  model  is  two-dimensional,  and  its  section  consists  of  a 
flat  undersurface  and  a  circular  arc  top  surface  of  11.5  in, 
radius,  as  shown  in  Fig,  B-4,  The  air  is  supplied  to -^rhe  front 
and  rear  chambers  independently,  and  is  distributed  within  each 
chamber  by  means  of  the  distribution  tubes  which  contain  a  0.03 
in.  longitudinal  slit.  Two  pressure  taps  in  each  chamber  are 
located  in  the  base  plate  to  measure  the  stagnation  pressure  in 
the  chamber.  A  small  block  is  welded  to  the  base  plate  just 
'^downstream"  of  these  taps  to  ensure  that  the  flow  is  brought  to 
rest  at  the  tap.  The  leading-edge  and  trailing»edge  nozzle  gaps 
are  adjustable  but  were  maintained  at  0.05  ins.  for  the  present 
series  of  tests.  The  dimensions  of  the  model  are  given  in  Figs. 


B-4  and  B-5. 
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The  model  was  mounted  in  the  test  section  on  two  main 
struts  as  shown  in  Fig.  B-6c  The  angle  of  attack  was  0°  for 
all  tests.  The  model  tips  were  recessed  "into  the  sideboards 
which  were  carefully  cut  to  prevent  the  model  fouling  the 
boards.  The  sideboards  were  attached  to  the  tunnel  floor  by 
slotted  angle-iron  and  the  groundboard  was  mounted  between  the 
sideboards  by  four  brackets  as  shown  in  Fig.  B-7.  The  clearance 
between  the  model  and  the  groundboard  could  be  adjusted  from 
0-18  in.  The  dimensions  of  the  sideboards  and  groundboards  used 
during  the  tests  are  shown  in  Fig.B-8. 

Jet  Air  Supply 

High  pressure  air  was  supplied  from  a  reservoir  to  a  manual 
regulator  which  in  turn  supplied  the  air  to  the  front  and  rear 
distribution  tubes  by  means  of  the  air  supply  pipes  as  shown  in 
Fig.  B-9.  A  gate-valve  was  fitted  in  each  of  these  air  supply 
pipes  so  that  by  suitable  adjustment  of  these  gate  valves  and  the 
manual  regulator  any  desired  mass  flow  to  the  front  and  rear 
chambers  could  be  obtained  Independently.  A  venturi  meter  was 
fitted  to  each  supply  pipe  to  measure  mass  flow  supplied  to  the 
chambers . 


III.  EXPERIMENTAL  PROCEDURE 

Measurement  of  Momentum  Flux  and  Jet  Angle  Along  the  Nozzles 

In  consideration  of  the  possibility  of  flow  changes  due  to 
structural  deformation  under  pressure,  the  jet  angle  and  the 


distribution  of  the  momentum  flux  along  the  leading-edge  and 
tralllng-edge  nozzles  were  measured  for  various  chamber  pres¬ 
sures  using  the  apparatus  shown  in  Fig.  B-10.  The  jet  angle 
was  82°  for  both  nozzles  and  was  found  to  be  invariant  with 
chamber  pressure.  The  maximum  variation  of  jet  angle  along 
the  nozzle  was  i:  2°.  It  was  found,  for  both  nozzles,  that  the 
momentum  flux  at  the  midspan  position  was  a  maximum  of  5% 
greater  than  the  momentum  flux  at  the  outer  ends.  It  was  assur.i: 
that  this  variation  was  sufficiently  small  to  be  neglected  for 
the  purposes  of  this  series  of  tests. 

Determination  of  Chamber  Pressure 

The  rig  shown  in  Fig.  B-10  was  used  with  the  pitot  tube  at 
the  jet  exits  to  determine  the  variation  of  the  stagnation  prej 
sure  along  the  span  for  a  series  of  different  chamber  pressures. 
These  readings  were  compared  with  those  obtained  from  the  Interr 
pressure  taps  in  the  baseplate.  The  readings  agreed  within  4^ 
and  a  callbratlcn  curve  was  constructed  to  convert  the  chamber 
pressure  readings  into  the  everage  value  of  the  stagnation  pres¬ 
sure  along  the  noxzle. 

Wind  Tunnel  Tests 

After  the  preliminary  experiments  described  above,  the  mode 
was  mounted  in  the  wind  tunnel  as  shown  in  Fig.  B-7.  The  lift  a 

drag  on  the  model  were  measured  at  0°  incidence  for  different  al 
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speeds  in  the  test  section,  chamber  pressures,  heights  above 
ground,  and  groundboards .  During  each  run  the  air  speed  in  the 
test  section  was  changed  from  0  to  200  fps  and  the  mass  flow  to 
the  chambers  was  kept  constant  by  adjustment  of  the  regulator. 

An  equal  mass  flow  of  air  was  supplied  to  the  front  and  rear 
nozzles  In  those  tests  in  which  the  model  was  tested  as  a  ground 
effect  support.  The  test  conditions  are  summarized  in  Table  B-I. 

IV.  RESULTS  AND  DISCUSSION 

Fig.  B-11  shows  that  there  is  close  agreement,  both  in  trend 
and  magnitude,  between  the  experimental  results  and  those  pre¬ 
dicted  by  the  simple  momentum  theory  of  Reference  1.  The  graph 
also  shows  that  an  error  of  a  few  degrees  in  the  measurement  of 
the  jet  angle  can  result  in  appreciable  changes  in  the  value  of 
the  augmentation  ratio  and  it  is  believed  that  this  is  the  major 
cause  of  the  discrepancy  between  the  theoretical  and  experimental 
results. 

Fig.  B-12  shows  how  the  total  lift  acting  on  the  model  in¬ 
creases  with  forward  speed  in  the  ground  effect  support  and  jet 
flap  wing  configurations. 

Fig.  B-13  reveals  that  the  aerodynamic  lift  coefficient  in¬ 
creases  rapidly  with  increasing  forward  speed  in  the  subcritical 
regime,  reaches  a  maximum  around  the  critical  forward  speed,  and 
thereafter  decreases  slowly  in  the  supercritical  regime.  This 


graph  confirms  the  fundamental  correctness  of  the  chosen  forward 
speed  parameter.  It  will  be  noticed  that  for  Run  1  the  maximum 
value  of  Ci  Is  attained  at  a  value  of  the  forward  speed  para¬ 
meter  greater  than  1.0.  The  reason  for  this  Is  that  the  forward 
speed  parameter  Is  based  on  the  free-stream  velocity  but  the  model 
Is  In  the  boundary  layer  of  the  groundboard.  This  Is  borne  out  by 
a  comparison  of  Runs  1  and  2.  The  jet  velocity  Is  the  same  for 
both  runs  but  In  Run  2  the  model  Is  further  from  the  groundboard. 

From  Fig.  B-14  It  can  again  be  seen  that,  with  the  exception 
of  Run  6,  the  maximum  value  of  the  aerodynamic  lift  coefficient  xs 
attained  close  to  a  value  of  1,0  of  the  forward  speed  parameter. 

The  effect  of  the  broundboard  boundary  layer  is  very  pronouced  in 
Run  6,  which  was  performed  using  a  very  low  jet  velocity.  With 
this  low  jet  velocity,  the  wind  tunnel  airspeed  at  the  critical 
condition  is  relatively  small  and  thus  the  boundary  layer  on  the 
groundboard  is  fairly  thick.  This  can  be  seen  by  a  con^arlson  of 
Runs  2  and  8  In  which  a  higher  jet  velocity  was  used;  both  runs 
were  perfomad  with  the  model  at  the  same  distance  above  the  ground 
board. 

Fig.  B-15  shows  the  effect  on  the  aerodynamic  lift  coefficient 
of  changing  the  groundboara  size.  The  jet  velocity  and  the  dis¬ 
tance  of  the  model  from  the  groundboard  was  the  same  for  both  runs. 
The  use  of  the  larger  groundboard  with  Its  attendant  thicker 
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boundary  layer,  has  two  consequences  —  firstly,  the  maximum 
value  of  cl  Is  attained  at  a  higher  value  of  the  forward  speed 
parameter  for  the  reasons  given  above.  Secondly,  the  magnitude 
of  the  maxJjnum  value  of  Cl  Is  less  because  tha  total  lift  Is 
less  due  to  the  model  being  In  a  lower  velocity  boundary  layer. 

For  the  above  reasons  the  result  ’j,  obtained  using  the  large  and 
small  groundboards  are  presented  8eparatel> . 

Fig.  B-16  confirms  that  for  the  jet  flap  wing  configuration 
out  of  ground  effect  there  is  a  linear  relationship  between  the 
lift  coefficient  and  the  radical  of  the  jet  momentum  coefficient  — 
a  result  already  reported  in  Reference  2. 

» 

i 

For  the  ground  effect  support  configuration  in  the  super¬ 
critical  regime,  Fig.  B-17  reveals  a  similarly  unique  relationship 
between  the  aerodynamic  lift  coefficient  and  the  radical  of  the 
jet  momentum  coefficient,  Irrespective  of  the  actual  magnitude  of 
the  jet  velocity.  This  would  suggeat  that  the  aerodynamic  lift 
experienced  by  the  model  In  the  supercritical  regime  Is  derived 
in  the  same  way  as  the  lift  on  a  jet  flap  wing  out  of  ground 
effect  l.is.  mainly  from  the  low  pressure  region  on  the  top  surface 
of  the  model.  The  airflow  in  these  two  cases  Is  shown  in  Figs. 
B-l(d)  and  B-l(g)  for  the  ground  effect  support  and  for  the  jet 
flap  wing,  respectively.  The  above  is  only  a  tentative  suggestion 
and  would  need  to  be  checked  by  a  survey  of  the  pressure  distribu¬ 
tion  around  the  model  at  various  forward  speeds. 


Fig.  B-18  presents  typical  results  for  the  model  In  the 
ground  effect  support  configuration.  It  will  be  noticed  that, 
although  the  total  drag  Is  less  than  the  profile  drag,  the  net 
force  Is  still  a  drag.  In  none  of  the  tests  carried  out  was  c 
net  thrust  experienced  by  the  model.  The  reason  for  this  is  t 
high  value  of  the  profi'^.e  drag  resulting  from  the  relatively 
thick  section  and  blunt  trailing  edge  of  the  model. 

Fig.  B-19  reveals  a  rapid  increase  of  the  thrust  recovery 
coefficient  with  forward  speed  in  the  subcrltical  regime.  Thi 
is  probably  mainly  due  to  the  fact  that  as  the  forward  speed 
increased  more  of  the  jet  curtain  is  swept  under  the  model  ra'; 
than  escaping  over  the  top  surface  of  the  model.  As  can  be 
from  a  conparison  of  Figs.  B-l(b)  and  B-l(d)  the  possibility  c 
thrust  recovery  is  inherently  greater  for  the  smoother,  leas'd 
slpative  airflow.  It  can  be  seen  from  this  graph  that  is 
also  a  function  of  the  ratio  of  forward  speed  to  jet  velocity. 
It  will  be  noted  that  valjes  of  C/f  greater  than  1.0  are  obtai 
in  the  supercritical  regime.  It  is  suggested  that  this  is  duo 
the  airflow  over  the  top  of  the  model  being  entrained  by  the  r 
jet  as  shown  In  Fig*  B->l(d).  The  Increase  In  with  forward 
speed  In  the  supercritical  regime  Is  due  to  the  fact  that  the 
losses  associated  with  the  mixing  of  two  alratreams  decrease  ai 
the  difference  between  their  velocities  Is  decreased. 
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Fig,  B“20  exhibits  more  clearly  the  strong  dependence  of 
the  thrust  recovery  factor  on  the  ratio  of  the  forward  speed 
to  the  jet  velocity.  The  scatter  in  the  experimental  data  is 
due  to  the  smallness  of  the  forces  being  measured. 

V.  CONCLUDING  COMMENTS 

Depsite  the  complications  arising  from  the  interaction 
between  the  jet  curtain  and  the  groundboard  boundary  layer,  it 
is  believed  that  the  following  valid  conclusions  for  a  two- 
dimensional  model  may  be  drawn  from  the  experimental  data: 


1.  For  a  ground  effect  support  employing  a  thin 
jet  curtain,  the  lift  at  zero  forward  speed  is  satis¬ 
factorily  predicted  by  the  theory  of  Ref.  1. 

2.  For  constant  height  and  constant  jet  momen¬ 
tum  flux,  the  lift  increases  with  forward  speed  both 
in  the  ground  effect  support  and  jet  flap  wing  con¬ 
figurations. 

3.  The  forward  speed  parameter  which  is  developed 
in  the  Appendix  is  shown  to  be  significant. 

4.  In  none  of  the  test  runs  was  there  a  net 
thrust  experienced  by  the  model.  However,  the  feasi¬ 
bility  of  utilizing  the  momentum  flux  from  an  "aero- 
dynamlcally  clean”  ground  effect  support  as  a  signi¬ 
ficant  contribut  on  to  thrust  has  been  demonstrated. 

5.  For  max,. mum  thrust  recovery  the  ideal  operating 

condition  is  apparently  obtained  by  operating  near  the 
critical  point  and  arranging  for  the  jet  velocity  to  be 
equal  to  the  forward  speed  of  the  support.  This  conclu¬ 
sion  requires  verification  from  data  obtained  for  values 
of  Vo./yi  >  1.0  (see  Fig.  B-20). 
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6.  In  the  ground  effect  support  configuration 
tests,  the  flows  were  essentially  incompressible  as 
the  maximum  value  of  the  free~stream  and  jet  flow 
Mach  numbers  was  0.25.  However  it  should  be  realized 
that  high  forward  speeds  were  simulated  by  values  of 
\^o/Vj‘  up  to  1,5.  It  is  believad  that  the  derived 
forward  speed  parameter  would  be  significant  for  high 
subsonic  speeds  as  compressibilit r  effects  would  not 
change  the  fundamental  nature  of  :he  flow. 


’Hie  above  remarks  are  only  strictly  applicable  to  two- 
dimensional  models  and  to  the  plane  of  syrrinetr'y  of  three- 
dimensional  models.  However,  there  is  nc  reason  to  believe  that 
a  similar  approach  to  the  problem  of  a  three-dimensional  model 
at  high  subsonic  speeds  would  not  yield  qualitatively  similar 
results. 
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SUMMARY  OF  TEST  CONDITIONS 

Run  Ho. 

.Groundboard 

i  h  (ina) 

Chaniter 

(lb/ft2 

Pre&Bure 

Gaage) 

Rfcer 

Remarks 

1 

Small 

0.38  * 

56 

56 

'ground  If  fee  C  : 

2 

tt 

0.69 

56 

56 

ff 

3 

ff 

1.00 

47 

47 

ff 

4 

r» 

2.00 

120 

120 

fw 

5 

Large 

0.69 

56 

56 

rr 

6 

8t 

1.63 

17 

17 

tr 

7 

rt 

1.63 

52 

52 

ft 

8 

ff 

1,63 

78 

78 

ft 

9 

ft 

2.50 

120 

120 

ti 

10  ' 

■n 

4.60 

0 

197 

Jet  flap  Wing 

11 

tt 

- 

8.50 

0 

78 

ft 

12 

ft 

12.50 

0 

310 

ff 

13 

ft 

12.50 

0 

78 

tr 

14 

tt 

12.50 

0 

318 

rr 

15 

tf 

16.60 

0 

78 

ff 

16 

ft 

16.60 

0 

470 

ff 

17 

ff 

4.60 

0 

0 

* 

determination  of 

18 

tf 

16.60 

0 

0 

ft 
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APPENDIX 

The  nature  of  the  airflow  around  the  model  is  character¬ 
ized  by  a  ''forward-speed  parameter'*  ^  defined  as  the  ratio 

^  _  height  of  model  above  ground _  „ 

^  critical  height  oi  model  above  ground  “  ^ 

It  is  assumed  that  at  the  critical  forward  speed  the  pressure 
acting  on  the  upstream  side  of  the  front  jet  curtain  is  the 
freest  ream  stagnation  pressure  .  It  is  further  assumed 

that  the  free  stream  static  pressure  He  acts  on  the  down¬ 
stream  side  of  the  rear  jet  curtain. 

FREE  Qf 

STREW*"  -  T-t  ■ 

♦  I 

Applying  the  momentum  theorem  to  the  control  volume  shown  and 

assuming  Invlacld  flow; 

*  tif!  ( I  -  cos  9ff)  +  f^R  ( I  *■  cos  &„)  ^  fLhc 


1h«  forward  ivpoad  paraoMitor  is  then 


\^(  pu  - _ 

cat  £(,)  +  HrO  ♦  cos  Bij 

•  I 

There  are  two  special  cases  to  be  considered: 

(1)  The  Ground  Effect  Support  Configuration 
In  this  case  the  momentum  flux  issuing  from  the  front  and 
rear  jets  is  the  same  l.e.  =*  Hr  •  The  front  and  rear 

jet  angles  are  also  the  sane  0^  ^  Bn 


$  * 


h(Pi  -  pj 
tif  * 


(2)  The  Jet  Flap  Wing  Configuration 


In  this  case  there  la  no  front  Jet  l.e.  Mp  *  0, 
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SUMMARY 


Preliminary  tests  have  been  carried  out  to  investigate 
the  pressure  distribution  on  the  outside  surface  of  a  vehicle 
moving  at  high  speed  within  a  tube.  The  vehicle  was  simulated 
by  a  hollow  cylinder  mounted  in  a  tube  of  circular  cross 
section.  The  contribution  of  these  surfaces  to  the  lift  and 
pitching  moment  on  the  vehicle  are  determined.  The  contribution 
to  pitching  moment  is  found  to  be  destabilizing. 

These  tests  are  supplemented  by  flow  visualization  tests 
at  low  flow  velocities. 
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SYMBOLS 


c  length  of  models .  For  all  models  =  7” 

d  outside  diameter  of  models.  For  all  models  ■  3" 

S  model  planform  area  cd,  (see  Appendix) 

local  lift  coefficient 

Ci  =  ^Cfy(x/c)cc<<;/> 


overall  lift  coefficient;  average  found  from 
vs  ^  curves 


^\.E. 


Iccal  pitching  rsionient  coefficient 


\.z. 


overall  pitching  moment  coefficient  about  leading 

edge.  Average  found  from  Cm  vs  (p  curves 

ii.  E. 

pressure^  coefficient  at  position  x 

c.  =  /  -  ^ 

total  lift  due  to  external  pressure  distribution. 

Here  lift  is  defined  to  act  in  that  plane  which  contains 
the  axis  of  the  cylinder  both  before  and  after  any 
change  of  disturbance  in  angle  of  attack. 

total  moment  about  the  leading  edge  due  to  external 
pressure  distribution. 

flow  Mfffh  number 


Ptank  manometer  reading  when  connected  to  atmosphere  at  one 
end  and  to  settling  tank  at  the  other. 
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manometer  reading  when  connected  to  atmosphere  at 
one  end  and  a  model  pressure  tap  at  the  other 

flow  velocity 

1  2 

dynamic  pressure,  ^  C  * 

See  Figure  C-7 

angle  of  attack  in  pitch,  see  Figure  C-7 

compound  angle  of  atta  ,  i.e.,  effective  angle 
of  attack  involving  be-  ^  pitch  and  yaw,  see  Appendix 

angle  of  attack  in  yaw,  see  Figure  C-7 

circumferential  position  being  simulated  by  placing 
model  at  a  compound  angle,  see  Figure  C-7 
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I.  INTRODUCTION 

In  a  report  on  a  proposed  high  speed,  self-propelled,  air- 
breathing  vehicle  in  a  tube,  Foa  (Ref.  1)  has  suggested  that 
ground  effect  supports  could  be  used  for  suspension  of  the  vehicle. 
West  (Ref.  2)  has  investigated  further  the  feasibility  of  this 
system  to  provide  both  suspension  and  propulsion  of  the  vehicle  at 
high  forward  speeds.  One  criterion  in  determining  the  required 
strength  of  the  ground  effect  support  is  the  stability  of  the 
vehicle  when  disturbed  from  its  axi-symmetric  position  within  the 
containing  tube. 

This  report  gives  the  results  of  a  preliminary  program  to 
investigate  the  stability  of  the  vehicle.  At  the  time  of  writing, 
the  shape  of  the  vehicle  was  not  clearly  defined.  Hence  a  very 
simple  shape,  i.e.,  a  hollow  cylinder,  was  used.  It  is  with  the 
pressure  distribution  on  the  outside  surface  of  the  vehicle  that 
this  report  is  concerned. 

II-  DESCRIPTION  OF  APPARATUS 

Two  separa  a  series  of  experiments  were  made. 

1.  Pressure  tap  models. 

Two  models  were  tested;  one  with  a  blunt  leading  edge  (Model 

I)  and  the  other  with  a  round  leading  edge  (Model  II).  The  models 

* 

and  apparatus  are  shown  in  Figs.  C-1  thru  C-7.  Fig.  C-1  shows 
Model  I  and  the  location  of  its  six  pressure  taps.  Fig.  C-2  shows 


Model  II  and  the  accompanying  mechanism  for  incidence  adjustmen 
and  measurement  in  pitch  and  yaw.  Figs.  C-3  thru  C-5  show  the 
detailed  configuration  and  dimensions  of  the  two  models.  Fig.C- 
shows  the  model  mounted  centrally  within  its  containing  tube. 
This  tube  had  a  bellmouth  inlet  open  to  the  atmosphere  and  an 
airflow  induced  by  operating  a  compressor  downstream  of  a  settl 
tank.  The  manometers  were  external  to  the  rest  of  the  test  c:.-^ 
Since  there  must  necessarily  be  some  clearance  between  the  stin 
and  the  tube,  an  inflow  of  air  and  resultant  distortion  of  t'  e 
main  airflow  over  the  model  was  of  some  concern.  Originally, 
tests  were  run  in  a  cell  in  which  a  compressor  was  mounted  up¬ 
stream  of  the  model.  As  a  result  the  static  pressure  at  the  fi'c 
(which  was  fairly  close  to  the  duct  exit)  was  equal  to  the  arobi 
static  pressure  and  hence  leakage  losses  about  the  model  sting 
were  negligible.  Mechanical  failure  of  this  facility  necessita* 
ted  Installation  of  the  model  in  the  present  cell.  Consequentlj 
the  static  pressure  in  the  region  near  the  model  was  no  longer 
at  the  ambient  static  value  and  leakage  became  appreciable,  glv: 
rise  to  the  need  for  a  sealing  shield,  as  indicated  in  Fig.  C-6. 

2.  Flow  visualization  model. 

This  model  was  made  from  a  transparent  plastic  tube  as  shov 
in  Fig.  C-8;  it  was  the  sau.u  size  as  the  pressure  tap  model  and 
also  mounted  in  a  6”  diameter  containing  tube.  Provision  was  ma 
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to  alter  the  angle  of  attack  in  pitch  only.  The  whole  apparatus 
was  mounted  in  the  Rensselaer  Polytechnic  Institute  low  speed, 
low  turbulence  smoke  tunnel.  A  flow  velocity  of  20  ft/sec  was 
used  in  all  flow  visualization  tests. 

III.  EXPERIMENTAL  PROCEDURE 
1.  Pressure  tap  model. 

A  traverse  was  made,  across  the  working  section  of  the  tube 
where  the  model  would  be,  to  measure  static  and  stagnation  pres¬ 
sure.  Once  out  of  a  reasonably  thin  boundary  layer  these  para¬ 
meters  were  constant  thus  ensuring  a  uniform  Mach  number. 

We  wish  to  measure  the  pressure  on  the  outside  surface  of 
the  model  at  as  many  places  as  possible.  On  the  other  hand  model 
complication  and  strut  interference  increase  as  the  number  of  taps 
is  increased.  It  was  possible,  however,  to  obtain  a  complete  sur¬ 
vey  of  pressures  around  the  model  at  any  given  effective  angle  of 
attack  by  adopting  the  following  procedure; 

a)  the  model  was  .set  at  ckg  in  pure  pitch;  then  6  more 

readings  were  obtained  by  using  —  pure  pitch. 

b)  the  model  was  set  at  in  pure  yaw.  Due  to 

symmetry  thivS  gives  the  same  result  as  —  in  yaw. 

Thus  the  pressure  at  six  lengthwise  stations  can  be  determined  at 

4  circumferential  positions. 

Now,  due  to  the  fact  that  a  plane  of  symmetry  for  the  model 
and  the  containing  tuhiJ  can  always  be  made  to  exist  through  the 


axis  of  the  model,  it  is  possible  by  the  use  of  compound  rota¬ 
tions  to  measure  the  pressure  ac  any  circumferential  position. 

All  rotations  should  be  about  the  geometric  center  of  the  model. 
For  pitch  rotations  this  was  the  case  but  lor  yaw  rotations  a 
small  error  was  incurred  due  to  the  fact  that  the  axis  of  rota¬ 
tion  could  not  be  made  to  pass  through  the  geometric  center. 
Simple  geometric  formulae  were  derived  which  gave  the  required 
pitch  and  yaw  angles  for  any  given  effective  angle  of  attack 
and  circumferential  position,  (j>  .  (See  Appendix). 

Tests  were  made  at  two  different  Mach  numbers,  and  for  two 
different  effective  angles  of  attack,  for  both  the  blunt  and  the 
rounded  leading  edged  models. 

Some  distortion  of  the  flow  field  results  from  the  presence 
of  the  sting  in  the  airflow.  This  distortion  was  kept  to  a 
mlntmum  by  suitable  sheplng,  see  Figs.  C-1  thru  C-5. 

2.  Flow  vlsuallaatlon  model. 

The  crux  of  the  problem  being  investigated  Is  that  the  flow 
field  around  and  through  the  model  is  influenced  by  the  contain¬ 
ing  tube  in  a  manner  which  is  not  qualitatively  predictable  on 
the  ba^is  of  existing  theories  or  previous  observations. 

To  get  some  idea  of  the  disturbed  airflow  pattern,  the 
transparent  flow  visualization  model  was  constructed  to  enable 
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photographs  to  be  obtvained.  Liquid  titanium  tetrachloride  was 
introduced  in  small  quantities  into  the  recess  of  a  small,  smooth 
■ake  some  distance  upstream  of  the  model.  The  air  passing  over 
the  rake  caused  a  steady,  white  stream  of  gaseous  titanium  tetra¬ 
chloride  to  be  drawn  downstream.  The  vertical  position  of  the 
rake  relative  to  the  model  was  adjustable  and  the  airflow  over 
several  parts  of  the  model  could  be  visualized. 

IV.  RESULTS  AND  DISCUSSION 
L.  Pressure  survey  tests. 

In  order  to  calculate  the  lift  and  pitching  moment  on  the 
aodel  due  to  the  pressure  distribution  on  its  outside  surface  a 
)lot  of  Cp  vs  X/C  is  made.  As  we  are  dealing  with  an  essentially 
;hree  dimensional  flow  a  scheme  must  be  adopted  to  produce  over- 
ill  values  of  and  from  the  many  Cp  values. 

There  must  be  a  stagnation  point  at  or  extremely  near  to 
he  leading  edge  of  the  model  around  tha  whole  circumference, 
bus  assuming  an  isentropic  flow  for  the  stagnation  streamtube 
c  can  fix  Cp  ■  1  at  the  leading  edge  throughout.  Secondly, 
ssume  that  the  pressure  inside  the  model  is  uniform  at  any  length^* 
ise  station  over  the  whole  cross  section;  'no  restriction  need  be 
ut  on  the  lengthwise  variation  of  pressure.  This  is  believed  to 
e  a  valid  assumption  over  most  of  the  length  of  the  tube,  and  is 
urther  supported  by  Fig.  C-9a,  where  an  internal  streamline  is 
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seen  to  become  parallel  the  walls  almost  instantaneously  on 
entering  the  model  and  to  keep  in  this  position  throughout.  In 
particular,  the  pressure  at  the  trailing  edge  of  the  model  is 
Uniform.  Now,  with  the  abc-ve  conditions,  pairs  of  Cp  vs  X/C 
curves  can  be  plotted  together  to  form  closed  curves.  This  has 
been  done  in  Figs.  C-11  and  C-12  using  Cp  values  on  those  longi¬ 
tudinal  strips  having  a  common  projection  on  a  plane  bisecting 
the  cylinder  into  its  upper  and  lower  surfaces;  i.e.,  results 
for  ^  j  and  (j)  =■  (or  for  one  value  of  Z/C)  are  plotted 

together  (see  Appendix).  The  area  enclosed  by  these  curves  gives 
the  sum  of  the  forces  normal  to  the  surface  at  the  two  points 
chosen.  The  resultant  local  lift  coefficient  is  easily  found  fr':'' 
this  ^urve. 

Figs.  C-13  and  C-14  show  these  values  versus  radial  dis¬ 
tance  measured  from  the  center  of  the  tube  along  the  bisecting 
plane  (sec  Pig.  C-7).  The  overall  value  is  found  as  the  aver¬ 
age  The  Included  table  suanarltes  the  values  so  obtained. 

Fig.  C-15  shows  vs  for  both  models.  As  would  be 

expected  the  rounded  leading  edge  model  generated  more  lift  than 
the  blunt  leading  edge  model  throughout.  Both  curves  appear  tc 
be  slightly  nonlinear,  which  is  not  surprising  in  view  of  the 
essentially  three  dimensional  nature  of  the  flow.  This  is  ftirthet 
emphasized  by  Figs.  C-13  and  C-14  that  reveal  that  a  fundamental 


change  in  the  external  flow  pattern  must  take  place  around  the 
models  between  >  5°  and  o<^=  10*^  in  order  to  produce  such 

a  change  in  the  shape  of  the  loading  or  values. 

The  value  of  Cjjj  is  determined  as  the  moment  of  the  area 
under  the  Cp  vs  X/C  curve  with  respect  to  a  suitable  position 
along  the  length  of  the  model;  the  result  being  factored  by 
cos  (j).  It  is  convenient  to  use  the  leading  edge  as  a  datum  for  the 
measurement  of  moment  arms. 

Figs,  G-16  and  C-17  demonstrate  the  variation  of  C  with 

“l.e. 

radial  position  from  the  center  of  the  tube.  Again  the  fundamental 
changes  between  =  5°  and  =  10°  are  seen. 

Fig.  C-18  shows  vs  for  both  models.  In  accord- 

Ti .  F . 

ance  with  expectation.  Model  II  everywhere  develops  larger  moments 
than  Model  I. 

Finally,  Figs.  C-19  and  C-20  show  that  the  pressure  distri¬ 
bution  on  the  external  surface  of  the  model  has  a  destabilizing 
effect.  The  internal  flow  may  modify  this  situation.  It  should 
be  noted  that  the  model  with  the  rounded  leading  edge,  although 
more  acceptable  aerodynamics lly,  could  also  be  more  unstable,  at 
least  as  far  as  the  external  pressure  field  is  concerned,  for 
certain  locations  of  the  center  of  gravity. 

Peripheral  pressure  distributions  represented  in  the  form 
of  polar  plots  of  Cp  vs  ^  are  shown  in  Figs.  C-22  thru  C-29. 
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Results  are  Included  for  tests  run  at  M  =  0.17  and  M  •  0.26  with 
oCg  .  5°  and  10°,  for  both  Models  I  and  II. 

2.  Flow  Visualization  Tests 

The  results  are  shown  as  Figs.  C-9  and  C-10.  Fig.  C-9  shows 
three  smokelines  with  the  model  set  at  o<.  =  7°.  The  auxiliary 
sketches  show  the  position  of  the  smokeline  before  it  is  influ¬ 
enced  by  the  model,  as  seen  from  upstream. 

In  Fig.  C-9a  the  smokeline  was  divided  by  allowing  it  to 
impinge  on  the  leading  edge.  The  early  and  continuing  alignment 
of  the  Internal  streamline  parallel  to  the  model  is  noteworthy; 
separation  occurs  at  the  blunt  leading  edge  as  would  be  expected. 

Fig.  C-10  shows  extreme  angles  of  oc  s  17°  and  18°  res¬ 
pectively,  and  consequently  large  displacements  of  the  smokeline. 
Note,  however,  that  the  sookellne  reverts  to  being  parallel  with 
the  containing  tube  very  such  sooner  at  this  large  vc  than  for 
»c  .  7°. 

As  noted  in  the  introduction,  these  experiments  were  related 
to  a  proposed  vehicle  moving  along  a  tube.  Our  simulation  here 
is  not  quite  complete  as  we  have  the  whole  airflow  having  a  high 
velocity  relative  to  the  tube  whereas  in  the  proposed  system  only 
the  air  in  the  region  near  the  vehicle  has  a  significant  velocity 


relative  to  the  tube. 
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V.  CONCLUSION 

The  results  of  these  tests  clearly  show  the  vehicle  to  be 
unstable  when  displaced  from  its  equilibrium  position.  This 
conclusion  is  a  result  of  tests  concerning  the  external  flow  of 
the  vehicle  only.  It  is  not  obvious  at  the  time  of  writing 
whether  the  resultant  of  all  moments  acting  on  the  vehicle,  con¬ 
sidering  both  its  internal  and  external  flows,  is  stabilizing  or 
not. 


The  quantitative  evaluation  of  the  instability  of  the 
/ehicle  will  provide  useful  criteria  for  the  design  of  the 
support  system  (proposed  to  be  of  the  ground  effect  type),  which 
Ls  itself  stabilizing. 

Nevertheless,  in  view  of  the  inherent  assumptions  and  ap- 
>roximations  used  in  obtaining  the  results  reported  here,  these 
■esults  should  be  accepted  as  quantitative  information  only  in 
■elatJ.on  to  the  idealized  model  by  which  the  vehicle  was  simulated 
nd  to  the  specific  conditions  under  which  the  tests  were  perforro- 
id. 


At  the  time  of  writing,  tests  are  being  conducted  on  a 
ehicle  of  the  type  discussed  in  this  report,  in  which  forces  and 
oments  are  being  measured  directly.  Early  results  are  in  good 

greeinent  with  those  obtained  in  this  report. 


viu 


Further  tests  will  have  to  be  performed  on  vehicles  of 
different  geometries,  with  the  Inclusion  of  internal  fairings 
and  gates  to  study  the  effects  of  blockage  and  the  like,  before 
a  full  evaluation  can  be  made  of  the  stability  of  an  air- 
breathing  vehicle  traveling  within  a  tube. 
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APPENDIX 


Imagine  the  cylinder  split  into  a  number  of  longitudinal 

strips,  i.e., 


Each  of  the  longitudinal  strips  is  treated  as  a  2>d  flat  plate 
airfoil,  along  which  the  pressure  distribution  is  recorded.  Co: 

sidar  now  two  typical  strips,  chosen  at  the  angles  <p,  /r.id  it  ^ 
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respectively,  such  thft  both  strips  corresponding  to  the  same 
value  of  Z  have  a  common  projection  on  an  imaginary  plane  bi¬ 
secting  the  cylinder  into  its  upper  and  lower  surfaces: 


Because  the  strips  were  chosen  such  that  they  are  both  disposed 
from  the  imaginary  bisecting  plane  by  equal  amounts,  the  loading 
curves  may  be  factored  by  cos  (j)^  and  plotted  together  to  form  a 


closed  curve  from  which  lift  and  momaint  coefflcisr.tb  may  be 


Thus, 

C-  =fc.  Cl*  </(■%) 

c,.  =  c? 

L  =*  t  s 

M.  =  -i-  5  c 


whara  S  s  projoctsd  ares  of  modal  on  .bisecting  plane  ■  cd, 
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Establishing  the  Equivalent  Angle  of  Attack 

fiq^ivalexit  angle  of  attack 
©  S  angle  of  yaw 
(p  S  angle  of  roll 
cK  =.  angle  of  pitch 

At  the  outset  we  have  four  unknowns;  the  procedure  follow¬ 
ed  is  to  stipulate  o<c  and  cj)  .  We  require  two  equations  in  the 
two  unknowns,  0  and  oC  .  Having  stipulated  »<£,  and  cj>  these 
equations  are  solved  for  0  and  o<  and  the  model  adjusted  accord¬ 
ingly. 


From  the  geometry  of  the  figure, 

C2xi  ^  —  R  Ca<. 

Making  a  small  angle  approximation  (cos  o<  5^  1.0)  to 

facilitate  computation  at  no  appreciable  loss  in  accuracy  there 
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results: 


Similarly, 


Rl^o(^J^(l)  =  Rj^c( 


It  la  worthy  of  note  that  the  small  angle  approximation  is 
equivalent  to  saying  that  the  order  of  finite  rotations  may  be 
Interchanged  for  small  angles. 
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SECTION  D 


MEASUREMENT  OF  FORCES  AND  MOMENTS  ACTIInG  ON 
A  VEHICLE  IN  A  TUBE 

by 

John  W.  Burr 


INTRODUCTION 


The  question  as  to  whether  a  vehicle  traveling  in  a  tube 
would  be  aerodynamically  stable  or  unstable  has  been  examined 
in  the  light  of  experimental  pressure  surveys  in  Section  C  of 
this  report.  The  present  section  deals  with  tests  for  the  re¬ 
examination  of  the  same  question  in  the  light  of  direct  measure¬ 
ments  of  forces  and  moments . 

Because  of  a  series  of  failures  in  the  facilities  used  for 
the  tests  reported  here,  the  plan”:ed  experimeni.o!  f^'^eram  could 
not  be  completed.  The  present  experimental  rig,  however.  Is 
satisfactory  and  will  be  used  in  further  tests.  Only  initial 
results  are  presented  here. 

EXPERIMENTAL  APPARATUS 
Facility 

The  tests  reported  here  were  conducted  in  one  of  the.  test 
cells  of  the  Department  of  Aeronautical  Engineering  and  Astro¬ 
nautics  of  Rensselaer  Polytechnic  Institute  Contained  within 
this  cell  is  a  Rolls  Royce  Merl.in  supercharger  driven  by  a 
variable  speed  dynamometer.  The  supercharger  is  connected  to  a 
large  reservoir  (see  Fig.  D-1)  and  can  evacuate  this  reservoir 
to  pressures  ranging  from  atmospheric  to  twenty  inches  of  vacuum, 
depending  upon  the  mass  flow  delivered.  While  running,  the 
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supercharger  can  hold  the  pressure  in  the  reservoir  steady  to 
within  0.2  of  an  inch  of  water. 

Mpdel 

The  tube  of  the  tube-transport  system  is  simulated  by  a 
stationary  6"  diam.  tube  connected  to  the  reservoir  at  its 
discha^S^d  endj  and  with  a  bellmouth  opening  to  the  atmosphere 
at  its  inlet  end  (see  Fig.  D-2).  The  model  of  the  vehicle  is 
held  within  the  tube  with  its  center  about  fifteen  inches  down¬ 
stream  of  the  inlet. 

The  model  of  the  vehicle,  shown  in  Fig.  D-3,  consists  of 
a  tube  seven  inches  long,  2  and  7/8  inches  in  diameter  with  a 
3/16  inch  thick  wall.  The  model  has  interchangeable  leading 
edges . 


The  model  is  supported  at  mid  station  with  a  1/2  inch  thick 
sting  I'hich  is  in  turn  connected  to  a  framework  totaly  supported 
by  strain  beam  members.  Since  the  sting  must  go  through  the  wall 
of  the  tube  without  touching  it,  a  space  must  be  left  around  the 
sting.  This  would  peirrait  leakage  of  atmospheric  air  into  the 
tube  if  a  seal  ware  not  provided.  To  provide  this  seal  the  en¬ 
tire  supporting  structure  <ding  the  strain  beams  is  enclosed 
in  a  chamber  which  was  in  turn  sealed  (see  Fig.  D-4). 
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Ins  t r umenta t ion 

The  strain  beam  members  consist  of  two  drag  beams  and  two 
lift  beams  connected  directly  above  and  below  the  model.  A  third 
drag  beam  is  connected  in  the  horizontal  plane  to  provide  yaw 
stability.  This  beam  is  connected  to  the  framework  by  two  ball 
joints  to  prevent  the  beam  from  providing  any  lift  or  side  motion 
support.  The  force  on  all  three  drag  beams  must  be  added  to  give 
the  total  drag,  the  upper  and  lower  drag  beams  must  be  subtracted 
and  the  result  multiplied  by  the  distance  between  the  pin  joints 
to  give  the  pitching  moment,  and  the  force  on  the  two  lift  beams 
must  be  added  to  give  the  total  lift.  Each  of  the  three  drag 
beams  is  instrumented  with  four  Baldwin  strain  gages  which  are 
connected  into  a  Wheatstone  bridge.  Each  of  the  two  lift  beams 
is  instrumented  with  two  Baldwin  strain  gages  and  these  four 
strain  gages  are  connected  to  form  a  single  Wheatstone  bridge, 
he  output  from  these  four  Wheatstone  bridges  is  read  on  a  Brown 
ferticle-Drum  Millivoltmeter.  The  instrument  rack  containing  the 
lillivoltmeter  and  the  balancing  and  control  resistor  is  shown 
,n  Fig.  D~5. 

■rocedure 


The  system  was  calibrated  by  attaching  a  long  rod  to  the 
lodel  and  hanging  weights  from  the  rod  in  such  a  manner  as  to 


simulate  lift,  drag  and  pitching  moment  and  the  deflection  of 
the  millivoltmeter  was  recorded.  These  readings  were  then  con¬ 
verted  to  calibration  curves.  The  calibration  curve  for  lift 
is  shown  on  Fig.  D-6. 

During  the  running  the  >Jach  number  in  the  tube  was  deter¬ 
mined  from  the  atmospheric  total  pressure  and  the  pressure  in 
the  reservoir  tank.  Since  the  reservoir  was  only  about  four  tube 
diameters  downstream  of  the  model,  the  tank  pressure  was  assuiuid 
to  be  equal  to  the  static  pressure  in  the  tube. 

Discussion 

A  plot  of  the  model  lift  coefficient,  defined  as 

c,  = 

where  A  is  the  total  external  wetted  area  of  the  model,  is  pre¬ 
sented  on  Fig.  D-7.  The  effect  at  Mach  number  is  seen  to  be  smal 
over  the  range  of  the  Mach  numbers  tested.  Also,  changing  the 
shape  of  the  leading  edge  from  a  rounded  edge  to  a  square,  blunt 
edge  does  not  seem  to  affect  the  lift  coefficient  appreciably. 

Comparison  of  these  results  with  the  results  of  the  pressure 
tests  presented  in  Section  C  indicates  good  agreement  (Fig.  7). 
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the  Experimental  Apparatus 


Fig.  D-2.  Installation  of  the  Test  Rig 
(sealing  shroud  removed) 


Fig.  D-3.  Test  Model  showing  Calibration  Rig 

and  Sealing  Chansber 

(model  moved  outside  sealing  chamber  for  calibration) 


Fig.  D~5,  Recording  Millivoltmeter 
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FIG:  D«-6:  Calibration  Curve 
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Section  E 


WALL  EFFECTS  ON  THE  DRAG  OF  A  BODY  MOVING  IN  A  TUBE 

by 


Peter  N.  Szucs 


WALL  EFFECTS  ON  THE  DRAG  OF  A  BODY  MOVING  IN  A  TUBE 


Summary 

Drop  tests  have  been  conducted  to  determine  how  the  drag  of 
a  body  moving  inside  a  tube  is  affected  by  the  tube  walls.  Measure¬ 
ments  have  been  made  for  various  end  conditions,  Reynolds  numbers, 
body-to-duct  diameter  ratios,  and  body  shapes.  Only  steady-state 
conditions  have  been  considered. 

Introduction 

The  effect  of  walls  on  the  motion  of  a  sphere  falling  in  a 
viscous  medium  within  a  vertical  tube  was  studied  for  low  Reynolds 
numbers  by  Ladenburg  (Ref.  1).  The  present  investigation  covers 
higher  Reynolds  numbers  and  a  wider  variety  of  shapes. 

As  in  Ref.  1,  drag  measurements  were  made  on  bodies  falling 
within  vertical  tubes.  The  tubes  were  filled  with  water.  In  some 
of  the  tests  the  lower  end  of  the  tube  was  open  (Fig.  E-la),  where¬ 
as  in  others  it  was  closed  (Fig.  E-lb), 

Apparatus  and  procedure 

The  experimental  apparatus  used  in  these  tests  is  shown  in 
Fig.  E-2.  The  test  tube  was  held  in  a  vertical  position  within  a 
Incite  tank  of  square  (1’  k  I’)  cross  section  and  was  topped  by  an 
electromagnetic  mechanism  for  the  support  and  release  of  the  body. 
The  motion  of  the  falling  body  was  followed  and  recorded  through  an 


array  of  vertically  arranged  light  sources  and  photo-electric 
cells  and  a  Hathaway  recording  oscillograph. 

Four  test  tube  diaaieters  were  used:  1-1/4",  1-3/4",  2"  and 

2-1/2". 

Fig.  E-3  shows  the  bodies  used  in  the  tests.  The  weights  of 
the  streamlined  bodies  could  be  adjusted  by  means  of  mercury 
ballast. 

The  "drag  coefficients"  C^,  defined  as 


(where  F  is  the  drag  force  on  the  body,  Ag  its  maximum  cross 
sectional  area,  p  the  density  of  the  medium,  and  Vg  the  velocity 
of  the  body  relative  to  the  walls),  were  calculated  from  the 
measured  terminal  velocities,  at  which  F  was  equal  to  the  weight 
of  the  body  in  water. 

In  some  tests,  an  aluminum  powder  suspension  was  used  for 
the  purpose  of  visualizing  the  particle  path  lines  In  the  tube- 
fixed  frame  of  reference  (see,  e.g.,  Fig.  E-4). 

Results  and  Discussion 

■iniMi*.  »>■ . »->•■»  ..  ■  •  -  -  — 

The  measured  drag  coefficients  are  plotted  in  Fig-v  E-5  for 
the  spherical  bodies  and  in  Fig.  E-6  for  the  streamlined  bodies. 
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The  range  of  Reynolds  numbers  for  the  conditions  of  these  tests 

lx  4 

is  f>*om  1  X  10  to  4  X  10  . 

It  will  be  noted  that,  although  the  drag  coefficient  is  al¬ 
ways  higher  in  the  closed-end  than  in  the  open-end  condition,  the 
difference  decreases  as  D/d  increases  and  becomes  negligible  for 
D/d  greater  than  about  2.5.  Furthermore,  the  drag  coefficient 
becomes  substantially  independent  of  D/d  (indicating  that  the  wall 
effect  has  substantially  vanished)  when  D/d  is  greater  than  about 
4.  This  result  is  in  sharp  contrast  with  the  low-Reynolds  number 
situation  described  in  Ref.  1,  where  it  is  pointed  out  that,  when 
viscous  effects  are  predominant,  significant  wall  effects  are  ob¬ 
served  even  with  D/d  in  excess  of  100, 

From  Fig.  E-7  it  can  be  seen  that  the  asynjptotic  (free  flight) 
value  of  Cj)  is  approached  more  rapidly  -  hence  that  the  wall  effect 
decreases  more  rapidly  -with  increasing  D/d  when  the  moving  body  is 
a  sphere  than  when  it  is  a  streamlined  bo|^y.  This  result  is  under¬ 
standable,  in  view  of  the  fact  that  viscous  effects  are  more  Im- 
portant  In  the  latter  than  In  the  former  case. 

It  should  be  noted  that  the  drag  coefficient  Is  referred 
here  to  the  velocity  of  the  body  relative  to  the  tube  (Vg)  rather 
than  to  the  velocity  of  the  body  relative  to  the  fluid.  The  latter 
velocity  may  be  taken  to  be  Vg  *  V,  where  V  Is  the  volumetric  flow 
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rate  divided  by  the  cross-sectional  area  of  the  tube.  It  is 
unlikely  that  even  a  drag  coefficient  based  on  the  relative 
velocity  Vg  -  V  will  he  totally  independent  of  the  end  conditions, 
but  this  may  be  assumed  to  be  the  case,  in  first  approximation. 
Then,  letting  subscripts  o  and  c  denote  open-end  and  closed-end 
conditions,  respectively,  one  has,  for  any  given  body  in  any  given 
tube, 

whence 

y  =  1 

v,  '  ICq,) 

The  drag  in  the  open-end  condition  is 

where  is  the  resultant  of  the  frictional  forces  exerted 

by  the  flow  on  the  walls  over  the  entire  length  of  the  tube,  Aj. 
being  the  cross-sectional  area  of  the  tube  and  the  loss  of 

static  pressure  produced  by  these  forces.  It  follows  that 
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For  example,  for  the  conditions  of  the  tests  reported  in 
Fig.  E-5  with  =  2  (i.e..  D/d  =  1.414),  the  drag  coefficients 

Cvv  and  are  2.7  and  1.55,  respectively,  and  Eqs.  (1)  and  (2) 
yield 

V  =  .242  V„ 

D 

AiA^.f=  .9245  F 

Thus,  the  resultant  of  the  frictional  forces  on  the  tube  walls 
is,  under  these  conditions,  equal  to  about  92.5%  of  the  drag  of 
the  body.  As  the  length  of  the  tube  is  increased,  V  decreases, 

Cjj  approaches  ,  and  the  frictional  force  on  the  tube  wall 
increases  asymptotically  to  100%  of  F.  Similarly,  the  fact  that 
Cn  /C-.  approaches  1.0  as  D/d  is  increased  finds  its  obvious  ex- 

UO  Uq 

planation  in  the  decrease  of  V/Vg  that  is  associated  with  increas¬ 
ing  area  ratios  A. /A_. 
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FIG.  E-3  BODIES  USED  IN  TESTS.  SPHERE  DIAMETERS 
(1  to  r):  1",  7/8",  3/16",  3/4",  11/16",  1/2", 
13/32".  MAXIMUM  TEARDROP  DIAMETERS  ARE  1"  FOR 
LARGER  BODY  AND  3/4"  FOR  THE  SMALLER  ONE. 
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PROPULSION  OF  A  VEHICLE  IN  A  TUBE 

by 

Joseph  V.  Foa 


Alio  Issued  as  TR  AE  6404 


SUMMARY 


This  study  concerns  the  flow  disturbances  and  energy 
transformations  that  a  veh3.cle  traveling  in  a  tube  produces 
in  the  adjacent  columns  of  air.  The  analysis  accounts  for 
compressibility,  friction,  and  heat  exchanges.  The  results 
are  applied  to  the  calculation  of  the  powt  r  required  in  rela¬ 
tion  to  blockage  ratio,  vehicle  speed,  an(  mode  of  propulsion. 
At  the  higher  speeds,  the  power  required  is  found  to  depend 
very  critically  on  the  mode  of  propulsion. 


ii 


SYMBOLS 


a  speed  of  sound,  fps 

A  cross-sectional  area  of  the  tube,  sq.  ft. 

K  rnajtimuin  cross-sectional  area  of  the  vehicle,  sq.ft. 

B  =  A^/A  (blockage  ratio) 

Cp  specific  heat  at  constant  pressure,  ft~lb/slug  °R 
c^  nondimensional  friction  factor 
D  = 

F  =  -h  mu.  » (stream  force) 

G  = 

ho  specific  stagnation  enthalpy,  ft-lb/slug 

m  =  pAu  (mass  flow  rate  in  vehicle-fixed  frame  of  reference), 
slugs /sec 

M  Mach  number  in  vehicle-fixed  frame  of  reference 

Mo  -  “i/A 

N  .  /i^(i  -t  rM*} 

p  Static  pressure,  psf 

p°  stagnation  pressure  in  vehicle-fixed  frame  of  reference,  psf 

P  j  power  factor  «  (power  required) /pqA  u^ 

tt 

q  heat  input  per  unit  mass,  ft'lb/slug 

r  radius  of  tube,  ft. 

R  gas  constant,  ft-ib/slug  °R 

T  static  temperature,  °R 


iii 


T°  stagnation  temperature  in  vehicle-fixed  frame  of 
reference,  or 

-t. 

T°  stagnation  temperature  in  tube-fixed  frame  of  reference, 

u  flow  velocity  relative  to  the  vehicle,  fps 

velocity  of  the  vehicle  relative  to  the  tube,  fps 

u^  =  -u^  (velocity  of  the  tube  in  vehicle-fixed  frame  of 
reference) ,  fps 

X  distance  from  vehicle  along  tube,  ft. 

y  =  (u^-u)/u^  (ratio  between  velocity  of  the  flow  and  velocity 
of  the  vehicle,  as  measured  in  tube-fixed  frame  of  reference) 

Z  =  N^D  (See  tabulation  in  the  Appendix,  for  Y  -  1.4) 

jjY  ratio  of  specific  heats 

0  thrust  required,  lb. 

^  density  of  gas  in  the  tube,  slugs/cu.  ft. 

Subscripts  (see  Fig.  1): 

0  conditions  outside  the  tube,  assumed  to  be  uniform 

1  nearest  station  ahead  of  the  vehicle  where  the  flow  can  be 

treated  as  one -dimensional 

2  station  at  tail  end  of  vehicle 

3  terminal  station  of  energy  addition  in  internal  propulsion; 
upstream  end  of  wake 

r\*k 

4  wake  station  downstrcsam  of  which  is  constant 

e  external  propulsion 

i  internal  propulsion 

t  throat  of  fluid  transfer  passage 
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INTRODUCTION 


A  vehicle  may  be  propelled  through  a  tube  in  an  infinite 
variety  of  ways,  but  these  may  be,  grouped,  for  convenience  of 
analysis,  in  three  main  categories.  The  first  category,  which 
includes  conventional  wheel  traction  and  linear  induction  motor 
drives,  is  that  in  which  thrust  is  generated  as  the  reaction  to 
a  force  exerted  by  the  vehicle  on  a  stationary  external  structure. 
Here  the  fluid  in  the  tube  merely  acts  as  a  resisting  medium.  In 
the  second  category  (typified  by  pneumatic  dispatch  systems),  the 
fluid  is  driven  through  the  tube  by  a  pump  or  by  a  high-pressure 
supply,  and  the  vehicle  is  carried  along  by  the  fluid  in  this 
motion.  In  the  third  category,  the  vehicle  still  derives  its  pro¬ 
pulsion  from  displacements  of  the  fluid  columns  within  the  tube 
but  these  displacements  are  produced  by  the  forces  which  are  exert¬ 
ed  on  the  fluid  by  the  vehicle  itself.  The  vehicle  propels  itself 
in  this  case  by  the  fore-to-aft  transfer,  throui;b  itself  and/or 
through  the  space  which  separates  it  from  the  tube  wall,  of  the  air 
or  other  fluid  that  fills  the  remainder  of  the  tube. 

The  present  study  deals  with  propulsion  modes  of  two  of  these 
three  categories;  the  first  and  the  third.  For  lack  of  better  de¬ 
signations,  the  former  will  be  referred  to  as  "external  propulsion" 
and  the  latter  as  "internal  propulsion".  Combinations  of  internal 
and  external  modes  of  propulsion  will  also  be  considered. 


The  situation  studied  here  is  that  of  a  straight,  horizon¬ 
tal  and  infinitely  long  tube,  within  which  a  vehicle  is  traveling 
at  constant  speed.  The  external  air  in  the  vicinity  of  the  tube 
is  at  rest  and  in  uniform  state  throughout,  and  the  tube  is  open 
to  the  surrounding  atmosphere  at  IT'OC  ,  The  present  investiga¬ 
tion  concerns  the  flow  disturbances  that  the  moving  vehicle  pro¬ 
duces  inside  the  tube,  the  associated  losses,  and  the  effect  that 
these  losses  have  on  the  power  required  to  maintain  the  motion  of 
the  vehicle. 

Similar  situations  have  been  treated,  for  incompressible 
fluids,  by  Tollmien  (Hef.  1).  One  of  the  consequences  of  neglect¬ 
ing  compressibility  Is  that,  the  power  available  being  finite,  tu 
calculated  flow  velocities  must  become  everywhere  vanishingly  smal 
as  the  tube  length  considered  is  increased  indefinitely.  In  con¬ 
trast,  when  compressibility  is  taken  into  account,  the  vehicle  is 
found  to  be  capable  of  generating  large  flow  disturbances  over 
long  stretches  of  the  tube  even  when  the  length  of  the  tube  is  in¬ 
finite.  As  a  result  of  these  disturbances  and  of  the  attendant 
thermodynamic  transformations,  each  particle  of  the  air  in  the 
tube  undergoes  a  thermal  cycle  as  it  travels  from  -  oo  to  +  ooin 
tha  vehicle-fixed  frame  of  reference.  With  external  propulsion, 
tha  work  output  in  this  thermal  cycle  is  negative,  whereas  with 
internal  propulsion  the  cycle  is  utilized  for  the  propulsion  of 
the  vehicle  Itself. 
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It  is  precisely  with  such  motions  and  cycles  that  this  study 
is  concerned,  and  special  attention  is  given  to  their  dependence 
on  blockage  ratio,  mode  of  propulsion  and  vehicle  speed,  and  to 
their  effects  on  the  required  power.  The  dependence  of  the  power 
required  on  the  mode  of  propulsion  is  found  to  be  very  critical, 
particularly  at  the  higher  speeds. 

General  governing  equations  for  the  flow  outside  of  region  1-3 . 

In  the  vehicle-fixed  frame  of  reference  the  flow  is  steady. 
Therefore,  tne  governing  equations  for  the  flow  outside  of  region 
1-3  can  be  derived,  in  this  frame  of  reference,  in  the  same  manner 
as  the  familiar  equations  of  steady  one-dimensional  flow  with 
friction,  heat  exchanges  and  chemical  reactions.  The  distinguish¬ 
ing  feature  of  the  present  situation  is  that  the  walls  of  the  tube 
are  moving  at  the  uniform  velocity  in  the  direction  of  the 

relative  flow,  which  is  taken  to  be  the  positive  x  direction.  The 
velocity  is,  of  course,  equal  and  opposite  to  the  velocity  of 
the  vehicle  relative  to  the  ground. 

From  the  definitions  of  Cy  and  F, 

F  ^ pA  (t  + 


Therefore, 


JF.JJSL  Sx.  0!  Jx 

F  I*  r 


(1) 
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From  the  energy  equation 


mc/A*  •  rhJ^ 


and  neglecting  changes  of  specific  heats,  there  follows,  afte 


some  manipulation. 


Jt*  y-/ 

r*  ~ /-. 


(  oi^  I  Jo 

■3^ 


whence 


sl^*m  £r  jf  lal  ,  Jo 

T*  ;7w^  ^  fW  ^ 

The  continuity  equation  may  be  written  in  the  form  (Ref, 


f  -^f  - 

#  ’/If /TT^  7^  +  r  ^ 


Since  A  is  constant 


^^4  -  . 


,  Jf-.. 
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or 


^  j-  dittjr^ 


From  Che  definitions  of  G  and  N 


f.  rn‘f' 


f- 


Therefore . 


A  TV  Li  AAt  - 


Finally, 


(4) 

(5) 


/  </7 
2  7-0 


a 

y- 


/, 


-  l+IrM'^  c/M 

■- Z  /-A)^  ^ 

Eqs.  1,  2,  3,  3',  and  6  yield 


(6) 


4 


/y* 


XvMhiFry, 


A  special  case:  adiabatic  flow. 

The  governing  equations  for  the  adiabatic  case  are  eqs.  1 
through  6'  with  dq  .  0  throughout.  When  u^  •  0,  they  reduce  to 
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the  Fanno  flow  equations . 

Adiabatic  flows  are  of  limited  interest  in  relation  to  the 
problem  at  hand.  The  following  discussion  concerns  only  the 
main  trends  of  the  transformations  in  such  flows. 


Eq.  3'  shows  that 


if  y  is  in  the  range 


Q<  ^  < 


r-t/ 


or,  for  air  (  Y  =1.4), 


<9  <  y  < 


.957 

/  -f- 


When  y  is  within  this  range,  the  local  flow  Mach  number 

o 

(relative  to  the  vehicle)  increases  or  decreases  in  the  direct £■ 
of  the  flow  depending  on  whether  the  local  flow  velocity  (relat 
to  the  vehicle)  is  supersonic  or  subsonic.  The  converse  is  trut 
if  y  falls  outside  of  the  range  indicated  above. 


and 


are,  as  expected,  positive  or  negative  de¬ 


pending  on  whether  y  is  positive  or  negative. 


Eq.  4  reveals  that  the  stagnation  pressure  (relative  to  th 
vehicle)  increases  or  decreases  in  the  direction  of  the  flow, 
regardless  of  whether  the  local  flow  is  subsonic  or  supersonic, 
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depending  on  whether  y  falls  within  the  range 

o  <  j/<  hi- 

or  outside  of  it. 

Finally,  Eq.  6’  shows  that,  if  y  falls  within  the  range 
0  <  y  <  ,  then  the  flow  velocity  increases  or  decreases  in  the 

direction  of  the  flow,  depending  on  whether  the  flow  is  supersonic 
or  subsonic.  The  converse  is  true  when  y  is  outride  this  range. 

Qualitative  information  of  this  kind  is  useful  in  the  quick 
detection  of  solutions  that  are  incompatible  with  the  boundary 
conditions,  as  will  be  shown  in  the  discussion  of  somewhat  differ¬ 
ent  situations  in  the  following  sections. 

Selection  of  an  analytical  model  for  the  present  investigation. 

In  the  application  of  the  general  equations  developed  above 
to  the  study  of  the  flow  in  tube-transport  systems,  special  atten¬ 
tion  must  be  paid  to  the  evaluation  of  the  heat -exchange  deriva¬ 
tives.  These  derivatives  depend  not  only  on  the  thickness  and 
thermal  conductivity  of  the  tube  walls  but  also,  at  each  station 
of  the  tube  and  at  each  instant,  on  the  local  time  derivatives  of 
the  enthalpy  and  velocity  of  the  internal  flow,  hence  on  the  local 
gradients  of  these  quantities  in  the  vehicle-fixed  frame  of  refer¬ 
ence  and  on  the  velocity  of  the  vehicle. 


F-8 


It  is  obvious,  however,  that  the  actual  program  of  heat 
exchange  between  the  internal  flow  and  the  surrofa’'iing  atmo&phe 
must  lie  between  two  extremes.  At  one  extreme  is  the  case  of  ze 
conductivity  -  the  adiabatic  flow.  This  case  is  of  no  more  th-^ 
quantitative  interest  in  relation  to  the  problem  at  bind.  At  t 
other  extreme  is  the  case  of  infinite  conductivity  whereby  thf 
internal  boundary  layer  is  kept  ever3rwher£  at  the  temperature  o 
the  extern:sl  air.  For  moderate  velocities  of  the  vehicle,  the 
latter  extreme  is  believed  to  be  very  olo^e  to  the  sttu^l  sicurt 
tion. 

Ifeder  the  conditions  that  are  to  be  coiisiderec  In  this  pi>.. 
leraj  the  flew  is  fully  turbulent  in  those  regions  where  it  un”^- 
goes  its  most  significant  transformations.  Furthermore,  as  wll 
be  shown  presently,  all  these  transformations,  thr<>ughout  misi- 
the  subsonic  range  of  vehicular  speeds,  take  place  very  gr^Ui  m' 
and  over  great  distances,  both  ahead  of  the  vehicle  ann  behind  i 
except  in  its  immediate  proximity.  These  observacion<-  suggest 
simple  analytical  model,  whereby  the  tube  wall  is  viewed  as  act* 
like  a  very  large  heat  exchanger,  through  which  the  beat  flows 
such  a  manner  that  the  stagnation  temperature  of  the  flow 
relative  to  the  tube  wall  is  made  to  remain  ever3rwh£'re  (except 
around  and  immediately  behind  the  vehicle)  substantially  equal  t 
the  external  temperature  Tq. 
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A  more  exact,  but  necessarily  only  numerical,  analysis  of 
the  same  problem  is  carried  out  in  Ref.  3,  where  Eqs.  1  through 
6  of  the  present  paper  are  applied  with  the  stipulation  that 

^  •  -Kir*-  Tj 

where  K  is  a  constant.  The  analytical  model  proposed  here,  al¬ 
though  less  exact,  is  believed  to  describe  the  actual  physical 
situation  in  fair  approximation  over  a  wide  range  of  speeds.  It 
also  has  the  advantage  that  it  leads  to  closed-form  solutions 
which  lend  themselves  to  clear  and  informative  interpretation. 
Nevertheless,  its  plausibility  is  limited,  as  has  been  noted 
above,  to  those  situations  in  which  the  tube  can  be  viewed  as 
performing  the  function  of  a  "very  large"  heat  exchanger  -  very 
large,  that  is,  in  relation  to  the  heat  transfer  rates  that  it  is 
called  upon  to  provide.  As  the  speed  of  the  vehicle  is  increased, 
the  rate  at  which  energy  is  dissipated  in  the  internal  flow  in¬ 
creases  and  the  gradients  of  all  flow  properties  ahead  of  the 
vehicle  become  steeper  and  steeper.  At  very  high  speeds,  the 
constant-T*^*  condition  cannot  be  maintained,  and  the  path  of  the 
transformation  must  approach  that  which  is  finally  established 
when  all  the  flow  changes  ahead  of  the  vehicle  take  place  through 
one  or  more  shocks  in  the  neighborhood  of  station  1.  The  condition 
that  is  approached  is  that  of  an  isoenergetic  process  in  the  vehicle 


fixed  frame  of  reference. 
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The  flow  trantfonaatlon  ahead  of  the  vehicle  will,  accord. 
Ingly,  be  analysed  aa  a  constant-T°*  procesa  for  xooderate  aub- 
sonlc  speeds,  and  as  a  constant>T°  process  for  supersonic  speed 
In  the  transonic  range,  the  flow  process  may  be  expected  to 
follow  an  Intermediate  path. 

Behind  the  vehicle,  the  condition  T®*  «  constant  is  tenabl 
under  all  conditions  downstream  of  a  station  A  which  is  Itself 
some  distance  -  as  yet  unspecified  >  downstream  of  station  3 
(Fig.  1).  The  region  between  stations  1  and  3  is  one  in  which 
rapid  transformations  take  place  and  energy  Is  added  to  the  fir 
in  some  modes  of  Internal  propulsion.  The  distance  from  3  to  £■ 
is  the  distance  required  for  enough  heat  to  be  extracted  from  t 
flow  to  bring  T®*  to  a  aubatantially  etationary  value. 

The  flow  ahead  of  the  vehicle. 

The  conetent-T°*  condition  will  be  anelysed  first. 

From  the  definitions  of  T®*  and  T*^  there  follow*? 

^  T*'^ (f- 0 

whence,  if  T®*  *  To* 
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f, 

I 

1 

I 

I 


and 


7"°-  Z 


(9) 


Eq.  9  reveals  the  important  (and  potentially  useful)  fact 
that  the  stagnation  temperature  of  the  internal  flow,  relative 
to  the  vehicle,  is  lower  or  higher  than  the  stagnation  tempera¬ 
ture  that  would  be  measured  in  the  free  stream  at  the  same  IJach 
number  Mq,  depending  on  whether  u/u^  is  smaller  or  greater  than 
1.0  (i.e.,  depending  on  whether  y  is  positive  or  negative). 

When  y  a  .5,  the  internal  stagnation  temperature  relative  to  the 
vehicle  is  equal  to  the  external  static  temperature.  This  seem¬ 
ing  paradox  is  explained  by  the  above-noted  heat  exchanger  action 
of  the  tube. 


From  Eqs.  8  and  9, 


whence 


^ 

Finally,  Eqp.  1,  6,  11,  and  tb i  continuity  equation 

1 


(11) 


fi/  Z  T' 


F' 


yield 


0-^U!fl  f 


rc^  JUj 


Therefore,  when  consideration  is  given  to  a  flow  region  in 
which  changes  of  can  he  neglected, 


C12) 


Eq.  12  permits  the  calculation  of  the  change  of  velocity 
from  a  station  a  to  a  station  ]b  in  any  region  upstream  of  station 
1  or  downstream  of  station  4.  On  the  right-hand  side  of  the  equa 
tion,  the  ” -f*  ”  sign  applies  when  y  >  0  and  the  sign  when 

y  <  0. 


If  the  Reynolds  number  is  found  to  undergo  large  changes 
within  the  region  considered,  the  calculation  should  be  repeated 
on  "ubregions  small  enough  that  the  changes  of  in  each  may  be 
neglected  or  adequately  accounted  for  by  the- use  of  suitable  aver 

ages. 
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The  stream  force  can  be  calculated  as  a  function  of  y  from 
the  continuity  equation  in  the  form 


F^N  constant 


which,  together  with  Eqs.  9  and  11,  yields 


Fit -a) 

tL  ~  td-nf 


-  constant 


At  ^  Oo  ,  y  =  0  and  F  =  PjA(l  -ff*  I'Mq  ).  Therefore, 

^ 


(13) 


As  has  been  noted  earlier,  the  assumption  of  a  constant  T° 
becomes  untenaole  at  speeds  approaching  the  speed  of  sound  in 
the  atmosphere  surrounding  the  cube,  and  at  supersonic  speeds  it 
becomes  more  plausible  to  treat  the  flow  ahead  of  the  vehicle  as 
isoenergetic  in  the  vehicle-fixed  frame  of  reference.  When  all 
flow  changes  take  place  through  shocks  around  the  nose  of  the 
vehicle,  conditions  at  station  1  are  obviously  the  same  as  at  in¬ 
finity,  and  the  flow  from  station  1  to  the  throat  can  be  treated 
like  the  inflow  through  a  supersonic  intake.  If,  on  the  other 
hand,  changes  occur  ahead  of  station  1,  then  the  flow  at  1  is 
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subsonic,  M-j^  is  determined  by  conditions  further  downstream,  and 
the  continuity  equation  gives 

A  _ 


M.  i,  .  WJ* 
H 


(15) 


An  examnle. 


To  illustrate  the  application  of  Eqs.  9  through  14,  consides 


the  flow  ahead  of  a  vehicle  traveling  at  550  fps  within  a  tube 
surrounded  by  air  at  Po  =  2000  psf  and  at  Tq  =  500  °R.  At  this 
temperature  the  speed  of  sound  is  a^  ■  1100  fps.  Therefore, 

Mq  =  .5.  Let  A  =  100  sq.  ft.,  and  let  the  blockage  ratio  and  th« 
mode  of  propulsion  be  such  that  y  »  .5  just  ahead  of  the  vehicle 
With  the  origin  (x  =  0)  at  station  1,  the  equations  above  provi;?* 
the  following  description  of  the  flow; 

TA9LE  X 

average  distai 
ahead 


y 

a 

fps 

F 

PqA 

M 

Po 

Po 

T© 

^0 

xlO-^ 

vehic 
mi  lei 

0 

550 

1.350 

.500 

1.00 

1.186 

1,05 

.1 

495 

1.425 

.450 

1.11 

1.275 

1.04 

2.25 

24 

.2 

440 

1.528 

.401 

1.249 

1.395 

1.03 

2.0 

11 

.3 

385 

1.667 

.351 

1.442 

1.55 

1.02 

1.8 

6 

.4 

330 

1.860 

.301 

1.651 

1.76 

1.01 

1.75 

2 

.5 

275 

2.145 

.252 

1.971 

2.065 

1.0 

0 
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The  calculated  Mach  number  just  ahead  of  the  vehicle, 

-  .252,  is  the  Mach  number  that  would  be  maintained  by  a 
well  streamlined  body  of  blockage  ratio  5— .58,  i.e.,  by  a 
vehicle  of  maximum  diameter  equal  to  about  76/i  of  the  diameter 
of  the  tube. 

The  stagnation  pressure  ahead  of  the  vehicle  is  seen  to  be 
considerably  higher  than  the  free-=stream  stagnation  pressure  at 
the  same  Mach  number.  On  the  other  hand,  the  stagnation  tempera¬ 
ture  relative  to  the  vehicle  is  everywhere  lower  than  the  stagna¬ 
tion  tenperature  that  would  be  measured  in  the  free  stream  at  the 
same  Mach  number,  in  accordance  with  the  general  observations 
made  in  the  discussion  of  Eq.  9. 

Inspection  of  Eq.  12  and  of  Table  I  reveals  a  strange  diffi¬ 
culty:  if  one  were  to  carry  out  a  similar  calculation  for  positive 
values  of  x,  starting,  e.g.,  with  y  a  .5  at  station  4,  y  would  be 
found  to  increase  monotonically ,  approaching  1.0  asymptotically j 
as  X  increases  indefinitely.  This  result  would  indicate  that  the 
region  behind  the  vehicle  is  an  infinite  wake  in  which  the  flow 
velocity  relative  to  the  vehicle  continues  to  decrease  with  in¬ 
creasing  distance  from  the  vehicle.  Such  a  wake  would  clearly  be 
incompatible  with  the  boundary  condition  «  0*  becomes 

necessary,  therefore,  to  re-examine  the  solutions  in  order  to 
determine  which  flows  are  possible  in  the  tegion  behind  the  vehicle 
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and  which  are  not. 

The  flow  behind  the  vehicle. 

Because  of  the  boundary  condition  only 

possible  flows  behind  the  vehicle  are  those  in  which  dy  ^  0 

dx 

for  y  ^  0,  respectively. 


Eq.  12  yields 

c/& 


where  the  "  +  "  sign  applies  when  y  >  0  and  th^*,  ~  sign  when 

y  0. 


The  numerator  of  the  right-hand  member  of  Eq.  16  is  alwcy 
positive. 


Now,  if  y  >0,  the  "  ”  sign  applies  in  F.q»  16.  On  the 

other  hand,  where  y  >  0,  the  stability  of  the  wake  demands  that 
dy/dx  be  negative.  Therefore,  a  wake  with  y  >  0  is  compatible 
with  the  boundary  conditions  only  if  the  denominator  of  the 
right  hand  member  of  Eq.  16  is  negative,  i.e.,  if 


This  condition  must  be  satisfied  everywhere,  including 
station  4,  where  y  will  have  a  higher  value  than  at  any  statioi 


further  downstream. 
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Where  y  <0,  the  sign  applies  in  Eq.  16  but  the  flow  is 

stable  only  if  ;>.0.  Therefore,  the  compatibility  condition 
dx 

requires  again  that  the  denominator  of  the  right-hand  member  of 
Eq.  16  be  negative.  In  this  case,  however,  the  highest  valfae  of 
y  anywhere  downstream  of  station  4  is  zero.  It  follows  then  from 
Eq.  17  that,  where  y  ^0,  the  flow  is  compatible  with  the  bound¬ 


ary  conditions  wherever  > 


/ 


VT 


.  The  line  of  demarcation  be¬ 


tween  stable  and  unstable  wake  conditions  at  station  4  is  plotted 
as  curve  a  in  Fig.  2.  It  should  be  noted  that  the  condition 


z  0  represents  a  singularity,  because  it  is  a  condition  that 
is  compatible  with  the  boundary  conditions  for  all  values  of  Mq: 
when  y4  3  0,  the  flow  is  at  rest  relative  to  the  tube  everywhere 
downstream  of  station  4. 


Curve  b  in  Fig.  2  is  a  plot  of  the  equation 

/V.*- 

which  represents  the  vanishing  of  the  denomination  of  the  right- 
hand  member  of  Eq.  16.  If  y^  and  Mq  are  such  that  the  initial 
point  lies  in  the  negative  quadrant  below  this  curve  (like  point 
P  in  the  figure),  y  will  cl  ange  downstream  of  station  4  in  the 
manner  indicated  by  the  arrow  until  curve  ^  is  reached.  At  that 
point  y  becomes  discontinuous,  and  no  solution  is  found  for  the 
flow  downstream  of  the  station  at  which  this  discontinuity  occurs. 
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The  only  physical  meaning  that  can  be  attributed  to  the 
"unstable  wakes"  is  that  they  cannot  exist  except  as  nonsteady 
transients  through  which  the  downstream  boundary  conditions 
eventually  impose  on  the  flow  the  condition  y^  =  0.  Therefore, 
the  only  steady-flow  condition  that  is  possible  behind  the 
vehicles,  when  Mq  is  lower  than  the  critical  value  for  wake 
stability,  is  one  in  which  the  flow  is  at  rest  (relative  to  the 
tube)  everywhere  downstream  of  station  4. 

This  result  is  of  great  importance  in  the  determination  of 
the  drag  of  the  vehicle  and  of  the  power  required  to  maintain 
its  motion  with  various  means  of  propulsion. 

According  to  Eq.  16,  downstream  of  station  4  the  derivati''  • 
dy/dx  can  only  vanish  as3miptotically  or  be  identically  zero 
everywhere.  Therefore,  y  in  this  region  can  only  vary  monotoni- 
cally  with  x  or  remain  constant. 

The  flow  transformation  from  station  3  to  station  4  can  be 
analyzed  as  follows; 

The  energy  equation  between  these  two  stations,  in  the 
vehicle-fixed  frame  of  reference,  may  be  written  in  the  form 

Tj)  =  ^m(F^  ^^3-4 


(1 
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where  ^3,4  is  the  heat  received  by  the  unit  mass  of  the  fluid 
between  the  two  stations: 

<1.-*  =  (rr-  rr) 

~  <T (To- 

Equs.  7  and  9  may  be  re-written  as 

(7-) 

and 

2' <s') 

Furthermore, 

(19) 


and,  from  Eq»  13, 


These  equations  yield 


A 


/ 


i.e.,  (by  virtue  of  Eq.  14),  P3  =  P4. 


(13‘) 


(20) 


Finally,  from  Eq.  20  and  the  continuity  equation  in  the 


form 
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•3 


one  obtains 


or 


/ 


[f  -  -7^0- 


(21? 


When  74-0  (i.e.,  when  is  below  the  critical  value  for 
stability  of  the  wake),  P3  -  P4  -  Po  and 

t  “  /-y.  “  ^ ^-/) 

-f  /hu^O-t/,)  ^23) 

•  ^  ^  ^  -/J  (2^) 

All  other  flow  properties  at  station  3  are  then  readily 
calculated  for  this  case.  If,  on  the  other  hand,  Mq  is  above 
the  critical  value,  the  wake  is  stable  and  74  is  no  longer 
uniquely  determined  by  the  boundary  conditions.  The  flow  at 
station  3  can  then  be  specified  independently  of  these  conditions 


these  results  will  now  be  applied  to  calculation  of  the  power 
required  to  maintain  the  steady  level  motion  of  simple  vehicles  by 


c 
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various  means  and  under  various  conditions. 


Power  required. 

The  following  calculations  concern  only  the  power  required 
for  propulsion  and  do  not  account  for  the  drag  of  supporting  and 
stabilizing  devices;  nor  do  they  account  for  the  possible  utili¬ 
zation  of  boundary  layer  air  in  internal  modes  of  propulsion,  or 
for  the  efficiencies  of  the  components  of  machinery  through  which 
the  required  power  is  delivered  from  its  prime  source  to  the 
vehicle  or  to  the  flow.  The  vehicle  is  assumed  to  be  in  each 
case  a  body  of  revolution,  with  its  axis  in  line  wi'th  the  axis 


of  the  tube. 


The  various  modes  of  propulsion  are  evaluated  in  terms  of  a 
nondimens ional  power  factor  P,  which  is  defined  as  the  ratio  be¬ 
tween  the  power  required  and  the  quantity  PqAu^,. 

In  the  case  of  external  propulsion,  the  power  required  is 
the  product  of  the  drag  and  the  velocity  of  the  vehicle.  Thus, 

P  ■  s  ^  Z/^  (25 

When  internal  propulsion  is  considered,  the  power  required 
is  calculated  as  the  increment  of  stagnation  enthalpy  (relative 
to  the  vehicle)  that  must  be  Impacted  to  the  flow  in  the  unit  of 
time  in  otder  to  maintain  the  steady  motion  of  the  vehicle.  The 
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transfer  of  energy  to  the  flow  is  assumed  to  take  place  in  ev 
case  between  stations  1  and  3.  Thus, 

p.  . 


A 

t-f  Tc 


The  following  conditions  and  inodes  of  propulsion  will  be 
considered: 


(1)  Mq  a  0.8,  with 

(I, a)  external  propulsion,  and 

internal  propulsion,  by  means  of  a  far. 
ahead  of  station  t. 

(2)  My  =  2.5,  with 

(2, a)  external  propulsion, 

(2,b)  internal  propulsion  by  means  of  a  crypto- 
steady  energy  exchanger  ahead  of  station 
t ,  and 

(2,c)  internal,  ramjet-type  propulsion 

Condition  1.  Mg  ^  0.8.  The  diameter  of  the  vehicle  consider 
Is  equal  to  68.5%  of  the  diameter  of  the  tube.  Thus,  B  =  .46 
and  A^/A  ■  .531.  The  relative  flow  is  sonic  at  station  t.  Th 
being  a  borderline  situation,  the  power  factor  will  be  calcul 
first  for  a  constant-T°*  flow  and  then  also  for  a  constant-T° 
flow  ahead  of  the  vehicle. 
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With  external  propulsion,  the  relative  flow  is  iso- 
energetic  from  station  1  to  station  3.  Assuming  the  flow  to  be 
isentropic  from  station  1  to  station  6,  the  continuity  equation 
gives  Dj^  -  .531D(1)  -  .3077,  hence  >  .328.  Then,  Eqs„  11„  13, 
and  9  yield 

»  .6' 

^  'j'O 

-  .37^-f 

'o 

Since  the  flow  between  1  and  3  is  isoenergetic , 

Furthermore,  Mq  is  subcritical.  Therefore,  P3  =  Po,  and  Eqs.  22 

2 

and  24  give  y^  =  ,125  and  M3  =  .56.  Finally, 

and 


/O  .  -.955 

It  should  be  noted  that  in  this  condition  the  flow  receives 
no  energy  and  is  not  otherwise  disturbed  between  stations  2  and 
3.  Therefore,  station  3  effectively  coincides  with  station  2. 
Since  the  cross-sectional  area  at  2  is  the  same  as  ac  1,  the 
stagnation  pressure  ratio  is  obtained  from  the  continuity  equation 
for  the  isoenergetic  flow  between  these  two  stations  as 


/=?  P, 


» .  566 
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The  large  loss  of  stagnation  pressure  is  due  primarily  to 
a  shock  between  stations  t  end  2. 


For  the  condition  of  internal  propulsion,  the  thrust 
generator  considered  in  this  case  is  a  fan  located  ahead  of 
station  t.  With  y  =  1.4,  and  neglecting  the  entropy  incremer. 
through  the  fan,  the  continuity  equation  between  gtatioi^.s  1  ar^ 
t  may  be  written  as 

Furthermore,  P2  =  P3  =.•  Po  because,  as  before,  Mq  is  sub- 
critical  and  stations  2  and  3  coincide.  Also,  since  the  flow 
is  isoenergetic  from  t  to  3,  -  ^2°  •  Finally,  the  net 

horizontal  force  acting  on  the  vehicle  must  be  zero,  hence 

these  conditions,  Eqs .  9,  11,  13,  14,  22, 

24  <nd  27  yield; 


y/  -  .235 
*  .615 

r, 

% 

Tzy 


1.0679 


-1.25 


y?**  .95 
/^*.84 
1.984 
/£i-.533 


It  ts  interesting  to  note  that,  whereas  the  total  pressure 
loss  frrii  statiott  t  <0  station  2  is  higher  in  this  case  than 
with  external  propul e Ion,  Lhht  over/ilJ  entropy  rise  between 
st;itions  X  and  2  is  aefeuaily  307.  lower. 
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Finally, 

Pc 

Thus,  the  ratio  between  the  power  required  with  this  mode 
of  internal  propulsion  and  the  power  required,  for  the  same 
vehicle  and  at  the  same  Mach  number  with  external  propulsion, 
is  Pi/Pg  =  .668, 

If  '1q  =  500  ‘^R,  hence  =  1100  fps^  the  velocity  of  the 
vehicle  in  the  condition  just  considered  is  880  fps.  The  maxi¬ 
mum  velocity  of  the  air  ahead  of  the  vehicle,  relative  to  the 
tube,  is  529  fps  with  external  propulsion  and  207  fps  with  in¬ 
ternal  propulsion.  The  maximum  velocity  of  the  air  in  the  wake, 
relative  to  the  tube,  is  110  fps  with  external  propulsion  and 
83,5  fps  (in  the  direction  opposite  to  the  direction  of  motion 
of  the  vehicle)  with  internal  propulsion. 

Intermediate  situations,  involving  both  internal  and  ex¬ 
ternal  propulsion,  may  be  investigated  as  follows: 


For  various  values  of  yj^,  Eq.  11  provides  hence 
Then,  Eqs.  13,  27,  and  9  provide  the  corresponding  values  of 


and 


z 


j  Eq.  24  gives  M^”,  and  F3/P0A  is 
To 


obtained  as  1  -f  •  The  difference  between  and  F3  is  the 


thrust  that  must  be  provided  by  external  forces  to  maintain  the 
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steady  inotion  of  the  vehicle.  The  results  are  plotted,  for  tt 
coridition  under  consideration,  in  Fig.  3.  Values  of  inter¬ 
mediate  between  .235  and  .6  are  produced  when  propulsion  is  oj 
mix^d  type,  i.e.,  both  internal  and  external,  the  latter  compc 
of  the  thrust  being  Fj  -  F3.  The  power  factor  P  decreases  me: 
tonically  as  the  mode  of  propulsion  is  changed  from  purely  ex> 
(yx  =  .6)  to  purely  internal  (y^  =  .235;.  The  condition  yi 
occurs  when  the  internal  thrust  exceeds  the  value  required  to 
tain  steady  motion. 

If  the  constant-T'-*  condition  is  replaced  by  the  assumpt’ 
of  an  isoenerg  ic  flew  ahead  of  station  1,  the  following  re>> 
are  obtained. 

With  external  propulsion,  =  *328  as  before,  but  Eq.  1! 
yields 

Furthermore,  since  13°  «  q  Tq®  3  1.128  Tq,  Eqs.  22  and  24 
m  .64,  hence  (remembering  that  P3  «  p^,  because  Mq  is  sub^ 
critical) 

and 
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For  the  case  of  internal  propulsion,  the  condition  =  F2^ 
and  Eqs.  15,  17,  and  24  yield  =  .621,  13°  >  1.1745  Tq  M3^  -  .739, 
and  Fj^/PqA  =  F3/PQA  =  2.03.  Then,  since  >  1.128  Tq, 

Thus ,  the  modified  analytical  model  yields  P^/Pe  =  • 662 . 

This  result  is  very  close  to  that  obtained  through  the  use  of  the 
constant-T°*  assumption  (P^/Pe  »  •668).  It  appears,  therefore, 
that  the  latter  assumption  can  be  used  with  confidence  for  all 
subsonic  speeds  up  to  at  least  =  .3.  However,  the  calculated 
power  required  is,  at  this  Mach  number,  between  8  and  10  percent 
higher  when  the  flow  is  assumed  to  be  isoenergetic  than  when  the 
constant-T°  assunqjtion  is  used. 

Condition  2.  Mq  =  2.5.  The  diameter  of  the  vehicle  is  72.%  of 
the  inner  diameter  of  the  tube.  The  nose  of  the  vehicle  is  de¬ 
signed  to  produce,  through  a  train  of  oblique  shocks,  a  stagnation 
pressure  ratio  (Pt.°/P(^°)  of  *80  at  Mq  =  2.5  in  the  absence  of  flow 
disturbances  ahead  of  station  1.  The  stipulated  diameter  ratio 
corresponds  to  B  =  .525  and  Aj./A  •  .475. 

With  external  propulsion,  the  flow  is  choked  at  t  and  sub¬ 
sonic  at  1.  If  the  flow  in  this  condition  is  assumed  to  be  isen- 
tropic  between  1  and  t,  the  continuity  equation  gives 


F-28 


Dj  -  .475  (Dl)  «  .275,  hence  .  .28?  and  .  .2609.  Then 
Eq.  15  yields  Fj^/pqA  «  14.39. 

Assuming  P2°/Pi°  ■  *95,  one  obtains 

D2  =  D]^/.95  =  .29,  hence  M2  =  2.194,  N2  =  .3972 

and  finally 

/?  -  •  ^-55 

The  possibility  that  the  flow  may  be  subsonic  at  2  must  '!■' 
excluded.  In  fact,  with  the  supersonic  solution  (M2  s  2.194) 
selected  above. 

If  the  flow  at  2  were  subsonic,  the  pressure  there  would  1 
at  least  five  times  higher  and  would  increase  further  from  the: 
In  the  direction  of  the  flow.  This  situation  would  be  incompa' 
ible  with  the  boundary  conditions.  Therefore,  the  flow  at  2  it 
supersonic  and  underexpanded.  It  is  possible  that  the  flow  wi. 
adjust  itself  to  a  lowei*  pressure  further  downstream  through  a 
pattern  of  expansion  waves  between  stations  2  and  3. 

With  Internal  propulsion,  if  a  cryptosteady  energy  exchang 
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is  used  to  energize  the  flow  between  1  and  t,  the  disturbances 
ahead  of  1  can  be  ’'swallowed”.  Then  the  flow  at  1  is  undisturbed, 
and  Fi  =  Fo,  =  M^,  and  =  Tq®. 


The  input  of  energy  Is  assumed  to- take  place  is ent topically 
behind  the  oblique  shocks  surrounding  the  nose,  through  which  the 
stagnation  pressure  ratio  retains  the  stipulated  value  .8.  A 
further  5%  loss  of  stagnation  pressure  is  assumed  to  take  place 


_  Pi^D,  r  p:a,\ 

^  (  W'J 

becomes 


On  the  other  hand,  the  momentum  equation  (F2  =  Fp)  may  be 
written  in  the  form 


Eqs^  28  and  29  yield 


76  =  Zo  =  7l0£»/O-‘ 
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whence 


Z2  -  535  /O' 


*  2.3ZO 


.33/7 


Therefore,  ^2  ^  <.  3917/. 3846)^  =  1.03  73,  and  finally 


7=r  tt 


The  ratio  between  the  power  required  with  this  mode  of 

■■  Af 

ternal  propulsion  and  the  power  required,  for  the  same  vehicL* 
at  the  same  Mach  number  with  external  propulsion,  is  in  fcj- 
case  P;j^/Pe  =  •0594.. 


It  can  be  shown  that  the  flow  is  always  underexpanded  at 
with  internal  propulsion  if  the  supersonic  solution  for  M2 
lower  than  Mq.  Indeed,  from  the  condition  F2  s  F3,  one  obtair 


This  confirms  that  in  the  example  just  treated  no  shock 
occurs  between  stations  t  and  2.  However,  as  before,  it  is  I3 
that  expansion  waves  will  be  present  between  stations  2  and  3, 
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These  waves  will  not  affect  the  power  required. 

If  the  mode  of  propulsion  is  of  the  ramjet  type,  the  required 
power  factor  is  somewhat  higher,  although  still  one  order  of  magni¬ 
tude  lower  than  with  external  propulsion.  It  should  be  noted  that 
the  use  of  purely  thermal  propulsion  in  the  type  of  vehicle  under 
consideration  would  not  be  practical,  not  only  because  it  would 
require  the  assistance  of  auxiliary  devices  for  acceleration  to 
the  design  speed,  but  also  because  it  would  not  provide  for  sta¬ 
bility  of  the  flow  ahead  of  the  throat.  This  stability,  as  well 
as  some  of  the  thrust  required  at  the  lower  speeds,  could  be  pro¬ 
vided  by  a  fan  or  by  an  ejector  pump  ahead  of  station  t.  In  the 
present  analysis  of  cruising  operation  at  Mq  «  2.5  it  is  assumed, 
for  simplicity,  that  whereas  most  of  the  energy  Input  is  in  the 
form  of  heat  addition  to  the  flow  between  stations  t  and  2,  the 
work  performed  by  the  pumping  device  ahead  of  t  is  just  enough  to 
make  up  for  the  loss  of  stagnation  pressure  between  t  and  2,  The 
calculated  power  factor  will,  of  course,  account  for  the  total 
power  that  is  required  for  this  combined  operation. 

Under  these  conditions,  the  loss  of  stagnation  pressure  from 
1  to  2  is  again  207«,  and  the  momentum  equation  (^2  =  con¬ 
tinuity  equation  (p  constant)  become 

.s(Sg  -  Go 
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whence  M2  =  2.223  and  12° /Tq^  x  1.0601.  Then,  since  To°/To  = 
there  follows  T2°/Tq  x  2.385  and 

This  power  factor  is  less  than  one  tenth  of  the  power  f-‘:r 
for  external  propulsion.  It  should  be  noted,  however,  that  th 
calculated  value  of  for  thermal  propulsion  is  very  sens  it iv 
changes  in  the  losses  associated  with  heat  addition  and  in  the 
tribution  of  mechanical  and  thermal  input. components.  Even  a 
moderate  increase  of  stagnation  pressure  loss  downstream  of  f 
not  conpensated  through  an  increased  allocation  of  power  to  th 
mechanical  input  device  ahead  of  t,  may  result  in  a  substantia 
increase  of  power  required. 


The  input  of  thermal  energy  may  be  effected  either  through 
combustion  of  a  fuel  or  by  transfer  of  heat  from  a  suitable  so', 
(e.g.,  a  nuclear  reactor).  The  increase  of  stagnation  tempera' 
that  is  required  for  thermal  propulsion  is  relatively  small  ant 
should  be  obtainable  with  heat  exchangers  of  acceptable  siae  ar 
drag.  In  some  cases  it  may  be  advisable  to  carry  out  the  addit 
of  heat  subsonically  (in  order  to  reduce  the  attendant  loss  of 


stagnation  pressure)  and  to  re-ac celerate  the  flow  to  supersonic 
speed  through  a  converging-diverging  passage  downstream  of  the 
heating  region. 

Internal,  ramjet-type  propulsion  with  heat  addition  behind  the 
vehicle. 

Under  certain  conditions  it  may  be  desirable  to  add  heat  to 
the  flow  between  stations  2  and  3  rather  than  upstream  of  station 
2.  When  this  is  done,  thrust  is  generated  in  much  the  same  manner 
as  it  is  in  most  conventional  ramjets,  except  for  the  fact  that, 
in  the  case  of  the  vehicle  in  the  tube,  the  structure  of  the 
"combustion  chamber"  is  part  of  the  stationary  tube  and  not  of  the 
moving  vehicle.  Of  course,  no  flow  process  downstream  of  station 
2  can  have  any  effect  on  the  vehicle  unless  the  flow  at  this 
station  -  and  for  some  distance  upstream  of  it  -  is  subsonic. 

It  has  been  shown  that  when  the  flow  at  2  is  subsonic  the 
pressure  there  is  too  high  to  be  supported  by  a  stable  wake  under 
the  conditions  examined  in  the  preceding  examples.  This  diffi¬ 
culty  is  eliminated,  however,  when  the  subsonic  flow  condition  at 
2  is  sustained  by  a  Rayleigh  process  between  2  and  3,  because  the 
rate  of  heat  input  that  is  capable  of  maintaining  this  condition 
is  always  above  the  critical  value  that  produces  thermal  choking 
at  the  downstream  end  of  the  heating  region.  With  M3  =  1,  the 
upper  limit  of  the  pressure  upstream  of  station  3  is  no  longer 
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controlled  by  the  downstream  boundary  conditions. 


The  calculation  of  the  power  factor  is  particularly  sim 
with  this  mode  of  propulsion  because,  from  the  conditions 
F3  =  F2  =  Fq  and  M3  =  1,  one  obtains 


whence 


For  Mq  =  2.5  in  air,  =  3.2  7.  This  power  factor  is  v’ 
above  the  value  calculated  above  for  the  case  of  heat  additi 
tween  t  and  2,  mainly  because  of  the  large  entropy  ir\prement 
is  produced  by  the  shock  through  which  the  flow  becomes  subi- 
upstream  of  2.  Nevertheless,  this  mode  of  propulsion  has  se* 
attractive  features  that  make  it  worthy  of  consideration  desj 
its  higher  power  requirements.  Its  most  important  advantage 
that  it  permits  the  input  of  propulsive  energy,  in  the  form  ( 
heat,  from  sources  which  are  activated  by  the  passing  vehicle 
fixed  to  the  tube. 


Discussion  of  results,  and  concluding  remarks . 

The  simple  analysis  developed  in  this  study  reveals  that 
flow  disturbances  that  are  set  up  by  a  moving,  vehicle  in  a  ti 
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are  largely  determined  by  the  mode  of  propulsion.  The  same  can 
be  said,  of  course,  of  the  attendant  flow  losses  and  of  the  power 
required.  In  other  words,  the  thrust  that  is  required  to  main¬ 
tain  the  motion  of  a  vehicle  in  a  tube  depends  on  the  manner  in 
which  the  thrust  itself  is  generated.  The  dependence  is  parti¬ 
cularly  critical  when  the  speed  of  uhe  vehicle  is  high.  This 
result  is  illustrated  by  the  relative  magnitudes  of  the  power 
factors  calculated  in  the  examples  above  for  a  few,  randomly 
selected  modes  of  propulsion. 

It  is  generally  possible  to  find,  for  each  configuration  of 
the  vehicle  and  for  each  cruising  speed,  an  internal  mode  of  pro¬ 
pulsion  that  produces  smaller  velocity  fluctuations  in  the  adjacent 
fluid  columns  -  and,  accordingly,  require.®  less  power  -  than  any  ex¬ 
ternal  mode.  The  potential  superiority  of  internal  over  external 
modes  is  small  at  low  speeds  but  increases  with  increasing  cruising 
speed  and  becomes  impressive  at  high  subsonic  and  at  supersonic 
speeds. 

Some  merit,  from  the  standpoint  of  power  required,  may  be 
found  under  special  conditions  in  mixed  schemes,  involving  both 
internal  and  external  propulsion. 

During  acceleration  with  internal  propulsion,  the  air  just 
ahead  of  the  vehicle  is  induced  to  flow,  relative  to  the  tube^  Jp, 
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the  direction  of  the  acceleration  of  the  vehicle  or  in  the  oppc 
direction,  depending  on  the  speed  of  the  vehicle  and  on  the  exc 
of  thrust  over  drag.  An  internally  propelled  vehicle,  acceler*. 
from  rest  to  cruising  speed,  generetes  forward- facing^  axpansioi 
coiapressior  waves,  in  this  order  of  succession.  The  net  effect 
these  waves  ahead  of  the  vehicle,  in  the  final  steady  state,  is 
most  cases  a  forward-facing  flow  if  the  cruising  speed  is  subt** 
In  general,  it  is  only  at  supersonic  speeds  that  an  i  ternal  ti 
generator  can  eliminate  all  disturbances  ahead  of  the  vehieJe  ' 
out  producing  an  excess  of  thrust  over  drag. 

Nevertheless,  even  at  subsonic  speeds,  the  flow  disturbanc 
ahead  of  the  vehicle  can  be  reduced,  by  appropriate  means  of  ii 
nal  propulsion,  to  a  much  lower  amplitude  than  would  be  produce 
external  propulsion.  As  a  consequence,  the  increment  of  stream 
force  (relative  to  the  vehicle)  that  is  imparted  by  wall  fricti 
to  the  disturbed  flow  ahead  of  the  vehicle  is  normally  smaller 
internal  than  with  external  propulsion. 

The  flow  behind  the  vehicle  undergoes  for  some  distance  ce 
rapid  transformations.,  the  nature  of  which  has  been  examiried  in 
tail  by  Hagerup  in  a  recent  study  (Ref.  4).  The  subsequent  beh 
of  the  depends  >  i  the  speed  of  the  vehicle  and  on  the  mod 

thrust  generation.  ...  '  ough  there  is  much  in  the  behavior  of  tl 
wake  that  Is  still  incompletely  understood,  it  has  been  ascerta: 
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that  there  is  a  range  of  operating  regimes  within  which  the  bound¬ 
ary  conditions  demand  that  the  wake  flow  be  substantially  at  rest 
relative  to  the  tube  everywhere  except  in  close  proximity  to  the 
vehicle.  Outside  of  this  range  the  wake  is  "stable",  in  the  sense 
that  even  relatively  large  local  perturbations  within  it  may  be 
compatible  in  steady  state  with  the  boundary  conditions.  Of  course, 
if  in  the  vehicle-fixed  frame  of  reference  the  wake  flow  is  sonic 
or  supersonic  at  any  station,  the  downstream  boundary  conditions 
can  have  no  effect  on  the  flow  upstream  of  that  station. 

When  flow  disturbances  are  sustained  in  steady  state  over 
great  distances  ahead  of  the  vehicle,  the  stagnation  temperature 
of  the  air,  relative  to  the  vehicle,  is  lower  than  it  would  be  in 
flight  at  the  same  Mach  number  through  the  surrounding  atmosphere. 

At  high  travel  speeds,  at  least  a  portion  of  the  power  input 
for  internal  propulsion  must  be  supplied  ahead  of  the  throat  of 
the  air  transfer  passage.  The  purpose  of  this  arrangement  is  to 
reduce  the  amplitude  of  the  flow  disturbances  ahead  of  the  vehicle 
at  high  subsonic  speeds  and  to  prevent  the  propagation  of  such 
disturbances  altogether  at  suporsonic  speeds. 

At  very  high  speeds,  thermal  modes  of  internal  propulsion,  in 
which  the  input  of  energy  (except  for  the  relatively  small  fraction 
which  must  be  supplied  ahead  of  the  throat)  is  in  the  form  of  heat 


addition  to  the  flow,  promise  attractive  advantages  of  econonqr 
and  mechanical  simplicity.  These  modes  also  lend  themselves  to 
interesting  possibilities  in  the  utilization  of  heat  from  nucU 
reactors . 

Special  attention  should  be  given  to  thermal  modes  in  whii 
the  heat  sources  are  fixed  to  the  tube  and  successiveljr  activ'*’ ^ 
in  such  synchromism  with  the  passing  vehicle  that  they  produce 
substantially  steady  regime  of  heat  addition  in  the  vehicle-£i> 


frame  of  reference. 
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1.2606-03 

2.132 

1.2566-03 

2.13 

1.2526-03 

2-136 

1. 2486-03 

2.138 

1.2446-03 

2-140 

1.2406-03 

2-142 

1.2366-03 

2.144 

1.2326-03 

2.146 

1. 2286-03 

2.  148 

1.2256-03 

2-150 

1.2216-03 

2.152 

1.2176-03 

2.154 

1.2136-03 

2.156 

1.2096-03 

2.158 

1.2056-03 

2.160 

1.2026-03 

2.162 

1.1986-03 

2.164 

1.1946-03 

2.166 

1.1906-03 

2.168 

1,1876-03 

2.170 

1.1836-03 

2.172 

1.1796-03 

2.174 

1.1756-03 

2.176 

1.1726-03 

2.178 

1. I  686-03 

2.180 

1.1646-03 

2.182 

1.1616-03 

2.184 

1.1576-03 

2.186 

1.1536-03 

2.188 

1.1506-03 

2.190 

1.1466-03 

2.192 

1.1436-03 

2.194 

1. 1396-03 

2.196 

1.1356-03 

2.196 

1.1326-03 

2.200 

1. 1286-03 

M 

Z 

2.202 

1.1256-03 

2.204 

2.206 

1.1186-03 

2.208 

1.1146-03 

2.210 

1.1116-03 

2.212 

1.1076-03 

2.214 

1.1046-03 

2.216 

1. 1006-03 

2.218 

1.0976-03 

2.220 

1.0935-03 

2.222 

1.0906-03 

2.224 

1.0876-03 

2.226 

1.0836-03 

2.228 

1.0806-03 

2.230 

1.0766-03 

2.232 

1.0736-03 

2.234 

1.0706-03 

2.236 

1.0666-03 

2.238 

1.0636-03 

2.240 

1.0606-03 

2.242 

1.0566-03 

2.244 

1.0536-03 

2.246 

1,0506-03 

2.248 

1.0476-03 

2.250 

1.0436-03 

2.252 

1.0406-03 

2.254 

1.0376-03 

2.256 

1.0346-03 

2.258 

1.0306-03 

2.260 

1.0276-03 

2.262 

1.0246-03 

2.264 

1. 0216-03 

2.266 

1.0186-03 

2.263 

1.0146-03 

2.270 

1.0116-03 

2.272 

1.0086-03 

2.274 

1.0056-03 

2.276 

1.0026-03 

2.278 

9.9876-04 

2.280 

9.9566-04 

2.362 

2.284  9.694£-“04  2.364 

2.286  9.863E-04  2.366 

2.288  9.032e~04  2.368 

2.290  9.802t-04  2.370 

2.292  9.77ir— 04  2.372 

2.294  9.7416-04  2.374 

2.296  9.niE-04  2.376 

2.298  9.6816-04  2.378 

2.300  9.6516-04  2.380 

2.302  9.6216-04  2.382 

2,304  9.5916-04  2.384 

2.306  9.5616-04  2.386 

2.308  9.5326-04  2.388 

2.310  9.5026-04  2,390 

2.312  9.4736-04  2,392 

2.314  9,4436-04  2.394 

2.316  9.4146-04  2.396 

2.315  9.3856-04  2.398 

2,320  9.3566-04  2.400 

2.322  9.3276-04  2.402 

2.324  9.2986-04  2.404 

2.326  9. 2696-04  2.406 

2.328  9,2416-04  2.408 

-  2.330  9.2126-04  2.410 

2.332  9.1846-04  2.412 

2.334  9. 1556-04  2.414 

2.336  9.1276-04  2.416 

2.338  9.0996-04  2.418 

2.340  9.0716-04  2.420 

2.342  9.0436-04  2.422 

2.344  9.0156-04  2.424 

2.346  8. 9876-04  2.426 

2.348  8.9596-04  2.428 

2.350  8.9326-04  2.430 

2.352  0.9046-04  2.432 

2.354  8.8776-04  2,434 

2.356  8.8506-04  2.436 

2.350  8.8226-04  2.438 

2.360  3.7956-04  2.440 


8.7686-04  2.442  7.7576- 

e. 7416-04  2.t44  7,7346-( 

8. 7146-04  2.446  7.710E-( 

8.6886-04  2.448  7.687E-( 

8.6616-04  2.450  7.6636-^ 

8. 6346-04  2.452  7.640E-{ 

8. 6036-04  2.454  7. 6176-'. 

8.581E-04  2.456  7. 594E-i 

8.  555E-04  2.458  7.571E-'. 

8. 5296-04  2.46C  7.540- 

8.5036-04  2.462  7.5256- 

8.4766-04  2.464  7.5026-1 

8. 4506-04  2.466  7.480E-) 

8.4256-04  2.468  7.457E-C 

8.3996-04  2.470  7.434t 

8.3736-04  2.472  7. 412E-' 


8.3476-04 

2.474 

7. 3896- 

8.3226-04 

2.476 

7.  367E-t' 

8.2966-04 

2.478 

7,3456-^ 

8.2716-04 

2.480 

7.323E-C 

8. 2466-04 

2.482 

7.301E-C 

8.2216-04 

2.484 

7.279E  -1 

8.1956-04 

2.486 

'7.257E-.‘ 

8.1706-04 

2.488 

7.235E-t 

8.1456-04 

2.490 

7.2136-0 

8.1216-04 

2.492 

7.1916-0 

8,0966-04 

2.494 

7. 1696-0 

8.0716-04 

2.496 

7.148E-Q 

8.0466-04 

2.498 

7. 126E-0 

8.0226-04 

2.500 

7. 105E-0 

7.9976-04 

2.502 

7.083E-0 

7.9736-04 

2.504 

7.062E-0 

7.9496-04 

2,506 

7.0416-0 

7.925E-04 

2.508 

7.020E-0 

7.900E-04 

2.510 

6.999E-0 

7.8766-04 

2.512 

6.978E-C 

7.8526-04 

2.514 

6.957E-0 

7.8286-04 

2.  516 

6.9366-0 

7.8056-04 

2.518 

6.9156-0 

7.7816-04 

2.520 

6.894E-0 

y 

Z 

2»5^^ 

6.a73E--04 

^.bZ4 

6..853E-04 

2.526 

6,832E-04 

2.528 

6.812E-04 

2.530 

6.7916-04 

2.532 

6.771E-04 

2.534 

6.7S1E-04 

2.536 

6.731E-04 

2.538 

6.710E-04 

2.540 

6.690E-04 

2,542 

6.670E-04 

2.544 

6.650E-04 

2.546 

6.631E-04 

2.  548 

6.611E-04 

2.550 

6.591E-04 

2.552 

6.571E-04 

2.554 

6.552E-04 

2.556 

6.532E-04 

2.558 

6.5136-04 

2.560 

6.4936-04 

2.562 

6.4746-04 

2.564 

6.4556-04 

2.566 

6.43Dt-a4 

2.568 

6. 4166-04 

2.570 

6.3976-04 

2.572 

6.3786-04 

2.574 

6.3596-04 

2.576 

6.340E-04 

2.578 

6.321E-04 

2.580 

6.303E-04 

2.582 

6.284E-04 

2.584 

6.2656-04 

2.586 

6.2476-04 

2.588 

6.228E-04 

2.590 

6.210E-04 

2.592 

6.igiC-04 

2.594 

6. 173E-04 

2.596 

6.155E-04 

2.598 

6. 136E-04 

2.600 

6.118E-04 

W 

Z 

2.602 

6.100E-04 

2.604 

6.082E-04 

2.606 

6.064E-04 

2.608 

6.046E-04 

2.610 

6.028E~04 

2.612 

6.011E-04 

2.614 

5.993E-04 

2.616 

5.975E-04 

2.618 

5.957E-04 

2.620 

5.940E-04 

2.622 

5.922E-04 

2.624 

5.905E-04 

2.626 

5.888E-04 

2.628 

5.870E-04 

2.630 

5.853E-04 

2.632 

5.836E-04 

2.634 

5.818E-04 

2.636 

5.801E-04 

2.638 

5.784E-04 

2.640 

5.767E-04 

2.642 

5.  750E-04 

2.644 

5.733E-04 

2.645 

5.717E-04 

2.648 

5.700E-04 

2.650 

5.6836-04 

2.652 

5.666E-04 

2.654 

5.650E-04 

2.656 

5.633E-04 

2.658 

5.617E-04 

2.66C 

5.600E-04 

2.662 

5.584E-04 

2.664 

5.568E-04 

2.666 

5.551E-04 

2.668 

5.535E-04 

2.670 

5.519E-04 

2.672 

5.503E-04 

2.674 

5.487E-04 

2.676 

5.471E-04 

2.678 

5.455E-04 

2.680 

5.4J9E-04 

M 

Z 

2.682 

5.423E-04 

2.684 

5.4075-04 

2.686 

5.3916-04 

2.688 

5.376E-04 

2,690 

5.360E-04 

2.692 

5,344E-04 

2.694 

5.329E-04 

2.696 

5.313E-04 

2.698 

5.298E-04 

2.700 

S.283E-04 

2.702 

5.267E-04 

2.704 

5.252E-04 

2.706 

5.237E-04 

2.708 

5.221E-04 

2.710 

5,206E-04 

2.712 

5.191E-04 

2.714 

5.176E-04 

2.716 

5,161E-04 

2.718 

5.146E-04 

2.720 

5.131E-04 

2>722 

5. 116E-04 

2.724 

5.102E-04 

2.726 

5.087E-04 

2.728 

5.072E-04 

2.730 

5.057E-04 

2.732 

5.043E-04 

2.734 

5.028E-04 

2.736 

5.014E-04 

2.738 

4.999E-04 

2.740 

4.985E-04 

2.742 

4.971E-04 

2.744 

4.956E-04 

2.746 

4.942E-04 

2.748 

4.928E-04 

2.750 

4.913E-04 

2,752 

4.899E-04 

2.754 

4.805E-04 

2.756 

4.871E-04 

2.758 

4.857E-04 

2.760 

4.843E-04 

M2  HZ  M 


2.762 

4.829E-04 

2.842 

2.764 

4. 815E-04 

2.844 

2,766 

4.a02E-O4 

2,846 

2,  768 

4. 788E-04 

2.848 

2.770 

4.774E-04 

2.850 

2.772 

4.760E“04 

2.852 

2.774 

4.747E-04 

2.854 

2.776 

4.733E~04 

2.856 

2.778 

4.720E-04 

2.858 

2.780 

4. 706E-04 

2.860 

2.782 

4.693E-04 

2.862 

2.784 

4.679E-04 

2.864 

2.786 

4.666E-04 

2.866 

2.788 

4.652E-04 

2.868 

2.790 

4.6J9E-04 

2.870 

2.792 

4.626E-'04 

2.872 

2.794 

4.613E~04 

2.874 

2,796 

4,600E-04 

2.876 

2.798 

4.586E-04 

2.878 

2.800 

4.573E-04 

2.880 

2.802 

4.560E-04 

2.882 

2.804 

4.547E-04 

2.884 

2.806 

4.534E-04 

2.886 

2.808 

4.521E-04 

•s  tyt%n 

£.  •  obo 

2.810 

4.509E-04 

2.890 

2.812 

4.496E-04 

2.892 

2.814 

4,483E-04 

2.894 

2.316 

4.4706-04 

2.896 

2.818 

4.458E-04 

2.898 

2.820 

4.445E-04 

2.900 

1.8  22 

4.432E-04 

2.902 

2.824 

4.420E-04 

2.904 

2.826 

4.407E-04 

2.906 

2.828 

4o395E-04 

2.908 

2.830 

4.382E-04 

2.910 

2.832 

4.370E-04 

2.912 

2.834 

4.357E-04 

2.914 

2.836 

4.345E-04 

2.916 

2.338 

4.333E-04 

2.918 

2.840 

4.320E-04 

2.920 

4. 308E-04 

2.922 

3.850E-0 

4.296E-04 

2.924 

3.839E-L 

4.284E-04 

2.926 

3.829E-0 

4.272E-04 

2.928 

3.8I.8E-0 

4,260E-04 

2.930 

3.807E-0 

4,248C-04 

2.932 

3.797E-*v 

4.236E-04 

2.934 

3.786E"a 

4.224E-04 

2.936 

3. 776C-C 

4.212E-04 

2.938 

3.765E-0 

4.200E-04 

2.940 

3.755E-« 

4. 188E-04 

2.942 

3.744E-.. 

4.  176E-04 

2.944 

3.734E-0 

4.165E-04 

2.946 

3.724E-0 

4.153E-04 

2.948 

3.713E-C 

4. 141E-04 

2.950 

3.703''-0- 

4.130E-04 

2.952 

3.6935;~f 

4.  118E-04 

2.954 

3.683E-C 

4.106E-04 

2.956 

3.672E-C 

4.095E-04 

2.958 

3.662E-C 

4.083E-04 

2.960 

3.652E-0 

4.072E-04 

2.962 

3.642E-0 

4.060E-04 

2.964 

3.632E'’^ 

4.049E-04 

2.966 

3.622E-0 

4 ,038E-04 

2.968 

3.612E-0 

4.026E-04 

2.970 

3,602E-0 

4.015E-04 

2.972 

3.592P-0- 

4.004E-04 

2.974 

3.582E-0 

3.993E-04 

2.976 

3k572E-0' 

3.981E-04 

2.978 

3.563E-0‘ 

3.970E-04 

2.980 

3.553E-0' 

3.959E-04 

2.982 

3.543E-0* 

3.948E-04 

2.964 

3.533E-0 

3.937E-04 

2.986 

3.523E-0- 

3,926E-04 

2,988 

3,514E-0« 

3.915E-04 

2.990 

3.504E-0' 

3.904E-04 

2.992 

3.495E-0 

3.893E-04 

2.994 

3.485E-0 

3.882E-04 

2.996 

3.475E-0 

3.872E-04 

2.998 

3.466E-0< 

3.861E-04 

3.000 

3.<»56E-0' 
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SECTION  G 


) 


ANALYSIS  OF  THE  WAKE  OF  A 
SELF-PROPELLED  VEHICLE  IN  A  DUCT 
WITH  CONSIDERATION  OF  HEAT  TRANSFER  AND  FRICTION 


by 


Duane  E«  Cromack 


and 

)  ^  * 

I  Leo  F.  Gijsels 


This,  report  coocei^na  tho  flow  in  the  of  ^  seif-, 

propelled  vehicle  traveling  In  an  Infinite  tube#  The  analysle 
Includes  the  effects  ^f  conipresslblUty,  friction  and  heat 
transfer  and  makes  use  of  the  Reynolds  analogy  to  relate  the 


latter  two  effects. 


SYMBOLS 


a 

Cr 

D  . 
F 

g 

h 

K 

in 

M 

P 

pO 

% 

r 

R 


speed  of  sound  (fps) 

specific  heat  at  constant  pressure  (ft-lb/slug) 
non-dimensional  friction  factor 

M  (i  +  v-1 

Stream  force  relative  to  the  vehicle  (lb) 

2 

gravitational  constant  (ft/sec  ) 

convective  heat  transfer  coefficient  (Btu/ft  sec  °R) 
heat  transfer  factor  in  eq.  3  (ft-lb/slug  °R) 
mass  flow  rate  relative  to  the  vehicle  (slugs/sec) 
Mach  number  of  the  flow  relative  to  the  vehicle 

s  const. 


^VgRTo  J 


VgRT; 

static  pressure  (psf) 
stagnation  pressure  (  *) 

heat  input  in  combustion  process  per  unit  mass  (ft-lb/slug) 
duct  radius  (ft) 
gas  constant  (ft-lb/slug  ^R) 

T  static  temperature  (^R) 

T®  stagnation  temperature  relative  to  the  vehicle  (°R) 

stagnation  temperature  relative  to  the  wall  (°R) 

< 

a  velocity  of  the  flow  relative  to  the  vehicle  (fps) 

u^  velocity  of  the  wall  relative  to  the  vehicle  (fps) 

V  velocity  of  the  vehicle  relative  to  the  wall  (fps) 

3^  duct  axial  coordinate  (ft) 

C-11 
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o<-  duct  cross-sectional  area  (ft^) 

^  non-dimenaional  velocity  ratio,  u/u^ 

V  specific  heat  ratio 

non-dimensional  axial  coordinate,  x/r 

Subscripts  (see  Fig.  G-1) 

0  conditions  outside  the  tube,  assumed  uniform 

1  nearest  station  ahead  of  the  vehicle  where  the  flow  can  be 
considered  one-dimensional 

2  station  at  tail  end  of  vehicle 

3  terminal  station  of  energy  addition  in  internal  propulsion 
upstream  end  of  wake 
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INTRODUCTION 

The  general  concept  of  high  speed  vehicles  propelling  them¬ 
selves  through  ducts  was  presented  by  Foa  in  reference  1.  In  the 
calculation  of  the  power  required  for  such  systems,  the  flow  in¬ 
duced  by  the  vehicle  must  be  determined. 

This  problem  has  been  treated  in  different  ways  by  various 
authors  (Refs.  2  through  6).  Tollmien  (Ref.  2)  was  the  first  to 
study  the  flow  produced  by  a  train  in  a  tunnel  and  the  effect  of 
this  flow  on  the  power  required.  Reference  3  considers  the  flow 
induced  by  a  self-propelled  subsonic  vehicle  through  a  quasi  static 
analysis  in  which  viscous  stresses  are  neglected  and  the  thrust 
generator  is  repres^inted  as  an  actuator  disc.  Reference  4  also 
confjiders  an  incomiaressible  and  inviscid  fluid  in  a  simplified 
treatment  of  the  same  problem.  Reference  5  includes  viscous  effects 
and  follows  a  boundary- layer  approach.  Reference  6  develops  the 
general  equations  for  the  one -dimensional  ans lysis  of  the  flow  in 
the  frame  of  reference  of  a  vehicle  moving  in  a  tube  of  infinite 
length,  with  consideration  of  friction  and  heat  transfer.  Closed- 
form  solutions  are  obtained  for  regions  ahead  of  and  at  some  dis¬ 
tance  behind  the  vehicle,  where  heat  exchanges  are  assumed  to  take 
place  according  to  such  a  schedule  that  the  stagnation  enthalpy 
relative  to  the  walls  remains  constant. 


In  the  present  analysis,  the  equations  of  Ref.  6  are  applied 
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to  the  region  behind  the  vehicle  without  the  sinq)lifylng  a8Suaf>Cion8 
that  are  required  to  yield  closed-form  solutions. 


ANALYSIS 

The  present  analysis  considers  a  tube  of  infinite  length  in 
which  a  vehicle  is  moving  at  constant  speed.  In  a  vehicle-fixed 
frame  of  reference  (see  Fig.  G-1)  the  flow  is  steady.  The  walls 
are  moving  at  a  velocity  equal  and  opposite  to  the  velocity  of  the 
vehicle.  The  energy  and  momentum  equations  for  this  situation  art 


(Ref.  6); 


rh^CftaLT*"  ^  Uar  cLF  'f' 


where  ot  T 


dr  -h  udi 


dr  *  -^o(^{u-Uur)\U^U^^\dK 

The  wall  temperature  Is  assumed  here  to  remain  constant  and  equal 
to  the  external  ambient  temperature  T^.  Thus,  the  Instantaneous 
rate  of  heat  transfer  can  be  expressed  as 

-  -K(a)(:r°‘-7:)dx  (3) 

wh^re  T^*  is  the  stagnation  tenperature  of  the  flow  in  the  vicinity 
of  and  relative  to  the  wall.  Thus 


J-  ^  \^-Uur) 


i 
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K(u)  is  a  function  of  the  convective  film  coefficient  and  through 
the  Reynolds  analogy  becomes  a  function  of  the  friction  coeffici¬ 
ent  and  of  the  absolute  velocity  of  the  flow  relative  to  the  wall, 
|u  -  u^j  .  From  the  Reynolds  analogy,  the  Stanton  Number  for  the 
case  of  moving  walls  can  be  written  as 

3.  =  h _ ^ 

9/^Cplu-Uurl  2 

Thus  the  film  coefficient  becomes 

and  K(u)  becomes 

K(uJ  =  2rrr/7  -  Trt'CfCp^^\U-Uar\  . 

Substitution  of  equations  2,  3,  and  4  into  1  gives 

mgCpdT  i-  mudu.  - 

(5) 

-irfCrCpqp\u-ujk*f^^:}^-%)dx  . 

The  mass  flow  rate  m,  may  be  written  as 


and  differentiating  gives 


rh^RdT  a  dF  -h  Fdu  -  zm udu  .  (6) 


Eqs.  5  and  6  give 


^_i^^^uclF  +  F'ctu-2/Kudt^  +  udu  — 
--^‘^(u-t^u-u^dx  -  ffCiColu-ujfFu  _u‘  (u-u^f  _T-) 

^  r  u  \rt>9R  W  2  9  4,“  y 

Substituting  for  dF  from  eq.  2  and  rearranging,  one  obtains 


Differentiation  with  respect  to  x  and  the  elimination  of  F/m 
between  the  resulting  equation  and  eq.  7  leads  to  the  following 
second  order  non-linear  differential  Equation  in  non-dimensional 
form  with  *• 
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\Cr  4±i)' 


The  values  for  /j  and 


d/3 


at  O  i.e.,  immediately 


behind  the  vehicle,  are  obtained  in  the  following  manner:  Re- 

rJ  /S 

arranging  eq.  7  in  terms  of  -~=J. — -  gives 

(9) 


dA  _ \cr\A-i\  ll/v  (ibOV'i)  -¥  Jj\ I 

[jflUur  -*  L  J 

where  all  terras  can  be  eva'^uated  at  7^s?o  for  a  given  /<5(0). 
The  initial  value  of ^  must  be  expressed  as  a  relation  between 
appropriate  parameters  of  the  propulsive  device, 
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Equation  8,  with  the  appropriate  initial  values  of 


and 


d\ 


may  be  used  to  describe  the  flow  in  the  wake  of 


a  subsonic  or  a  supersonic  vehicle  traveling  in  a  tube. 
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coefficients  of  velocity  transform 

specific  heat  ratio 

initial  mass  flow  disturbance 

pressure  increment  relative  to  initial  value,  A||» 

initial  velocity  disturbance 

initial  enthalpy  disturbance 

arbitrary  function 

enthalpy  perturbation 

Modified  Bessel  functions  of  the  first  kind 
indices 

Modified  Bessel  functions  of  the  second  kind 
constant 

coefficients  of  enthalpy  transform 
MSch  number  upstream 
viscosity  perturbation 
pressure  perturbation 


h 

<• 

heat  transfer  function,  ^ 

K 

A.. 

Rx. 

f 


tube  radius 
radial  coordinate 


Reynolds  number,  s  g ^ 

r* 

density  perturbation 
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^  Prandtl  number 
W.  axial  velocity  perturbation 
U*  radial  velocity  perturbation 
jL  axial  coordinate,  positive  dowstream 
Laplace  transform  variable 

‘Ti’&N  shear  function,  (A  X  —  I 

''  IwJUA 

^  Laplace  transform  of  (.  ) 

(  ^  dlraens-fnnal  ph3'sical  quantity  conqolete 
(  dimensional  upstream  reference  value 
(  value  immediately  behind  disturbance 
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Summary 


Linearized  two-dimensional  analysis  indicating  some  general 
features  of  flow  induced  by  disturbance  propagating  at  constant 
velocity  along  slender  tube.  Formulation  applies  to  small  dis¬ 
turbances,  includes  effects  of  viscosity  and  heat  conduction  for 
gases  of  arbitrary  Prantl  number,  neglects  viscous  dissipation. 
Typical  subsonic  and  supersonic  solutions  have  been  calculated 
for  case  of  Prandtl  number  one. 
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A  new  transportation  system  in  which  a  vehicle  propels 
itself  at  high  velocity  through  a  long  tube  has  been  proposed 
(Ref.  1).  The  investigation  of  the  turbulent  flow  field  in  the 
tube  downstream  of  the  vehicle  forms  an  important  part  of  the 
complete  system  analysis.  The  purpose  of  the  present  discuesion 
is  to  explore  in  simpler  terms  some  aspects  of  the  induced  flow 
field  associated  with  a  local  disturbance  on  a  more  modest  scale, 
traveling  at  constant  velocity  along  an  essentially  infinite  tube. 

To  an  observer  stationary  in  the  tube,  the  passage  of  the 
disturbance  is  felt  as  a  sudden  change  in  velocity,  pressure  and 
temperature,  followed  by  a  gradual  return  to  the  undisturbed  state. 
To  an  observer  traveling  with  the  origin  of  the  disturbance,  it 
might  reasonably  be  expected  that  the  flow  would  appear  steady,  and 
the  conservation  equations  would  be  simpler  in  this  case;  however, 
the  assumption  of  steady  flovj  is  of  course  admissible  only  if  the 
induced  flow  field  predicted  on  this  basis  decays  asymptotically 
with  distance  downstream,  so  that  the  etid  conditions  are  automati¬ 
cally  satisfied. 

The  two-dimensional  analysis  reported  in  this  preliminary  note 
is  confined  to  the  laminar  version  of  the  steady  flow  problem.  The 
problem  situation  referred  to  the  moving  coordinate  system  is 
formalized  in  the  following  terms:  a  frontal  discontinuity  at  some 
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station  X~0  in  a  tube  of  circular  cross  section  imparts  to 
the  flew  immediately  behind  it  a  velocity  u!*s  parallel 

to  the  axis  and  an  enthalpy  ,  both  perturbations 

here  chosen  as  initially  constant  across  the  tube.  Other  dis¬ 
tributions  would  however  be  equally  admissible  within  the  two- 
dimensional  formulation.  The  walls  as  well  as  the  fluid  upstre?%' 
are  taken  to  be  at  conditions  and  though  the  up¬ 

stream  state  enters  the  mathematical  problem  only  as  a  reference 
state.  The  objective  of  the  analysis  is  to  describe  the  develop¬ 
ment  of  the  initial  block  perturbation  profiles  and  Iv,  B 

in  the  flow  field  downstream. 


The  problem  may  for  present  purposes  be  regarded  as  in  the 
nature  of  a  boundary  layer  situation,  with  x-derivatives  and 
r-conponents  of  velccity  generally  negligible  compared  to  r- 
derlvatives  and  x-components ,  and  x-derivatives  of  enthalpy 
generally  negligible  coiiq}ared  to  r-derivatlves .  Assuming  the 
flow  also  to  be  axi85rmmetric  and  without  swirl,  we  have  the  simpli 
fled  governing  system  of  equations  for  the  steady  flow  of  a  per fee 


gas  with  constant  Prandtl  number  in  the  form 
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The  boundary  and  initial  conditions  are 

A.'*  r  Tk  ^  ~ 

s  0  vr^rO 


Because  of  the  finite  velocity  at  the  wall,  the  local  velocity 
in  the  flow  field  need  nowhere  depart  significantly  from  some 
reasonably  defined  mean,  and  the  Oseen  approximation  suggests 
itself  as  a  promising  approach  to  the  problem.  The  equations 
are  accordingly  linearized  by  replacing  the  transport  derivative 
operator  as  follows 


where  is  an  unspecified  consent,  which  in  the  small  perturba¬ 
tion  analysis  to  follow  may  be  taken  as  1  •  Introducing  non- 

dimensional  variables  through  the  expressions 

Ka.  X**?:  A.*'* 


one  obtains  under  the  restriction  to  small  perturbations  the 
linearized  formulation  for  the  compressible  flow  development  in 
terms  of  the  normalized  system 


L  (s  iiA  -  •=  — 

A  iN  V  yn'-  JU 


where  it  should  be  noticed  in  particular  that  the  viscous  dis¬ 
sipation  term  in  the  energy  equation  has  been  neglected,  as  a 
higher  order  effect.  The  system  is  subject  to  the  ^nditions 

K  9  \  W  =  •  O 

j4*o 


The  nature  of  the  problem  suggests  now  the  use  of  the  Laplace 
transformation 

H  ^  WC  N.)  A ^ 
o 

Transforming  first  the  momentum  equation,  one  obtains  an  ordine» 
differential  equation 


^  WS.  r 


-  5 


which  has  the  solution  (Ref.  2) 


^  (?i^)  *  c,  -V  C.^ 

The  axis  N»0  is  within  the  flow  field,  so  must  be 

bounded;  hence  C^»0  .  The  velocity  boundary  condition  at  the 
wall  ^s|  gives 


But  the  p/jrtlcular  solution  is  by  inspection 

mmA  ^ 

*  **  . <4|l  (IImIiI  I^i  iw  ,&>k 

'  ^  y  H'*  ^ 
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with  the  initial  pressure  given  by  the  equation  of  state  as 
-  A -'S g.  where  A  is  a  term  allowing  for  the  possibility 
of  mass  addition  as  part  of  the  initial  disturbance  mechanism. 
Solving  for  the  tranj formed  pressure  one  obtains 

A  third  constraint  is  provided  by  the  constant  mass  flow  require¬ 
ment  for  steady  flow: 

R 

0 

Normalizing  this  and  taking  the  Laplace  transform,  one  obtains 
with  the  equation  of  state  the  mass  flow  constraint  in  the 
linearized  form 

o 

From  the  energy  equation,  next,  it  follows  by  transformation 

and  the  admissible  part  of  the  solution  is  again  of  the  form 


with 
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from  the  wall  temperature  boundary  condition.  The  particular 
solution,  using  the  pressure  as  already  obtained,  is 


Substitution  of  the  expressions  for  u.  ,  ,  and  ja  into  the 
mass  flow  constraint  provides  an  implicit  relation  for  c.^  an, 
upon  integration  leads  to  the  result  eventually,  neglecting 


perturbation  products: 

c  S  -  -1—  Cl-f  -  [£■»  I.Cvm) 

^  IJ'S),, 


The  analysis  is  capable  of  greater  accuracy,  but  for  the  presr 
exploratory  purpose  the  first  approximation  is  adequate. 


Rather  than  to  attempt  at  this  point  the  description  of  t 
entire  velocity  and  temperature  fields,  we  confine  the  follow! 
discussion  to  the  more  essential  parameters  of  wall  shear,  hea 
transfer,  and  pressure  ,as  functions  of  x.  The  normalized  shea; 
transform  is  given  by  the  expression 

'Kt\ 

The  theory  as  formulated  applies  to  a  gas  of  arbitrary  Prandtl 
number,  but  in  the  following  we  consider  the  special  case  of 
I  for  algebraic  convenience.  The  shear  transform  follow*. 
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then  from  the  general  expression  for  C,  in  the  relatively  sinq)le 


form: 


S-s. 


'fX  \  -  n 


I  sfA  un'  x«\\r?)J 

The  transform  of  normalized  heat  transfer  to  the  wall  is  given  by 


a  K 


•A  *  ( 


-  *■  r,  (viStf) 


and  it  follows  similarly  that 


vJa  \  -in'* 


ntJi-  ^  1 


The  transform  of  the  pressure  change  referred  to  the  initial 
condition  is  given  as 


v  — •  ^  A*-T'"£ 

■IP  =  r  -  ‘^-r— 


vn"[£  ♦  Xoi^^eJ 


and  takes  the  form 


'X^ 

VT^ 


S-£. 

l  -  fn»* 


T.OK)  ^ 

1.  ^  r,ty7J)i 

The  remaining  problem  is  to  establish  the  Inverse  transforms  of 
the  functions  ,  and  Af  0^)  . 


The  inversion  is  most  readily  accomplished  by  power  series 
expansion  and  by  residue  calculus  for  the  two  ranges  corresponding 
to  small  X  and  large  x  respectively.  For  the  initial  part  of  the 
flow,  the  asymptotic  representations  for  Bessel  functions  of  large 


arguinent  lead  to  the  expansion 


=  s  {4'  r*"  swi'  ■*■  C*  ^  * 


a.«c-i 

/ivfjj 


\  -M*- 


The  inverse  transform  representing  the  normalized  local  shear 
for  small  x  follows  from  this  in  the  form 

^  lx  +  +  e  (xv^i3 

For  the  heat  transfer  one  obtains  in  the  same  manner  the  result 

\cz  +  •<.  ^  0  (kw^ 

and  for  the  pressure  change  similarly: 

Ar(.«)=  -(?-«-) 


The  a8yin)totic  expansions  of  the  same  flow  properties  for  large 
values  of  x  are  obtainable,  with  a  little  caution,  from  the 
simple  expansion  foxrmula  (Ref.  3) 

‘•-1  w]  •-*:■{  ^5} 

when  the  singularities  of  4.(4)  are  simple  poles  at  . 

For  the  shear  one  obtains  for  large  x  the  result 
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where  the  location  of  the  singularities  follow  from  the  trans¬ 
cendental  equations 


T,  =0  ^ 

^  1  f  ■ 

1 0  ' 

The  heat  transfer  is  given  for  large  x  similarly  by 

^  vsiV  / 

-  E|cs-c)fii^]' _ ikiinl _ _]/i 

while  the  pressure  change  for  large  x  is  •  ^ 


in  terms  of  the  roots  of  the  second  transcendental  equation  above. 


Results  of  sample  calculations  of  shear,  heat  transfer,  and 
pressure  variation  along  the  tube  for  C  •  I  and  for  character¬ 
istic  velocities  upstream  corresponding  to  M  =  0,  .5,  1.5,  and  2.5 
are  presented  at  the  end  of  this  note.  The  first  figure  shows  the 
steady  flow  response  of  the  wall  shear  to  unit  initial  block  pro¬ 
file  perturbations  in  velocity  and  in  enthalpy  ©  ,  the  second 
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figure  shows  the  corresponding  heat  transfer  response,  and  the 
third  figure  shows  the  pressure  response  to  a  velocity /enthalpy 
Input  sufficient  to  give  a  unit  Initial  perturbation  In  pressure. 
The  results  are  Interesting  in  their  own  right,  being  directly 
pertinent  to  the  conq)re8slble  version  of  the  classical  Inlet 
problem  with  or  without  moving  walls.  This  and  related  applica¬ 
tions  will  be  discussed  In  detail  elsewhere. 

In  terms  of  the  problem  of  present  concern,  in  the  supersoffi 
case  the  linearised  solution  as  based  on  the  steady  flow  assumptk 
displays  the  proper  asynqjtotic  decay  and  is  thus  acceptable  for  tt 
vehicle-in-the-tube  problem  situation:  the  failure  of  the  end  co.? 
dltlon  to  cojiinunlcate  upstream  In  this  case  could  Indeed  have  beet 
anticipated  by  considering  the  signal  propagation  within  the  flow 

c 

field  under  supersonic  conditions.  The  results  shown  In  the  flnaJ 

figure,  representing  for  illustrative  purposes  the  response  of  the 

flow  at  M  s  1.5  to  Initial  perturbations  of  107.  velocity  decrement 

and  25%  enthalpy  rise,  teveal  some  Interesting  circumstances  of 

the  laminar  problem.  As  reflected  in  the  wall  shear  variation,  tt 

downstream  flow  field  turns  out  for  these  conditions  to  be  charSet 

ized  by  a  temporary  reversal  in  flow  direction  rather  than  by  the 

< 

Initially  expected  gradual  deoay.  A  significant  feature  Is  the 
finite  heat  transfer  at  the  station  whete  the  wall  shear  passes 
through  zeroj  obviously  the  use  Reynolds  analogy  on  the  problefi 

must  be  handled  with  per titular  tarei 
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In  the  subsonic  case,  on  the  other  hand,  the  steady  flow 
assumption  leads  to  a  response  which  according  to  the  dissipation¬ 
less  linearized  theory  in  general  diverges  exponentially  for  large 
X.  The  solution  and  one  or  more  of  its  underlying  assumptions  are 
therefore  fundamentally  inapplicable  to  the  subsonic  vehicle-in- 
the-tube  problem  situation,  except  when  (5  -  0-  0  .  According  to 
the  linearized  theory,  however,  in  this  particular  case  the  pressure 
perturbation  is  everywhere  zero,  the  divergent  terms  in  shear  and 
heat  transfer  cancel  out,  and  the  subsonic  solution  then  satisfies 
the  end  conditions  after  all.  This  special  situation  recalls  super¬ 
ficially  the  problem  of  the  propagation  of  a  laminar  flame  at  essen¬ 
tially  constant  pressure,  but  it  remains  to  be  seen  if  this  precarious 
balance  of  oppositely  divergent  terms  is  the  only  steady  subsonic 
solution  of  Interest, 

A  useful  next  application  of  the  linearized  theory  would  be 
to  the  initial,  unsteady  phase  of  the  problem  when  the  disturbance 
is  nearer  to  the  end  of  the  tube,  in  order  to  check  the  extent  to 
which  the  subsonic  asymptotic  behavior  for  large  time  would  re¬ 
flect  the  constant  pressure  solution  indicated  above.  The  role 
of  viscous  dissipation  in  the  problem  should  also  be  reconsidered, 
but  that  is  outside  the  scope  of  the  linearized  formulation.  An¬ 
other  problem  outstanding  is  the  much  more  difficult  question  of 
turbulent  effects.  The  quantitative  results  derived  above  apply 
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of  course  specifically  to  the  laminar  problem;  however,  the 
general  conjectures  suggested  are  basic  to  the  problem  situa¬ 
tion  and  would  appear  to  apply  to  the  turbulent  case  as  well. 
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SECTION  I 


ON  TWO  PROBLEMS  OF  TRANSONIC  NONSTEADY 
MOTION  OF  A  VEHICLE  IN  A  TUBE 


John  W.  Burr 


SYMBOLS 


^  -  speed  of  sound  (ft/sec) 

A  -  speed  of  sound  nondimenionalized  (  ) 

b  -  acceleration  of  vehicle  (ft/sec^) 
b  -  acceleration  of  the  air  (ft /sec  ) 

S  -  percent  mass  flow  blead  off 
D  -  Mach  number  function  (see  page  1-8) 
h  -  stream  force 

G  -  Mach  number  function  (see  page  1-8) 

0  -  acceleration  of  gravity 

^  -  Mach  number 

0 

-  mass  flow  rate 

-  Mach  number  function  (see  page  1-8) 

?  -  Riemann  parameter  on  right  running  characteristic 

P  -  pressure  (lb/ft2) 

Q  -  Riemann  parameter  on  left  running  characteristic 
R  -  gas  constant 
T*  -  temperature  (°R) 
r  -  time 

'U  -  velocity  (ft/sec) 

T7  -  velocity  nondimens ionalized  (  ^ 

X  -  distance  along  tube 
(X  -  cross  sectional  area 
y  -  specific  heat  ratio 
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Subscripts 


L  -  initial  values 
0  -  reference  state 

R  -  pertaining  to  the  flow  on  the  right  hand  side  of 
the  vehicle  relative  to  the  tube 

S  -  shock 

V  -  vehicle 

VR  -  pertaining  to  the  flow  on  the  right  hand  side  of 
the  vehicle  relative  to  the  vehicle 

3  -  region  in  tirae-distance  plane  (See  fig.  I-l) 

V  -  region  in  time-distance  plane  (See  fig.  I-l) 


srscripts 


0  -  stagnation  conditions 
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INTRODUCTION 


The  problem  of  accoi -“rating  through  the  transonic  region  jto 
reach  supersonic  speeds  has  always  been  a  difficult  problem  and 
is  even  more  difficult  for  the  case  of  a  vehicle  in  a  tube.  When 
a  vehicle  travels  at  supersonic  speeds,  either  through  a  tube  or 
in  open  air,  there  usually  exists  a  shock  wave  at  the  nose  of  the 
vehicle.  It  may  be  attached  to  the  vehicle  or  be  only  a  small 
fraction  of  a  body  length  ahead  of  the  vehicle,  but  there  is 
always  a  shock  wave  present.  When  a  self-propelled  vehicle,  which 
uses  the  air  in  the  tube  as  the  propulsive  medium,  accelerates 
from  rest  in  a  tube  it  usually  sends  out  an  expansion  wave  initi¬ 
ally.  Since  a  shock  wave  will  exist  ahead  of  the  vehicle  in  super 
sonic  flow,  there  must  be  some  point  during  the  acceleration  ’at 
which  the  wave  being  sent  upstream  changes  from  an  expansion  to  a 
compression  type  wave. 

This  acceleration  of  a  vehicle  in  a  duct  through  transonic 
speeds  was  first  investigated  in  1957  and  reported  in  RPI  report 
TR  AE  6101  (Ref.  1).  This  investigation  showed  that,  although  it 
may  he  possible  to  reach  supersonic  speeds,  many  unusual  and  inter 
esting  flow  situations  occur  during  the  transonic  acceleration 
period. 

The  problem  of  ob^  .initig  supersonic  speeds  within  a  tube  is 
very  similar  to  the  problem  of  starting  a  supersonic  inlet.  The 
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presence  of  the  vehicle  in  the  tube  form  a  "throat*'^  and  before  a 
supersonic  type  flow  is  cbtaino'i  the  normal  shock  which  forms  in 
the  tube  ahead  of  the  vehicle  must  be  ’‘swallov’'ed*'.  At  some  sub¬ 
sonic  Mach  number  the  flow  around  the  vehicle  will  reach  sonic 
velocity  relative  to  the  vehicle  and  will  become  "choked*^  As 
the  vehicle  accelerates  beyond  this  point,  the  vehicle  will  be 
forced  to  "push  along"  seme  of  the  air  cbead  of  it  and  the  distur¬ 
bances  which  are  sent  upstream  will  no  longer  be  expansion  waves 
as  in  the  initial  acceleration  but  rather  compression  waves.  ThuK' 
compression  waves  will  coalesce  and  form  a  normal  shock  some  dia- 
tance  ahead  of  the  vehicle  and  the  flow  around  the  vehicle  will 
remain  choked  until  the  vehicle  overtakes  the  shock.  This  can  on.y 
occur  after  the  vehicle  has  reached  or  exceeded  the  Mach  number, 
relative  to  the  duct,  at  which  the  vehicle  was  designed  to  pass  the 
norsul  shock. 

X£  there  Is  sn  appreciable  region  of  disturbed  flow  between  tb 
shock  and  the  vehicle,  a  great  deal  of  power  will  be  required  to 
push  this  slug  of  high  speed  air  through  the  tube.  The  first  patt 
of  this  report  deals  with  the  determination  of  the  distance  at  whi( 
the  shock  forms  ahead  at  the  vehicle  and  of  its  standoff  distance 
during  the  transition  time  before  it  is  overtaken  by  the  vehicle. 

Tbtt  second  part  of  this  report  deal!  with  the  problem  of  slow: 
down  the  shock  end  of  ihoffcanlhi  distance  and  time  required  to 


overtake  it.  There  are  several  ways  of  doing  this,  but  the  only 
technique  considered  here  is  that  of  the  "shock 'trap"  which  consists 
of  a  section  of  tube  with  perforations  through  which  the  overpressiire 
produced  by  the  shock  can  be  relieved. 


SHOCK  FORMATION  AM)  GROWTH 

Formation  of  a  Shock  by  a  Permeable  Piston 

The  technique  used  to  determine  the  time  and  position  of  the 
initial  formation  of  the  shock  is  similar  to  those  used  by  Foa  (Ref. 
2)  and  Rudinger  (Ref.  3).  Consider  a  vehicle  moving  in  an  infinite 
tube  at  a  velocity  such  that  the  flow  around  the  body  is  choked. 

If  the  ratio  of  the  area  of  the  tube  to  the  area  of  the  annular 
space  between  the  vehicle  and  the  tube,  is  fixed,  then  the  Mach 
number  of  the  flow  juat  ah3w]  of  the  vehicle,  relative  to  the 
vehicle, 

^  R 

will  remain  constant  as  the  vehicle  decelerates. 


The  vehicle  is  moving  at  Kyj  and  the  flow  velocity  just  ahead 
of  the  vehicle  is  u^j  (both  relative  to  the  tube)  and  the  Riemann 
Parameter,  Q  on  the  left  running  characteristic  is  constant  ahead 


of  the  vehicle.  At  time  t“o  consider  that  the  vehicle  is  accelerated 

steadily  at  a  rate  b  =  .  Since  the  flow  around  the  vehicle  is 

dt 


ch'>ked,  the  flow  will  also  be  pushed  along  or  accelerated  as  the 
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vehicle  accelerates.  However,  the  flow  will  not  be  accelerated 
at  the  same  rate  as  the  vehicle  but  at  some  slower  rate,  b  = 

oLv 

Now  the  velocity  of  the  flow  just  ahead  of  the  vehicle,  relative 
to  the  vehicle,  is 


o-  W, 


Hence 


=  lAy  -  Q.B 


Since  Mwa  is  constant 


ar-< 


Also,  since 


Chen 


end  Cherefere 


b'  -  b  -  riv« 

0  •  ,4,  . 


r  o>o.yTf', 


0  M  Ju  .  dv« 
dt  dt 


iiL»  .  ^  sL^ 

'  (Lr  I'-i  dt 

end  Che  ecceleretlon  of  the  flow  just  ahead  of  the  vehicle  is 

b 


b'  = 


I  /^VR 


If  we  now  consider  this  acceleration  as  the  acceleration  of  an 

e<iuivalent  piston,  and  esiploy  the  technique  used  by  Foa  (Ref.  2)  foi 


.  ‘-'n-  ^  ^ 
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the  formation  of  a  shock  by  a  uniformly  accelerated  piston,  we 
find  that  the  time  required  for  the  first  formation  of  the  shock 
is 

f  =  ('^ai  ~  '^¥1^  yvT  [5 


During  this  time  both  the  wave  and  the  vehicle  have  been  moving 
down  the  tube.  The  distance  traveled  by  the  first  wave  is 


b' 


(7) 


and  the  distance  traveled  by  the  vehicle  is 

Xv  ^  TAy/t  f-  b-|-  (8) 

Combining  these  with  eq.  (6)  the  distance  at  which  the  shock  first 
appears  ahead  of  the  vehicle  can  be  found  from  eqs.  7  and  8. 


t  U 


J, 

b' 


(9) 


Numerical  Example 

Consider  a  vehicle  designed  to  pass  a  normal  shock  at  Mach 
number  2.0.  With  this  condition  the  flow  velocity  behind  the 
normal  shock  is  577  ft/sec.  If  the  vehicle  is  designed  so  that  it 
can  accelerate  up  to  577  ft/sec.  before  the  annular  passage  around 
it  becomes  choked  then  it  will  be  able  to  travel  at  the  same  speed, 
relative  to  the  flow  behind  the  shock,  as  the  shock  itself.  In  this 
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way  It:  can  swallow  a  normal  shock  at  Mach  number  2.0. 

Also  at  this  point  let  1000  ft /sec.  and  let  tihe 

vehicle  be  accelerated  at  30  ft /sec.  Under  these  conditions  riyfl  • 
.577  and  eq.  (14)  gives  a  value  of  X,  *  12900  ft.  or  2,45  miles.  How¬ 
ever  the  distance  of  the  shock  ahead  of  the  vehicle  X,  -  X^,  is  onl}* 
2970  ft.  or  .56  miles. 

From  this  point  on  the  growth  of  the  shock  can  be  calculated 
by  means  of  one -dimensional  non-steady  characteristic  theory  (See 
Ref.  1).  Since  the  Mach  number  of  the  flow  just  ahead  of  the  vehiu. 
relative  to  the  vehicle,  remains  constant  then 

T7v  -  a.  A,  f  F, 

The  left  running  characteristic  reaches  the  vehicle  from  upstream 
after  passing  through  the  shock.  The  Rienann  paraineter  on  this 
characteristic  is 

Qr  =  ^  A,  -  m  <11 


Adding  th«i«  two  equations  we  have 
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Now,  substituting  this  into  eq.  (17)  we  have 

1  .  T-J- 

Vv  *  2  ^  ^  ^ 

Mv'fl 


=  '-T 


“  ft 


Now,  by  the  use  of  equations  (13)  and;  (14)  the  Riemann  parameter  on 
the  right  running  characteristic. 


which  originates  at  the  vehicle  can  be  determined,  and  the  chiracter- 

istics  net  can  be  constructed  between  the  vehicle  and  the  shock. 

The  results  of  a  calculation  of  this  type  are  shown  on  Fig.  I-l. 

Although  the  vehicle  in  the  case  considered  is  designed  to  swallow  a 

normal  shock  at  *  2.0,  the  shock  is  still  a  distance  ahead  of  the 

vehicle  at  the  time  the  shock  first  reaches  M  “2.0.  As  a  conse-^ 

s 

quence,  the  vehicle  continues  to  accelerate  and  strengthen  the  shock 
and  does  not  overtake  it  until  “  2.53  and  ■  2.A2. 


SHOCK  TRAP 

The  shock  trap  consists  of  a  section  of  the  tube  with  wall  per¬ 
foration  which  connect  the  tube  to  either  the  atmosphere  outsi(de  or 
to  some  large  sealed  chambers.  Before  the  shock  passes  throxigh  the 
perforated  section,  the  tube  is  under  the  same  pressure  as  either  the 
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outside  atnx>sphere  or  the  sealed  chambers.  ^Theu  the  shock  passes 
the  perforated  section,  the  high  pressure  behind  the  shock  causes 
some  of  the  air  in  the  tube  to  pas#?  out  through  the  perforations. 
Hence  the  mass  flow  rate  in  the  tube  behind  the  shock  is  less  down¬ 
stream  of  the  perforations  than  it  is  upstream  and  the  shock  is 
weakened . 

One  possible  advantage  of  having  the  perforations  connected 
to  large  sealed  chambers  instead  of  open  to  the  atmosphere  is  that 
the  air  which  passes  out  of  the  tube  will  be  stored  in  the  chambetc 
at  the  high  pressure  produced  by  the  shock.  After  the  vehicle  pass, 
through  the  perforated  section  the  air  will  be  forced  back  into  the 

tube  and  thereby  provide  additional  thrust  to  the  vehicle. 

o 

For  simplicity  in  analyzing  the  operation  of  the  shock  trap, 
perforations  will  be  assumed  to  exist  at  only  one  station  and  to  be 
arranged  and  sized  to  remove  a  fixed  percentage  B  of  the  flow  behind 
any  normal  shock.  Prior  to  the  arrival  of  the  shock  the  conditions 
in  the  tube  are  uniform  and  the  air  is  at  rest  throughout  the  tube. 
Now  consider  that  a  shock  is  propagating  through  the  tube  at  a  con¬ 
stant  Mach  number  (see  Fig.  1-2). 

After  the  shock  has  passed  the  perforations  the  continuity  and 
npmentym  equations  across  the  perforated  section  are 

(1 


m,  =.  rh^(l-e) 


1 
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F.  -  F, 

- 

(16) 

fcsing 

the  continuity  equation,  the  momentum  equation  reduces  to 

1 

r 

-  B 

(17) 

where 

the  mass  flow  rate 

is 

-  -B 

r  p'a  D 

t)  ^ 

1 


I  the  stream  force  is  p  ■=  lr>^  V.  t  Pu 

=  -Pa  6 

2 


and 


G  - 


D 


I  f  yM 

(/  f  fnj 
n 


h 


I±L 


The  momentum  equation  can  be  reduced  to 


1“ 

p; 


I  h  -® 


/•UTMJ 

and  the  continuity  equation  becomes 


A 


p; 


4 


0  -  e) 


(18) 


(19) 


til  »*■ 


I-IO 


Combining  these  two  equetlcns  gives 


lirtiere 


^3 


r 


y 

I 


■  '^'1 
/ 

( 

.N  » 

D 

G 

M  ( 1  +  ^  n 

1  -h 

Now 

E  * 

A 

2r-i 

•t 

X 

/  rsRV  \i. 

% 

a*-! 

l'+  T'lJ 

(2 


(2 


Likewise 


n  -m  A  1  I  ^9RTi!  y-/_\ 

'-'-•(TTfTifji'  T/1.) 


(2 


Now  since  there  it  tu>  energjr  added  or  removed  from  the  flow, 
T/  *  T;  and 


E. 

0, 


'  +  T^l-/ 


'  +  T 


(2 


/  Since  the  values  of  P.  and  Q.  are  known  equation  (20)  and  (23) 
can  be  solved  graphically  to  yield  the  Mach  number  on  either  Ride  t 
the  perforated  section^  This  solution  is  shown  on  Fig.  1-3*  By 
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using  an  iterative  procedure  to  determine  the  value  of  Q;  acrosa  the 
shock,  the  ahock  attenuation  resulting  from  any  given  flow  diversion 
can  be  determined. 

For  example,  if  a  shock  trap  is  so  located  that  the  shock 

p 

reaches  the  perforated  station  when  M  =  1.43  then  is  approxi- 

8  G?3 

P 

mately  1.24  and  if  10%  of  the  flow  is  removed  then  is  approxi* 

^3 

mately  1^22  and  Mg  immediately  downstream  of  the  perforated  wection 
is  approximately  1.38.  This  shows  that  a  shock  trap  may  be  useful 
in  reducing  the  strength  of  the  shock. 

•p 

For  shock  Mach  nxnnbers  greater  than  about  2.0  the  value  of  ^ 

Qs 

will  exceed  1.5  and  this  analysis  will  not  apply.  A  more  complicated 
analysis  which  includes  sonic  conditions  at  the  perforated  section 
must  be  developed. 
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SECTION  J 


PRELIMINARY  DESIGN  CONSIDERATIONS  FOR  A 
3 -FT  DIAMETER  TUBE  TRANSPORT  PILOT  MODEL 


ABSTRACT 

This  report  deals  with  the  design  of  a  pilot  model  of  the 
tube  transport  system.  The  means  of  propulsion  considered  here 
is  a  jet  engine,  placed  behind  the  vehicle  and  acting  as  the 
driver  of  an  ejector  pump  which  sucks  the  air  from  ahead  of  the 
vehicle  and  transfers  it  to  the  rear.  The  paper  only  covers 
the  selection  of  some  of  the  basic  parameters  involved. 


SYMBOLS 


oc  cross  sectional  area  (sq  ft) 

111  mass  flow  rate  (slugs/sec) 
p  pressure  (Ib/sq  ft) 

T  temperature  (°R) 

(5  density  (slugs/cu  ft) 
a  flow  velocity  (ft/sec) 

CL  speed  of  sound  (ft/sec) 
n  Mach  number 
f  streaia  force  (lb) 
drag  (lb) 

6  thrust  (lb) 

J  diameter  of  tube  (ft)  ^ 

S  entropy  (Btu/°R  lb) 

P  power  (HP) 

Gr 


J-iti 


SUBSCRIPTS 


1  undisturbed  ambient  conditions 

2  conditions  between  vehicle  and  engine 

3  conditions  after  vehicle  and  engine 

4  conditions  after  mixing 

5  final  conditions  after  heat  extraction 

p  primary 

s  secondary 

SUPERSCRIPTS 
o  stagnation  values 
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INTRODUCTION 


The  problem  considered  here  is  that  Qf  selecting  a  power 
plant  for  a  pilot  model  of  the  tube  transport  system.  It  is 
specified  that  the  vehicle  is  to  travel  at  300  mph.  (M  r  0.4) 
in  a  tube  of  3  ft.  diameter. 

The  analytical  model  used  is  (  <t  of  an  entropy  and  drag 
producer  (the  vehicle),  followed  by  a  thrust  generator,  an 
adiabatic  mixing  region  and  a  heat  extraction  region  (Fig.  J-1). 

The  equations  are  simplified  by  changing  the  frame  of 
reference  from  that  fixed  to  the  walls  to  one  moving  with  the 
vehicle.  In  the  latter  frame  of  reference  the  flow  is  substan¬ 
tially  steady  or  quasi  steady.  The  air  undergoes  a  comp.iete 
cycle  from  a  station  some  distance  upstream  of  the  vehicle  to  a 
station  5  eome  distance  do'fnstream  of  it. 

DRAG  FRODUCTION 

The  first  step  is  to  analyze  what  is  happening  between 
sections  1  and  2,  around  the  vehicle.  The  flow  conditions  at 
station  1  are  taken  to  be  T°  a  500°R  and  p^  =  2,120  psf,  hence 
p  =  0.0024  slugs/cu  ft;  u  =  30C  mpn  i  440  fps  (Mi  =  0.4). 
With  o<j=:  dL-.-  -  7.07  sq  ft,the  mass  flow  rate  is:  au,  =  7.6 

slugs/sec» 


The  flow  transformation  from  section  1  to  2  Is  Isoenergetl 
hence:  .  tJ  .  Sie^R, 

The  best  criterion  presently  available  for  the  design  of 

//  " 

the  air  passage  around  the  vehicle  is  to  follow  the  current 
practice  for  jet-engine  air  intakes.  If  this  is  done,  it  is 
permissible  to  assume  a  0.95,  hence  p§  s  2,240  psf. 

'  t 

The  other  flow  parameters  at  station  2  are  obtained  from 
the  conservation  equations  (Ref.  1).  „ 


Continuity  e^n. 

Pa  ^  ^ 

/W  ‘  fW 

(1> 

or 

F,h/^  _  %  Afe. 

ffi  '  fa 

Momentum  eqn.: 

^ =  /i* “t  -Pa  “'!* 

)! 

(3) 

Energy  eqn.: 

1 

T»  T®  - 
<»  "  - 

(4) 

0 

^f> 

One  easily  finds  from  the  continuity  eqn;  M2  =  0.426,  whll 
the  aKMMntun  eqn.  fives  the  drag  ^  ■  F2  -  Pi  «  >745  lb. 

The  entropy  production  is  given  by; 

AS=^C.l,^.1\SnJit-S3.S^o.SS=2.73 

A  /a 

THRL%ST  GENERATOR 

Since  the  three  conservation  eqns.  have  five  unknowns: 
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P3,  T3,  M2,  0  ,  ,  two  additional  relationships  are  required. 
One  of  those  requirements  is  that  the  thrust  B  must  exactly  bal¬ 
ance  the  drag  ,  or  B  --2)  1=.  745  lb. 


Now  the  eqns.  can  be  solved  for  T2,  P3,  ^2  with  ^  as 
parameter.  However,  the  values  for  ^  are  limited  in  both 
directions.  Due  to  the  second  law  of  thermodynamics,  the  total 
entropy  change  (i.e.  -d>S’  of  the  flow  plus  the  surroundings)  must 
be  positive  for  a  spontaneous  process.  Hence  a  minimum  value  for 
is  established.  On  the  other  hand,  from  the  continuity  eqn. 


^  (zjy  follows,  that,  for  the  subsonic  case  and  in 


absence  of  shocks,  is  limited.  Hence  from  the  energy  eqiu  (4) 
can  be  found. 

/a  />iax 


An  approximate  expression  for  the  power  required  is  found 
in  Ref.  2,  nl. : 

'P-  (5) 

where  p^  s  stagnation  pressure  ratio. 

Hence  P  «  2,120  x  7.07  x  440  x  (1  -  0.95)  x  1.818  x  10“^^  « 
600  HP. 

With  these  two  parameters,  nl.  the  thrust  B  «  745  lb  and 
the  power  P  s  600  HP,  the  two  determining  characteristics  for  the 

power  plant  are  fixed.  Furthermore  the  mass  flow  rate  through 


the  engine  should  be  less  than  =  7.68  slugs/sec  while  its 
cross  sectional  area  is  also  limited  (such  that  always  M  <  1). 

The  thrust  of  the  power  plant  and  also  its  cycle  efficiency 
can  be  increased  by  using  aii  inlet  diffuser  ahead  of  the  jet 
engine. 

The  available  catalogues  give  several  possible  power  plants, 
which  satisfy  the  above  stated  requirements  of  thrust,  power, 
diameter  and  mass  flow  rate  (see  e.g.,  Ref.  3). 

The  following  numerical  example  gives  an  idea  of  the  ord^r 
of  magnitude  of  the  different  parameters  involved.  The  jet 
engine  has  a  constant  cross  sectional  area:  oLp  =  1.67  sq  ft. 

The  parameters  at  stations  1  and  2  are  already  known. 

At  station  2,  the  flow  is  divided  into  two  parts: 

1.  the  primary  flow  (p)  through  the  engine  under¬ 
goes  a  Rayleigh  process  with  heat  addition. 

^  2.  the  secondary  flow  (s)  around  the  engine  under¬ 

goes  an  isoenSrgetic  process  with  a  drop  in 
stagnation  pressure. 

A  uniform  flow  is  assumed  at  the  inlet  over  the  whole  croas  sec- 

0 

tional  area,  hence  =  — rfip  a  1.815  slugs/sec  and 
a  5.865  slugs /sec. 


0 


J-5 


For  the  secot.dary  flow  process,  with  T^g  =  T2  -  and 

^3s  *  0,96  P2  =  2.150  psf,  the  conservation  eqns.  yield  s 
0.450  and  Rg  =  432  lb.  This  additional  drag  must  also  be 
balanced  by  the  thrust  of  the  primary  flow.  Hence  the  total 
thrust  is:  =  745  t  432  =  1,177  lb. 

For  the  primary  flow  process,  the  amount  of  heat  to  be  added 
is: 

a  ~  .Ji— r  — — -  Btu/lb  =  7.27  Btu/lb 

^  777.:5 

The  conservation  eqns.  give  further:  T^p  =  546.3*^R; 

F3p  =5,320  lb;  M3p  =  0.312  and  p3p  =  3,000  lb. 

The  flow  conditions  at  station  4  are  determined  by  applying 
the  equations  for  isoenergetic  mixing  in  a  constant  area  duct. 


From  Ref.  1  one  gets; 

/  - 

L 

hence  z  s  0.6825. 


Furthermore:  Ml  z:  — LZ-^^  yields  M,  r  0.403 

=  18>300  Ih,  gives  -  2,110  lb 
u  z:  — H - ft/sec 


m 


J~6 

=  4-/'  /  -  ~Lsl 

and  =  506.6  °R 

\ 

Reference.  4  gives  now;  p°  =  2^360  psf  and  =  523.3  °R.  \ 

// 

In  the  next  part  of  the  duct,  between  stations  4  and  5, 
heat  is  extracted  in  such  an  amount  that  the  net  result  betweer 
2  and  5  is  zero  energy  addition.  Hence:  f-,  ~  ^ 

or  Btu/lb. 

Furthermore  from  the  conservations  eqns.  T5  =  515.1  “^R 
«  0.399  and  P5  »  2^360  psf.  Except  for  some  minor  calcula¬ 
tion  errors,  the  final  state  is  the  same  as  the  initial  one. 
Hence  it  is  possible  for  the  flow  to  undergo  the  above  trans¬ 
formations  and  to  come  back  to  the  same  original  state.  However 
nothing  Is  said  about  the  relative  distances  involved.  Probably 
a  compression  wave  travels  ahead  of  the  vehicle^  so  that  the  un¬ 
disturbed  section  1  Is  far  upstream.  Maybe  the  compression  wav«t 
Is  weakened  by  an  expansion  caused  by  the  sucking  of  air  around 
the  vehicle  through  the  engine.  This  problem  is  still  under 
further  study.  On  the  other  hand,  the  length  of  the  mixing  regi( 
and  energy  extraction  region  can  be  very  large  too,  A  more  de¬ 
tailed  atudy  of  the  nechaniam  of  energy  decay  la  found  in  Ref.  5t 
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CONCLUSION 

S  •  I  / 

From  the  vehicle  under  consideration,  nwving  at  a  speed  of 
300  mph  (M  =  0.4)  in  a  duct  of  3  ft,  diameter,  a  jet  engine  is 
needed  which  can  produce,  in  steady  motion,  a  thrust  of  750  lb 
corresponding  to  a  power  of  600  HP.  However  during  accelera¬ 
tion,  the  required  values  of  9  and  P  are  higher  so  that  the 
capacity  of  the  power  plant  must  be  higher  too,  depending  upon 
the  desired  acceleration.  The  diameter  of  the  engine  must  be 
less  than  a  certain  critical  value  (to  keep  M  1). 

About  the  support  for  the  vehicle,  some  basic  calculations 
can  be  found  in  Ref.  2. 

e 

For  the  calculation  of  the  power  required,  the  following 
expression  from  Ref.  2  was  used;  P  =  f^) 

P  is  a  linear  function  of  u^,  hence  of  but  is  iiidependent 
of  the  weight  of  the  vehicle.  However  the  power  calculated  with 
eqn.  (5)  is  the  power  requited  for  .steady  operation  only,  and  is 
determined  rather  by  the  shape  of  the  vehicle  (determines  drag 
and  p°)  than  by  its  weight.  The  power  required  ftr  the  accelera¬ 
tion,  ho<7ever,  is  proportional  to  the  weight  of  the  vehicle. 
Furthermore,  for  an  aerodynamic  suspension  system,  the  additional 
pov?er  for  the  air  support  would  increase  with  the  weight.  Though 
this  p-Twer  is  very  low  because  of  "ground  effect"  it  has  been 
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shown  in  Ref.  6  that  most  of  this  energy  caii  be  recovered  for 
thrust  generation.  A  separate  power  source  in  the  nose  of  the 
vehicle  could  deliver  this  airstream  for  the  "aerodynamic 

y 

cushions"  (see  Ref.  2). 
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RECOMMENDATIONS 


A  radically  new  riiode  of  transportation  cannot  be  properly 
evaluated  in  the  light  of  today's  problems  and  needs;  rather  it 
must  be  judged  on  the  basis  of  its  adequacy  to  serve  the  demands 
of  tomorrow,  and  of  its  probable  influence  on  the  development 
and  use  of  the  land  and  on  the  patterns  of  living  for  generations 
to  come. 

A  tube-transport  system  of  the  type  discussed  in  the  preced¬ 
ing  sections  would  provide  a  mode  of  movement  characterized  by 
ylery  high  speed,  a  high  level  of  Inherent  safety,  and  complete 
independence  from  atmospheric  and  surface  conditions.  This  mode 
could  directly  penetrate  the  large  urban  centers  and  also  provide 
direct  access  to  large  regional  airports. 

The  immediate  benefit  would  be  the  provision  of  very  high 
speed,  all-f^eather,  direct  mass  transit  between  distant  cities  and 
within  large  metropolitan  areas.  Its  greatest  eventual  value,  how¬ 
ever,  will  in  all  likelihood  be  the  impetus  provided  to  urban  growth 
and  the  freedom  permitted  to  the  patterns  of  growth  as  the  res’trtc- 
tio^  of  our  present  forms  of  transportation  are  removed. 

It  is  believed  that  this  system  would  be  best  suited  for 
operation  at  speeds  from  250  to  2000  over  distances  where  the 


vise  of  such  speeds  would  be  economically  justifiable. 

The  first  step  in  its  implementation  should  probably  be 
a  pilot,  mail -carrying  line;  a  second  step  could  be  a  shuttle 
service  for  passengers  between  t:;o  cities  not  too  far  apart 
(e.g.,  New  York  and  Philadelphis  or  Washington)  or  between  one 
of  these  metropolitan  centers  and  its  relocated  **super  airport". 
The  operating  frequency  of  these  shuttle  services  could  be  very 
high,  because  of  the  intrinsic  safety  of  this  transportation 
system. 

o 

It  should  be  noted  that  the  engineering  know-how  for  the 
innilemenfcation  of  these  first  steps  (covering  speeds  up  to  300 
or  350  a^h)  is  already  available.  Thus,  from  the  standpoint  of 
technical  feasibility,  it  is  believed  that  these  steps  could  be 
undertaken  Inmediately.  The  utilization  of  the  system  over  longer 
distances  and  at  higher  speeds  would  have  to  be  a  gradual  develop¬ 
ment,  following  the  successful  solution  of  the  engineering  problen 
associated  with  the  transition  to  supersonic  speeds  and  with  the 

if 

propulsion,  support  and  stability  of  the  vehicle  at  scqrersonlc 
speeds.  An  effort  in  this  "direction  is  believed  to  be  highly  desi 
able  from  the  very  start,  no  matter  how  successful  the  "low-speed' 

f 

,  versions  may  turn  out  to  be,  because  all  indications  of  the  analys 
are  that  the  superiority  of  the  proposed  carrier  over  existing 

•  b 
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transport  systems  would  be  far  greater  at  supersonic  than  at 
subsonic  speeds. 

The  successful  development  of  a  transportation  system  based 
on  the  ideas  discussed  in  this  report  will  undoubtedly  take 
several  years.  Unfortunately,  the  novelty  of  the  concept  makes 
it  impossible  to  predict  all  the  problems  that  will  arise  in  this 
work  and  to  plan  the  developmental  program  very  far  in  advance. 
Recommendations  for  further  work  will,  therefore,  be  confined  to 
a  program  for  the  first  two  years. 

This  program  should  consist  of: 

(a)  engineering  analyses  and  research  aimed  at  the 
determination  of  quantitative  information  required 
for  a  study  of  the  economic  feasibility  of  the 
proposed  system; 

(b)  design  of  a  three-mile  pilot  model,  and 

(c)  a  preliminary  feasibility  study  of  the  proposed 
system. 

The  engineering  analysis  should  deal  with  the  following  areas: 

(1)  Propulsion.  Various  energy  sources  and  thrust -genera ting 
cycles  should  be  considered  and  analyzed,  including  -  for 
very  high  speeds  -  thermal  propulsion  systems  of  the  type 
diagrammatically  illustrated  in  Fig.  6. 

(2)  Transonic  acceleration.  Attention  should  be  given  to  non- 
steady,  flow  means  for  overcoming  the  choking  difficulties. 


and  to  the  uae  of  suitable  wall  slots  and  perforations! 
over  those  portions  of  the  tube  where  the  transonic 
acceleration  is  expected  to  take  place; 

(3)  Vehicular  support  and  stability.  The  study  of  stability 
and  ground-effect  support  should  be  extended  to  higher 
speeds  and  larger  models,  with  special  attention  to  new 
means  for  increasing  the  stability  derivatives  of  the 
proposed  system,  thereby  permitting  larger  vehicular  loads 
and  higher  speeds  in  f  '^rn.  Improvements  in  curtain  momentum 
recovery  for  thrust  should  also  be  sought,  through  a  detail¬ 
ed  study  of  the  flow  around  the  ground  support  pads.  Finally, 
alternative  systems  of  vehicular  support  should  be  explored. 

The  desj,gn  of  the  pilot  model  will  also  involve  a  considerabi 
amount  of  research.  The  problems  of  elastic  stability  and  structu 
al  dynamics  of  high-speed  tubed  vehicles  are  believed  to  be  siraila 
in  many  respects,  to  those  of  high-speed  aircraft.  Therefore,  the 
tools  for  the  design  end  analysis  of  vehicular  structures  are  al¬ 
ready  available  to  a  great  extent  in  aeroelasticity  and  in  aircraf 
structural  dynamics.  In  addition,  there  are  special  structural 
problems  that  are  associated  with  the  peculiar  pressure  distribute 
in  the  annular  transfer  space  around  the  vel  icle.  Other  design  pj 
lems  that  should  receive  special  attention  are  those  of  cabin  air 


conditioning  and  structural  cooling  at  the  higher  speeds,  of  air 
purging  in  the  tube,  and  of  automatic  control.  The  design  of  the 
tube  itself,  its  methods  of  support,  and  limitations  in  aligtspent 
should  also  be  studied  during  this  phase  of  the  project. 

The  economic  analysis  conducted  concurrently  with 

the  other  phases  of  the  proposed  resddreh.  The  economic  feasi¬ 
bility  study  must  of  necessity  be  of  a  preliminary  nature  until 
all  design  factors  and  operational  characteristics  are  resolved. 
Consideration  should  be  given,  however,  to  the  probable  public  ac¬ 
ceptance  and  use,  to  the  effects  on  cc!n5>lementing  and  con5)eting 
inodes  of  transportation,  and  to  eatinataa  of  capital  and  operating 


coats. 


